Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Small-scale signatures of primordial non-Gaussianity in k-nearest neighbour cumulative distribution functions

ABSTRACT Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbour cumulative distribution functions, kNN-CDFs, to the quijote-png simulations to explore the sensitivity of kNN-CDFs to primordial non-Gaussianity. An interesting result is that for halo samples with $M_\mathrm{ h}\langle 10^{14}$ M$_\odot$ $h^{-1}$, the kNN-CDFs respond to equilateral PNG in a manner distinct from the other parameters. This persists in the galaxy catalogues in redshift space and can be differentiated from the impact of galaxy modelling, at least within the halo occupation distribution (HOD) framework considered here. kNN-CDFs are related to counts-in-cells and, through mapping a subset of the kNN-CDF measurements into the count-in-cells picture, we show that our results can be modelled analytically. A caveat of the analysis is that we only consider the HOD framework, including assembly bias. It will be interesting to validate these results with other techniques for modelling the galaxy–halo connection, e.g. (hybrid) effective field theory or semi-analytical methods.

Coulton, William R. (ORCID:0000000212973673)↗
DOE OSTI2023

Measurement of Magnetic Field using Doppler-Free Saturation Spectroscopy (DFSS) in C-2W FRC Plasma

The Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and TAE Technologies (TAE) focused on implementing an existing mobile Doppler-free saturation spectroscopy (DFSS) diagnostic to determine the direction of the equilibrium magnetic field within TAE’s C-2W field-reversed configuration (FRC) fusion device. Non-perturbative measurements of magnetic fields in the core of FRCs can provide critical information needed to constrain equilibrium models and for the continued advancement of the fusion research program at TAE. DFSS is a well-established, laser-based absorption spectroscopy technique used to suppress inhomogeneous (Doppler) broadening and reveal the Zeeman fine structure of the target emission lines. The collaboration involved the use of ORNL’s recently constructed mobile DFSS diagnostic. The DFSS diagnostic was shipped and installed on C- 2W in February 2023. TAE’s expertise was used to design and construct the mechanic support structures required to mount the optical hardware on C-2W. Additionally TAE provide critical input to minimize the noise imposed on the DFSS signal by vibrations of the mechanic support structure during operations. DFSS data was obtained for over 500 C-2W shots and will continue until the expiration of this CRADA on September 30 th 2023. Over the course of this work, significant progress was made to maximize the signal to noise (SNR) ratio. Analysis of recently acquired data is on-going, however, initial indications are highly suggestive that a DFSS signal above background has been successfully acquired in C-2W.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microstrata

Abstract Microstrata are the non-extremal analogues of superstrata: they are smooth, non-extremal (non-BPS) solitonic solutions to IIB supergravity whose deep-throat limits approximate black holes. Using perturbation theory and numerical methods, we construct families of solutions using a consistent truncation to three-dimensional supergravity. The most general families presented here involve two continuous parameters, or amplitudes, and four quantized parameters that set the angular momenta and energy levels. Our solutions are asymptotic to the vacuum of the D1-D5 system: AdS 3 ×S 3 × 𝕋 4 . Using holography, we show that the they are dual to multi-particle states in the D1-D5 CFT involving a large number of mutually non-BPS supergravitons and we determine the anomalous dimensions of these states from the binding energies in supergravity. These binding energies are uniformly negative and depend non-linearly on the amplitudes of the states. In one family of solutions, smoothness restricts some of the fields to lie on a special locus of the parameter space. Using precision holography we show that this special locus can be identified with the multi-particle states constructed via the standard OPE of the single-particle constituents. Our numerical analysis shows that microstrata are robust at large amplitudes and the solutions can be obtained to very high precision.

Physics↗
DOE OSTI2020

Atmospheric teleconnection processes linking winter air stagnation and haze extremes in China with regional Arctic sea ice decline

