Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Probability 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 55 records · Page 3

All-electron molecular tunnel ionization based on the weak-field asymptotic theory in the integral representation

Tunnel ionization (TI) underlies many important ultrafast processes, such as high-harmonic generation andstrong-field ionization. Among the existing theories for TI, many-electron weak-field asymptotic theory (ME-WFAT) is by design capable of accurately treating many-electron effects in TI. An earlier version of ME-WFATrelied on an accurate representation of the asymptotic tail of the orbitals, which hindered its implementation inGaussian-basis-set-based quantum chemistry programs. In this work, we reformulate ME-WFAT in the integralrepresentation, which makes the quality of the asymptotic tail much less critical, hence greatly facilitating itsimplementation in standard quantum chemistry packages. The integral reformulation introduced here is thereforemuch more robust when applied to molecules with arbitrary geometry. Here, we present several case studies, amongwhich is the CO molecule where some earlier theories disagree with experiments. Here we find that ME-WFATproduces the largest ionization probability when the field points from C to O, as experiments suggest. Anattractive feature of ME-WFAT is that it can be used with various types of multielectron methods whether ofdensity functional or multiconfiguration types, this inturn facilitates tunnel ionization calculation in systems exhibiting a strong multireference character.

74 ATOMIC AND MOLECULAR PHYSICS

Transient anisotropic kernel for probabilistic learning on manifolds

PLoM (Probabilistic Learning on Manifolds) is a method introduced in 2016 for handling small training datasets by projecting an Itô equation from a stochastic dissipative Hamiltonian dynamical system, acting as the MCMC generator, for which the KDE-estimated probability measure with the training dataset is the invariant measure. PLoM performs a projection on a reduced-order vector basis related to the training dataset, using the diffusion maps (DMAPS) basis constructed with a time-independent isotropic kernel. In this paper, we propose a new ISDE projection vector basis built from a transient anisotropic kernel, providing an alternative to the DMAPS basis to improve statistical surrogates for stochastic manifolds with heterogeneous data. The construction ensures that for times near the initial time, the DMAPS basis coincides with the transient basis. For larger times, the differences between the two bases are characterized by the angle of their spanned vector subspaces. The optimal instant yielding the optimal transient basis is determined using an estimation of mutual information from Information Theory, which is normalized by the entropy estimation to account for the effects of the number of realizations used in the estimations. Consequently, this new vector basis better represents statistical dependencies in the learned probability measure for any dimension. Three applications with varying levels of statistical complexity and data heterogeneity validate the proposed theory, showing that the transient anisotropic kernel improves the learned probability measure.

Diffusion maps

Elastic cross section is entanglement entropy

We present universal relations between entanglement entropy, which quantifies the quantum correlation between subsystems, and the cross section, which is the primary observable for high-energy particle scattering, by employing a careful formulation of wave packets for the incoming particles. For 2-to-2 elastic scattering with no initial entanglement and subdividing the system along particle labels, we show that both the Rényi and Tsallis entropies in the final states are directly proportional to the elastic cross section in units of the transverse size for the initial wave packets, which is then interpreted as the elastic scattering probability. The relations do not depend on the underlying dynamics of the quantum field theory and are valid to all orders in coupling strengths. Furthermore, computing quantum correlations between momentum and nonkinematic data leads to entanglement entropies expressed as various semi-inclusive elastic cross sections. Our result gives rise to a novel “area law” for entanglement entropy in a two-body system. Published by the American Physical Society 2025

Low, Ian (ORCID:0000000275709597)

Pair wave function symmetry in UTe2 from zero-energy surface state visualization

Although nodal spin-triplet topological superconductivity appears probable in uranium ditelluride (UTe 2 ), its superconductive order parameter Δ k remains unestablished. In theory, a distinctive identifier would be the existence of a superconductive topological surface band, which could facilitate zero-energy Andreev tunneling to an s-wave superconductor and also distinguish a chiral from a nonchiral Δ k through enhanced s-wave proximity. In this study, we used s-wave superconductive scan tips and detected intense zero-energy Andreev conductance at the UTe 2 (0-11) termination surface. Imaging revealed subgap quasiparticle scattering interference signatures with a-axis orientation. The observed zero-energy Andreev peak splitting with enhanced s-wave proximity signifies that Δ k of UTe 2 is a nonchiral state: B 1u , B 2u , or B 3u . However, if the quasiparticle scattering along the a axis is internodal, then a nonchiral B 3u state is the most consistent for UTe 2 .

