Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion theory”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗
DOE OSTI2025

Infalling ultra-faint dwarfs as emissaries of the Axiverse

Recent discoveries of ultra-faint dwarf galaxies (UFDs) infalling onto the Milky Way, namely Leo K & M at r ≃ 450 kpc, considerably strengthens the case that UFDs constitute a distinct galaxy class that is inherently smaller and fainter, and metal-poorer than the classical dwarf spheroidals (dSph). This distinction is at odds with the inherent continuity of galaxy halo masses formed under scale-free gravity for any standard dark-matter (DM) model. Here, we show that distinct galaxy classes do evolve in cosmological simulations of multiple light bosons representing the “Axiverse” proposal of string theory, where a discrete mass spectrum of axions is generically predicted to span many decades in mass. In this context, the observed UFD class we show corresponds to a relatively heavy boson of 3 × 10 −21 eV, including Leo K & M, whereas a lighter axion of 10 −22 eV comprises the bulk of DM in all larger galaxies including the dSphs. Although Leo M is larger in size than Leo K, we predict its velocity dispersion to be smaller ( ≃ 1.7 km/s) than that of Leo K ( ≃ 4.5 km/s) because of the inverse de Broglie scale dependence on momentum. This scenario can be definitively tested using millisecond pulsars close to the Galactic center, where the Compton frequencies of the heavy and light bosons imprint monotone timing residuals that may be detected by the Square Kilometre Array (SKA) on timescales of approximately one week and four months, respectively.

dark matter↗
DOE OSTI2024

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quadrupole collectivity in the neutron-rich sulfur isotopes 38,40,42,44 S

Electromagnetic transition strengths in the even-even neutron-rich sulfur isotopes 38,40,42,44 S were measured using intermediate-energy Coulomb excitation at the National Superconducting Cyclotron Laboratory. By utilizing the sensitivity of the experimental technique to E2 excitations from the ground state, the evolution of the pattern of B(E2) strengths to several low-lying 2 + states was investigated at Z = 16 from near stability to the N = 28 island of inversion. The experimental results allowed a detailed comparison with predictions from shell-model calculations using the SDPF-MU Hamiltonian, which was designed to describe collectivity in this region of the nuclear chart. While the shell-model calculations succeeded in modeling transition strengths at N = 22, 24, 26, the experimental B(E2; $0$ $^{+}_{1}$ → $2$ $^{+}_{1}$) at N = 28 was smaller than the predicted value by about a factor of 2, similar to previous observations for chlorine and argon isotopes around N = 28. The dependence of this overprediction by theory on the choice of effective charges was explored.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
DOE OSTI2021

First-principles theory for Schottky barrier physics

We develop a first-principles theory for Schottky barrier physics. The Poisson equation is solved selfconsistently with the electrostatic charge density over the entire barrier using the density functional theory (DFT) electronic structure converged locally, allowing computation of a Schottky barrier entirely from DFT involving thousands of atomic layers in the semiconductor (SC). Here, the induced charge in the bulk consists of conduction and valence band charges from doping and band bending, as well as charge from the evanescent states in the gap of the SC. The Schottky barrier height (SBH) is determined when the induced charge density and the induced electrostatic potential reach self-consistency. Tests on the GaAs-graphene and Si/Al heterostructures yield SBH, width, along with depletion and inversion layers obtained self-consistently as functions of temperature and bulk doping.

36 MATERIALS SCIENCE↗
DOE OSTI2022

α-cluster structure of 18 Ne

In this work, we investigate α-clustering in 18 Ne and compare it with what is known about clustering in the mirror nucleus 18 O. The excitation function for resonant elastic scattering of 14 O + α was measured in inverse kinematics using the TexAT active-target detector system, covering excitation energies from 8 to 17 MeV. The data were analyzed using a multi-channel R-matrix approach. Detailed spectroscopic information—including excitation energies, spin and parity assignments, and partial α and total widths—was extracted from the R-matrix analysis and compared with theoretical models and previously available data. Correspondence between levels in 18 O and 18 Ne is established. An extensive shell-model analysis of the 18 O– 18 Ne mirror pair shows very good agreement between theory and experiment, especially in the distribution of clustering strength. The comparison reveals that certain states, particularly at higher excitation energies, exhibit significantly stronger α-cluster character than predicted. This suggests that these states have a collective structure aligned with the corresponding α-reaction channel.

6 ≤ A ≤ 19↗
DOE OSTI2023

Non-perturbative production of fermionic dark matter from fast preheating

We investigate non-perturbative production of fermionic dark matter in the early universe. We study analytically the gravitational production mechanism accompanied by the coupling of fermions to the background inflaton field. The latter leads to the variation of effective fermion mass during preheating and makes the resulting spectrum and abundance sensitive to its parameters. Assuming fast preheating that completes in less than the inflationary Hubble time and no oscillations of the inflaton field after inflation, we find an abundant production of particles with energies ranging from the inflationary Hubble rate to the inverse duration of preheating. The produced fermions can account for all observed dark matter in a broad range of parameters. As an application of our analysis, we study non-perturbative production of fermionic dark matter in the model of Palatini Higgs inflation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Deep learning of interface structures from simulated 4D STEM data: cation intermixing vs. roughening ∗