Abstract. Recent studies suggested significant impacts of boreal cryosphere changes on wintertime air stagnation and haze pollution extremes in China. However, the underlying mechanisms of such a teleconnection relationship remains unclear. Here we use the Whole Atmosphere Community Climate Model (WACCM) to investigate dynamic processes leading to atmospheric circulation and air stagnation responses to Arctic sea ice changes. We conduct four climate sensitivity experiments by perturbing sea ice concentrations (SIC) and corresponding sea surface temperature (SST) in autumn and early winter over the whole Arctic and three subregions in the climate model. The results indicate distinct responses in circulation patterns and regional ventilation to the region-specific Arctic changes, with the largest increase of both the probability (by 132 %) and the intensity (by 30 %) of monthly air stagnation extremes being found in the experiment driven by SIC and SST changes over the Pacific sector of the Arctic (the East Siberian and Chukchi seas). The increased air stagnation extremes are mainly driven by an amplified planetary-scale atmospheric teleconnection pattern that resembles the negative phase of the Eurasian (EU) pattern. Dynamical diagnostics suggest that convergence of transient eddy forcing in the vicinity of Scandinavia in winter is largely responsible for the amplification of the teleconnection pattern. Transient eddy vorticity fluxes dominate the transient eddy forcing and produce a barotropic anticyclonic anomaly near Scandinavia and wave train propagation across Eurasia to the downstream regions in East Asia. The piecewise potential vorticity inversion analysis reveals that this long-range atmospheric teleconnection of Arctic origin takes place primarily via the middle and upper troposphere. The anomalous ridge over East Asia in the middle and upper troposphere worsens regional ventilation conditions by weakening monsoon northwesterlies and enhancing temperature inversions near the surface, leading to more and stronger air stagnation and pollution extremes over eastern China in winter. Ensemble projections based on state-of-the-art climate models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) corroborate this teleconnection relationship between high-latitude environmental changes and midlatitude weather extremes, though the tendency and magnitude vary considerably among each participating model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗
DOE OSTI2020

Improved Frequency-Domain Steady-State Modeling of the Dual-Active-Bridge Converter Considering Finite ZVS Transition Time Effects

This article presents an improved analytical modeling (IAM) approach for the dual-active-bridge (DAB) converter in the frequency-domain. Specifically, finite transition times are incorporated into the improved modeling that are significant in DAB converters operated at high frequency. The IAM approach is first motivated by comparing traditional frequency-domain modeling against realistic simulation results in LTspice. Next, the IAM approach is mathematically developed in an iterative fashion, where the conventional frequency-domain approach is used as the initialization of the operational behavior of the converter. Subsequent iterations of the IAM algorithm incorporate finite rise- and falltimes of device output capacitances based on detailed switching behavior analyses, which may considerably perturb the desired modulation parameters (i.e., phase shift and zero states). While the initially proposed IAM approach demonstrates advantages over ideal modeling, it is still shown to exhibit errors when applied to device transitions with significant nonlinearity in parasitic output capacitance. To account for this, a three-slope approximation to the switching transition is then motivated and developed, which enables significantly improved steady-state prediction accuracy. Simulation and experimental results confirm unacceptable inaccuracies associated with the conventional frequency-domain modeling across a wide range of operating conditions, alongside enhanced accuracy provided by the proposed modeling approach.

14 SOLAR ENERGY↗
DOE OSTI2022

A Search for Neutron to Mirror Neutron Oscillation Using Neutron Electric Dipole Moment Measurements

Baryon number violation is a key ingredient of baryogenesis. It has been hypothesized that there could also be a parity-conjugated copy of the standard model particles, called mirror particles. The existence of such a mirror universe has specific testable implications, especially in the domain of neutral particle oscillation, viz. the baryon number violating neutron to mirror-neutron (n−n′) oscillation. Consequently, there were many experiments that have searched for n−n′ oscillation, and imposed constraints upon the parameters that describe it. Recently, further analysis on some of these results have identified anomalies which could point to the detection of n−n′ oscillation. All the previous efforts searched for n−n′ oscillation by comparing the relative number of ultracold neutrons that survive after a period of storage for one or both of the two cases: (i) comparison of zero applied magnetic field to a non-zero applied magnetic field, and (ii) comparison where the orientation of the applied magnetic field was reversed. However, n−n′ oscillations also lead to variations in the precession frequency of polarized neutrons upon flipping the direction of the applied magnetic field. Precession frequencies are measured, very precisely, by experiments searching for the electric dipole moment. For the first time, we used the data from the latest search for the neutron electric dipole moment to constrain n−n′ oscillation. After compensating for the systematic effects that affect the ratio of precession frequencies of ultracold neutrons and cohabiting 199Hg-atoms, chief among which was due to their motion in non-uniform magnetic field, we constrained any further perturbations due to n−n′ oscillation. We thereby provide a lower limit on the n−n′ oscillation time constant of τnn′/|cos(β)|>5.7s,0.36T′

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Enhancing nuclear data validation analysis by using machine learning