Gu, Qiangqiang [Cornell Univ., Ithaca, NY (United

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Bayesian inference of nuclear-matter density from proton scattering

Background: Proton elastic scattering at intermediate energy is widely employed as a tool for determining the matter radius of atomic nuclei. Here, the sensitivity of the approach relies on high-resolution measurements at small scattering angles and low-momentum transfer. Under these conditions, the Glauber multiple scattering theory accurately describes the proton-nucleus elastic cross section. Purpose: Investigate the sensitivity of the Glauber multiple scattering theory to uncertainties associated with input parameters such as the nuclear-matter density distribution and nucleon-nucleon data. Method: A joint Bayesian inference was performed using 12 angular distributions of elastic scattering at different energies on 58 Ni, 90 Zr, and 208 Pb targets. A Metropolis-Hastings algorithm was implemented to make an uncertainty quantification analysis for the input parameters used in the Glauber multiple scattering theory. Results: The experimental cross sections were fitted simultaneously using a joint Bayesian inference approach. Posterior probability density distributions of 42 input parameters were obtained from the analysis. A moderate correlation between the nuclear density parameters and the nucleon-nucleon cross sections was found. This correlation impacts the extraction of the nuclear-matter radius. Conclusions: The present analysis provided a consistent method for extracting the nuclear-matter density distribution of 58 Ni, 90 Zr, and 208 Pb from data across different incident energies. Due to the correlation of the nucleon-nucleon cross sections with the other input parameters, a constrained Bayesian inference using free nucleon-nucleon cross section data was performed. The nuclear-matter radii obtained from the analysis are in good agreement with multiple results reported in the literature.

190 ≤ A ≤ 219

Development of Parameters for the Particle Size Distribution of TATB

Laser light scattering (LLS), manual counting of scanning electron microscopy (SEM) images, and reduction of SEM images using ImageJ (open source program) to determine Feret (caliper) diameters were applied to determine particle size distribution (PSD) of four preparations of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB) and yttria‐stabilized zirconia (YSZ), an SEM certified standard. Mie theory was used to reduce the LLS data. The spherical nature of the YSZ made it a good candidate for LLS. Variations in n , the refractive index, and iκ, the imaginary component, produced very little change in the PSD. However, changing the carrier liquid from H 2 O to a 40% aqueous sucrose solution, thereby changing the carrier refractive index, n 0 , substantially affected the PSD. The Mie complex refractive indices for the YSZ were n = 2.200, iκ = 0.100, with a 40% aqueous sucrose solution, n 0 = 1.400. The triclinic crystal structure of TATB made refractive index determinations more difficult, so a study was conducted varying Mie parameters and comparing them to the same data reduced using the Fraunhofer theory. Changing the n and κ parameters produced PSD with a small concentration of particles less than 1 µm in size or none in this range. SEM images, Feret data, manual counting, and Fraunhofer data reduction indicate particles less than 1 µm are probably < 5% in concentration. The final selection of Mie parameters for TATB was n = 2.283, iκ = 0.1, and suspension medium, n 0 = 1.330. Finally, computations, using density functional theory produced similar parameters.

Feret diameter

Nonperturbative quantum gravity in a closed Lorentzian universe

We study how meaningful physical predictions can arise in nonperturbative quantum gravity in a closed Lorentzian universe. In such settings, recent developments suggest that the quantum gravitational Hilbert space is one-dimensional and real for each α-sector, as induced by spacetime wormholes. This appears to obstruct the conventional quantum-mechanical prescription of assigning probabilities via projection onto a basis of states. While previous approaches have introduced external observers or augmented the theory to resolve this issue, we argue that quantum gravity itself contains all the necessary ingredients to make physical predictions. We demonstrate that the emergence of classical observables and probabilistic outcomes can be understood as a consequence of partial observability: physical observers access only a subsystem of the universe. Tracing out the inaccessible degrees of freedom yields reduced density matrices that encode classical information, with uncertainties exponentially suppressed by the environment’s entropy. We develop this perspective using both the Lorentzian path integral and operator formalisms and support it with a simple microscopic model. Our results show that quantum gravity in a closed universe naturally gives rise to meaningful, robust predictions without recourse to external constructs.

AdS-CFT Correspondence

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations

Control-Affine Schrödinger Bridge and Generalized Bohm Potential

From a stochastic control perspective, the Schrödinger bridge is a density-valued continuous curve parameterized by time that connects a given pair of initial and terminal probability densities via minimum effort controlled Brownian motion. The control-affine Schrödinger bridge extends this idea to a generic control-affine Itô diffusion, possibly with an additive state cost. Here, in this letter, we recast the necessary conditions of optimality for the control-affine Schrödinger bridge problem as a two point boundary value problem for a quantum mechanical Schrödinger PDE with complex potential. This complex-valued potential is a generalization of the real-valued Bohm potential in quantum mechanics. Our derived potential is akin to the optical potential in nuclear physics where the real part of the potential encodes elastic scattering (transmission of wave function), and the imaginary part encodes inelastic scattering (absorption of wave function). The key takeaway is that the process noise that drives the evolution of probability densities induces an absorbing medium in the evolution of wave function. These results make new connections between control theory and non-equilibrium statistical mechanics through the lens of quantum mechanics.

Markov processes

Noncollinear ground states of solids with a source-free exchange correlation functional

In this paper, we expand upon the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and coworkers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+U+J functional. In the process, we compare the noncollinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground-state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin-polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analyses on probability currents, net torque, and constant terms draw additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and pave the way for future extensions and corrections to the SF corrected XC functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

On symmetry-resolved generalized entropies

Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Analytic Sensitivity Coefficients for Bethe's Solution of the Neutron Slowing Down Equation

Neutron slowing down theory is used to derive expressions for the sensitivity coefficients of the neutron collision density in a hydrogenous infinite medium with respect to the fixed source, scattering probability, and macroscopic nuclear cross sections. Analytic expressions for Bethe’s solution of the neutron slowing down equation are derived for the constant cross section approximation with a point, uniform, and gamma lethargy spectrum. Analytic expressions for the corresponding sensitivity coefficients are derived and used to verify Monte Carlo neutron transport calculations.

Analytic Benchmark

The impact of nickel concentration and stacking fault energy on deformation mechanisms in high-purity austenitic Fe-Cr-Ni alloys

Understanding how composition affects deformation mechanisms in austenitic stainless steels is essential for developing accurate predictive models of stress-induced failures and stress corrosion cracking. Nickel (Ni), an element classified as a critical element, plays a crucial role in these processes. It is important to examine how Ni concentration influences stacking fault energy (SFE) and, consequently, the deformation mechanisms of austenitic stainless steels. However, in commercial stainless steels, the effects of other alloying elements and impurities can obscure Ni's role, complicating efforts to isolate its impact. Here, in this study, we use two high-purity Fe-Cr-Ni alloys to investigate how Ni concentration and SFE interact to alter deformation mechanisms and induce martensitic transformation. By combining in situ synchrotron X-ray diffraction (XRD) tensile testing and post-mortem electron microscopy with density functional theory simulations, we gain precise insights into these phenomena. We find that the Fe18Cr10Ni (wt%) alloy, with its low SFE, exhibits higher stacking fault probability, deformation-induced martensitic transformation, and a lesser increase in dislocation density with plastic strain. In contrast, the Fe18Cr14Ni (wt%) alloy, with its higher SFE, shows enhanced deformation twinning and greater dislocation density with increasing strain. These findings from high-purity ternary alloys provide valuable insights that can guide the search for alternative elements to replace Ni while achieving similar effects on phase stability and deformation behavior.

36 MATERIALS SCIENCE

Efficient Unitary Designs from Random Sums and Permutations

A unitary k-design is an ensemble of unitaries that matches the first k moments of the Haar measure. In this work, we provide two efficient constructions of k-designs on n-qubits using new random matrix theory techniques. Our first construction is based on exponentiating sums of random i.i.d. Hermitian matrices and uses O(k2n2)-many gates. In the spirit of central limit theorems, we show that this random sum approximates the Gaussian Unitary Ensemble (GUE). We then show that the product of just two exponentiated GUE matrices is already approximately Haar random. Our second construction is based on products of exponentiated sums of random permutations and uses Õ(k poly (n)) many gates. The k dependence is optimal (up to polylogarithmic factors) and is inherited from the efficiency of existing k-wise independent permutations. Furthermore, replacing random permutations with quantum-secure pseudorandom permutations (PRPs), we also obtain a pseudorandom unitary (PRU) ensemble that is secure under nonadaptive queries. A central feature of both proofs is a new connection between the polynomial method in quantum query complexity and the large-dimension (N) expansion in random matrix theory. In particular, the first construction uses the polynomial method to control high moments of certain random matrix ensembles without requiring delicate Weingarten calculations. In doing so, we define and solve a moment problem on the unit circle, asking whether a finite number of equally weighted points can reproduce a given set of moments. In our second construction, the key step is to exhibit an orthonormal basis for irreducible representations of the partition algebra that has a low-degree large-N expansion. This allows us to show that the distinguishing probability is a low-degree rational polynomial of the dimension N.

algebra

Charge radii measurements of exotic tin isotopes in the proximity of $N=50$ and $N=82$

We report nuclear charge radii for the isotopes $^{104-134}$Sn, measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the $N=50$ and $N=82$ shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models fall short, however, of reproducing the magnitude of the known $B(E2)$ transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

FOS: Physical sciences

Four-Dimensional Quasielastic Neutron Scattering (4D-QENS) Analysis of Correlated Ionic Conduction in Sr⁢Cl2

Methods of elucidating the mechanisms of fast-ion conduction in solid-state materials are pivotal for advancements in energy technologies such as batteries, fuel cells, sensors, and supercapacitors. In this study, we examine the ionic conduction pathways in single crystal Sr⁢Cl2, which is a fast-ion conductor above 900 K, using four-dimensional quasielastic neutron scattering (4D-QENS). We explore both coherent and incoherent neutron scattering at temperatures above the transition temperature into the superionic phase to explore the correlated motion of hopping anions. Refinements of the incoherent QENS yield residence times and jump probabilities between lattice sites in good agreement with previous studies, confirming that ionic hopping along nearest-neighbor directions is the most probable conduction pathway. However, the coherent QENS reveals evidence of de Gennes narrowing, indicating the importance of ionic correlations in the conduction mechanism. This highlights the need for improvements both in the theory of ionic transport in fluorite compounds and the modeling of coherent 4D-QENS in single crystals.

Coles, Jared [Materials Science Division, Argonne

Operationally induced preferred basis in unitary quantum mechanics

The preferred-basis problem and the definite-outcome aspect of the measurement problem persist even if the detector is modeled unitarily, because experimental data are necessarily represented in a Boolean event algebra of mutually exclusive records whereas the theoretical description is naturally formulated in a noncommutative operator algebra with continuous unitary symmetry. This change of mathematical type constitutes the core of the 'cut': a structurally necessary interface from group-based kinematics to set-based counting. In the presented view the basis relevant for recorded outcomes is not determined by the system Hamiltonian alone; it is induced by the measurement mapping, i.e., by the detector channel together with the coarse-grained readout that defines an instrument. The probabilistic mapping is anchored in symmetry and measure theory: by Gleason-type uniqueness (Gleason for projections in $d>2$ and Busch's extension for Positive Operator-Valued Measures (POVMs) including $d=2$), the trace rule is the unique probability measure consistent with additivity over exclusive events and basis-independence of the unitary sector. A compact qubit--pointer model yields an induced unsharp POVM $E_\pm=\tfrac12(\id\pm η\,σ_z)$ with $η$ fixed by pointer resolution, displaying explicitly how the detector induces the relevant basis. Finally, nested-observer paradoxes are tightened into a non-composability lemma: joint assignment of outcome propositions is obstructed unless a joint instrument exists. This relocates the origin of randomness to the stochasticity of the transition rules.

Pronskikh, Vitaly [Fermilab] (ORCID:00000002518174