Abstract Interface structures in complex oxides remain an active area of condensed matter physics research, largely enabled by recent advances in scanning transmission electron microscopy (STEM). Yet the nature of the STEM contrast in which the structure is projected along the given direction precludes separation of possible structural models. Here, we utilize deep convolutional neural networks (DCNN) trained on simulated 4D STEM datasets to predict structural descriptors of interfaces. We focus on the widely studied interface between LaAlO 3 and SrTiO 3 , using dynamical diffraction theory and leveraging high performance computing to simulate thousands of possible 4D STEM datasets to train the DCNN to learn properties of the underlying structures on which the simulations are based. We test the DCNN on simulated data and show that it is possible (with >95% accuracy) to identify a physically rough from a chemically diffuse interface and create a DCNN regression model to predict step positions. We quantify the applicability of the model to different thicknesses and the transferability of the approach. The method shown here is general and can be applied for any inverse imaging problem where forward models are present.

42 ENGINEERING↗
DOE OSTI2021

Observation of spatially resolved Rashba states on the surface of CH 3 NH 3 PbBr 3 single crystals

Hybrid organic-inorganic perovskites (HOIPs) are prime candidates for studying Rashba effects due to the heavy metal and halogen atoms in their crystal structure coupled with predicted inversion symmetry breaking. Nevertheless, observation of the Rashba effect in cubic CH 3 NH 3 PbBr 3 single crystals that possess bulk inversion symmetry is the subject of extensive debate due to the lack of conclusive experiments and theoretical explanations. In this work, we provide experimental evidence that Rashba state in cubic CH 3 NH 3 PbBr 3 single crystals at room temperature occurs exclusively on the crystal surface and depends on specific surface termination that results in local symmetry breaking. We demonstrate this using a suite of spatially resolved and depth-sensitive techniques, including circular photogalvanic effect, inverse spin Hall effect, and multiphoton microscopy, that are supported by first principle calculations. Our work suggests using surface Rashba states in these materials for spintronic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Isospin 0 and 2 two-pion scattering at physical pion mass using all-to-all propagators with periodic boundary conditions in lattice QCD

A study of two-pion scattering for the isospin channels, 𝐼 = 0 and 𝐼 = 2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of 𝐿 = 4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the 𝐼 = 0 and 𝐼 = 2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for 𝐼 = 0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying 𝐾 → 𝜋⁢𝜋 decay, on physical-mass ensembles using the GEVP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Pd-promoted reduction and restructuring of an In 2 O 3 -based catalyst for CO 2 hydrogenation at room temperature

An unconventional reaction mechanism in an In 2 O 3 /Pd(1 1 1) inverse model catalyst for the CO 2 hydrogenation reaction has been uncovered: In 2 O 3 is partially reduced at room temperature in a reaction atmosphere as a result of its direct contact with Pd(1 1 1), which is an efficient H 2 splitter. The reduction induces changes in surface free energy, leading to a dynamical restructuring at the In 2 O 3 /Pd(1 1 1) interface via formation of InO x and outward diffusion of Pd, as revealed by ambient pressure X-ray photoelectron spectroscopy, X-ray absorption spectroscopy and density functional theory simulations. This dynamical restructuring eventually promotes the growth of 2D InPd y O x nanodomains as the catalytically active phase and the exclusive formation of methanol upon hydrogenation of CO 2 at room temperature. A comparable high selectivity toward CH 3 OH was found in more realistic bulk catalytic systems (2 wt% Pd/In 2 O 3 catalyst and commercial CZA catalyst). Scanning tunneling microscopy under ultrahigh vacuum and ambient pressure reaction atmospheres further reveals the structural dynamics at the InO x /Pd(1 1 1) interface, where we follow in situ the evolution of the InO x particles on Pd(1 1 1) and the mobility of the InPd y O x nanodomains in a CO 2 + H 2 environment. The present findings of the formation of a mixed oxide phase in a dynamically restructuring metal/reducible-oxide interface indicate further implications for other heterogeneous catalytic systems beyond the present CO 2 hydrogenation example and highlight the importance of in situ investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantum fluctuations in dense plasma simulations

Molecular dynamics (MD) simulations are a powerful tool for modeling warm and hot dense matter. Density functional theory (DFT) MD simulations are often preferred in dense plasmas in order to accurately model quantum electronic structure. However, DFT-MD simulations neglect interaction effects due to fluctuations in excited states. In this work, we present an MD approach that uses excited state method pseudoatoms to run dense plasma simulations with many different core-electron configurations at classical MD speeds. We also allow for transitions between different configurations in our simulations and find that these fluctuations are especially important for highly excited states. Our results suggest that finite configuration lifetimes that are comparable to the inverse ion plasma frequency need to be accounted for in order to accurately model ion distributions in dense plasma simulations. We also demonstrate that excited state fluctuations have a direct impact on ion plasma microfields, generate different plasma microfields for different excitation levels, and thereby induce absorption–emission line shape asymmetries even in steady-state plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Lorentzian inversion formula and the spectrum of the 3d O(2) CFT

We study the spectrum and OPE coefficients of the three-dimensional critical O(2) model, using four-point functions of the leading scalars with charges 0, 1, and 2 ( s , $\phi$, and t ). We obtain numerical predictions for low-twist OPE data in several charge sectors using the extremal functional method. We compare the results to analytical estimates using the Lorentzian inversion formula and a small amount of numerical input. We find agreement between the analytic and numerical predictions. We also give evidence that certain scalar operators lie on double-twist Regge trajectories and obtain estimates for the leading Regge intercepts of the O(2) model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