In this work, we demonstrate how machine learning (ML) techniques can be used as an auxiliary tool for nuclear data validation analysis. The output of the ML analysis can inform evaluators and validators on the quality, or lack thereof, of specific nuclear data and benchmark experiments with respect to simulating these experimental benchmark values. To this end, measured and simulated effective neutron multiplication factors, k eff , of 875 selected ICSBEP critical assemblies and the latter's sensitivities with respect to nuclear data as well as benchmarks' features (e.g., material nuclide, core geometry) are used as input for a random forest (RF) regression model. The RF is used to encode the complex inter-dependencies between thousands of nuclear data features (i.e., sensitivity profiles and aspects describing the measurements) and 875 simulated and experimental k eff values in order to understand what nuclear data features are most informative for predicting bias. The complexity of relationships and high-dimensional space is difficult-to-impossible to search fully using simply expert judgment. As a first proof-of-concept—step, it is demonstrated that this technique is able to correctly trace large discrepancies between simulated and experimental k eff back to fabricated shortcomings in nuclear data that were perturbed to simulated k eff values. In a second, real-case scenario, step, the RF algorithm is used to validate the ENDF/B-VIII.0 library in comparison to ENDF/B-VII.1 nuclear data. One case is showcased where the chosen ML algorithms highlighted nuclear data (the 19 F(n,inl) cross section from 0.4–0.9 MeV) that are shown to be problematic by comparing them to associated differential experimental data and nuclear data from other libraries. In addition to that, it is shown that the RF results point towards poor benchmark experiments and associated underestimated uncertainties (e.g., the PU-SOL-THERM-028 series). However, using the RF algorithm for validating nuclear data with respect to k eff is currently limited to pinpointing groups of questionable nuclear data due to the inherent correlations between features introduced by the nuclear data themselves and how k eff is simulated. Due to this, we recommend that the ML methods presented be used to augment—rather than replace—the expert knowledge of evaluators and validators.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2024

Self-consistent quasiparticle G W and hybrid functional calculations for Al/InAs/Al heterojunctions: Band offset and spin-orbit coupling effects

The electronic structure of surfaces and interfaces plays a key role in the properties of quantum devices. Here, we study the electronic structure of realistic Al/InAs/Al heterojunctions using a combination of density functional theory with hybrid functionals and state-of-the-art quasiparticle G W (QSGW) calculations. We find a good agreement between QSGW calculations and hybrid functional calculations, which themselves compare favorably well with angle-resolved photoemission spectroscopy experiments. Our paper confirms the need for well-controlled quality of the interfaces to obtain the needed properties of InAs/Al heterojunctions. A detailed analysis of the effects of spin-orbit coupling on the spin splitting of the electronic states shows a linear scaling in k space, related to the two-dimensional nature of some interface states. The good agreement by QSGW and hybrid functional calculations opens the door towards trustable use of an effective approximation to QSGW for studying very large heterojunctions. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of injected flux and current temporal phasing on self-organization in the HIT-SI3 experiment

The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗
DOE OSTI2021

Gravitational effective field theory islands, low-spin dominance, and the four-graviton amplitude

Abstract We analyze constraints from perturbative unitarity and crossing on the leading contributions of higher-dimension operators to the four-graviton amplitude in four spacetime dimensions, including constraints that follow from distinct helicity configurations. We focus on the leading-order effect due to exchange by massive degrees of freedom which makes the amplitudes of interest infrared finite. In particular, we place a bound on the coefficient of the R 3 operator that corrects the graviton three-point amplitude in terms of the R 4 coefficient. To test the constraints we obtain nontrivial effective field-theory data by computing and taking the large-mass expansion of the one-loop minimally-coupled four-graviton amplitude with massive particles up to spin 2 circulating in the loop. Remarkably, we observe that the leading EFT coefficients obtained from both string and one-loop field-theory amplitudes lie in small islands. The shape and location of the islands can be derived from the dispersive representation for the Wilson coefficients using crossing and assuming that the lowest-spin spectral densities are the largest. Our analysis suggests that the Wilson coefficients of weakly-coupled gravitational physical theories are much more constrained than indicated by bounds arising from dispersive considerations of 2 → 2 scattering. The one-loop four-graviton amplitudes used to obtain the EFT data are computed using modern amplitude methods, including generalized unitarity, supersymmetric decompositions and the double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Discovery of a Damped Lyα Absorber Originating in a Spectacular Interacting Dwarf Galaxy Pair at z = 0.026

We present the discovery of neutral gas detected in both damped Lyα absorption (DLA) and H I 21 cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an H I bridge connecting two interacting dwarf galaxies (log (M star /M ⊙ ) = 8.5 ± 0.2) that host a z = 0.026 DLA with log[N(H I )/cm –2 ] = 20.60 ± 0.05 toward the QSO J2339–5523 (z QSO = 1.35). At impact parameters of d = 6 and 33 kpc, the dwarf galaxies have no companions more luminous than ≈0.05L * within at least Δv = ±300 km s –1 and d ≈ 350 kpc. The H I 21 cm emission is spatially coincident with the DLA at the 2σ–3σ level per spectral channel over several adjacent beams. However, H I 21 cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with a spin temperature to covering factor ratio T s /f c > 1880 K). Observations with VLT-MUSE demonstrate that the α-element abundance of the ionized interstellar medium (ISM) is consistent with the DLA (≈10% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf–dwarf interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and H I 21 cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