Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlocal models”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonlocal nucleon matrix elements in the rest frame

Extracting parton structure from lattice quantum chromodynamics (QCD) calculations requires studying the coordinate scale 𝑧 3 dependence of the matrix elements of bilocal operators. The most significant contribution comes from the 𝑧 3 dependence induced by ultraviolet (UV) renormalization of the Wilson line. We demonstrate that the next-to-leading order perturbative calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice QCD rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. The residual discrepancies can be modeled by a leading effect from the structure of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Practical methods for simulating photo-induced dynamics of molecules often rely on a trajectory-description of the nuclear dynamics, incorporating approximate electron-nuclear feedback. Through approximating electron–nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. Here, a set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Analysis of Neural Networks as Random Dynamical Systems

In this report we present our findings and outcomes of the NNRDS (analysis of Neural Networks as Random Dynamical Systems) project. The work is largely motivated by the analogy of a large class of neural networks (NNs) with a discretized ordinary differential equation (ODE) schemes. Namely, residual NNs, or ResNets, can be viewed as a discretization of neural ODEs (NODEs) where the NN depth plays the role of the time evolution. We employ several legacy tools from ODE theory, such as stiffness, nonlocality, autonomicity, to enable regularization of ResNets thus improving their generalization capabilities. Furthermore, armed with NN analysis tools borrowed from the ODE theory, we are able to efficiently augment NN predictions with uncertainty overcoming wellknown dimensionality challenges and adding a degree of trust towards NN predictions. Finally, we have developed a Python library QUiNN (Quantification of Uncertainties in Neural Networks) that incorporates improved-architecture ResNets, besides classical feed-forward NNs, and contains wrappers to PyTorch NN models enabling several major classes of uncertainty quantification methods for NNs. Besides synthetic problems, we demonstrate the methods on datasets from climate modeling and materials science.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum order by disorder in frustrated spin nanotubes

Here we investigate quantum order by disorder in a frustrated spin nanotube formed by wrapping a J1–J2 Heisenberg model at 45° around a cylinder. Using Schwinger boson theory and density matrix renormalization group (DMRG), we have computed the ground-state phase diagram to reveal a Z2 phase in which collinear spin stripes form a right- or left-handed helix around the nanotube. We have derived an analytic estimate for the critical η c = J 1 /2J 2 of the Z 2 -helical phase transition, which is in agreement with the DMRG results. By evaluating the entanglement spectrum and nonlocal string order parameters we discuss the topology of the Z 2 -helical phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Engineering 3D Floquet Codes by Rewinding

Floquet codes are a novel class of quantum error-correcting codes with dynamically generated logical qubits arising from a periodic schedule of noncommuting measurements. We utilize the interpretation of measurements in terms of condensation of topological excitations and the rewinding of measurement sequences to engineer new examples of Floquet codes. In particular, rewinding is advantageous for obtaining a desired set of instantaneous stabilizer groups on both toric and planar layouts. Our first example is a Floquet code with instantaneous stabilizer codes that have the same topological order as the three-dimensional (3D) toric code(s). This Floquet code also exhibits a splitting of the topological order of the 3D toric code under the associated sequence of measurements, i.e., an instantaneous stabilizer group of a single copy of the 3D toric code in one round transforms into an instantaneous stabilizer group of two copies of the 3D toric code up to nonlocal stabilizers in the following round. We further construct boundaries for this 3D code and argue that stacking it with two copies of the 3D subsystem toric code allows for a transversal implementation of the logical non-Clifford controlled-controlled- Z gate. We also show that the coupled-layer construction of the X-cube Floquet code can be modified by a rewinding schedule such that each of the instantaneous stabilizer codes is finite depth equivalent to the X-cube model up to toric codes; the X-cube Floquet code exhibits a splitting of the X-cube model into a copy of the X-cube model and toric codes under the measurement sequence. Our final 3D example is a generalization of the 2D Floquet toric code on the honeycomb lattice to three dimensions, which has instantaneous stabilizer codes with the same topological order as the 3D fermionic toric code. Published by the American Physical Society2024

Physics↗
DOE OSTI2024

Giant spin transport anisotropy in magnetic topological insulators

Here, we report on exceptionally long spin transport and giant spin lifetime anisotropy in the gapped surface states of three-dimensional (3D) magnetic topological insulators (MTIs). We examine the properties of these states using the three-dimensional Fu-Kane-Mele Hamiltonian model in the presence of a magnetic exchange field. The corresponding spin textures of surface states, which are well reproduced by an effective two-band model, hint at a considerable enhancement of the lifetime of out-of-plane spins compared to in-plane spins. This is confirmed by large-scale spin transport simulations for 3D MTIs with disorder. The energy dependence of the spin lifetime anisotropy arises directly from the nontrivial spin texture of the surface states, and is correlated with the onset of the quantum anomalous Hall phase. Our findings suggest different spin filtering capabilities of the gapped surface MTI states, which could be explored by nonlocal spin valve measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Algebraic Hastatic Order in One-Dimensional Two-Channel Kondo Lattice

The two-channel Kondo lattice likely hosts a rich array of phases, including hastatic order, a channel symmetry breaking heavy Fermi liquid. In this work, we revisit its one-dimensional phase diagram using density matrix renormalization group and, in contrast to previous work, find algebraic hastatic orders generically for stronger couplings. These are heavy Tomonaga-Luttinger liquids with nonanalyticities at Fermi vectors captured by hastatic density waves. We also find a predicted additional nonlocal order parameter due to interference between hastatic spinors, not present at large N, and residual repulsive interactions at strong coupling suggesting non-Fermi-liquid physics in higher dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Flat Bands and Chiral Optical Response of Moiré Insulators

Here, we present a low-energy model describing the reconstruction of the electronic spectrum in twisted bilayers of honeycomb crystals with broken sublattice symmetry. The resulting moiré patterns are classified into two families with different symmetry. In both cases, flat bands appear at relatively large angles without any magic-angle condition. Transitions between them give rise to sharp resonances in the optical absorption spectrum at frequencies well below the gap of the monolayer. Owing to their chiral symmetry, twisted bilayers display circular dichroism, that is, different absorption of left and right circularly polarized light. This optical activity is a nonlocal property determined by the stacking. In hexagonal boron nitride, sensitivity to the stacking leads to strikingly different circular dichroism in the two types of moirés. Our calculations exemplify how subtle properties of the electronic wave functions, encoded in current correlations between the layers, control physical observables of moiré materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Low-Density Parity-Check Codes as Stable Phases of Quantum Matter

Phases of matter with robust ground-state degeneracy, such as the quantum toric code, are known to be capable of robust quantum information storage. Here, we address the converse question: given a quantum error-correcting code, when does it define a stable gapped quantum phase of matter, whose ground-state degeneracy is robust against perturbations in the thermodynamic limit? We prove that a low-density parity-check (LDPC) code defines such a phase, robust against all few-body perturbations, if its code distance grows at least logarithmically in the number of degrees of freedom, and it exhibits “check soundness.” Many constant-rate quantum LDPC expander codes have such properties, and define stable phases of matter with a constant zero-temperature entropy density, violating the third law of thermodynamics. Our results also show that quantum toric-code phases are robust to spatially nonlocal few-body perturbations. Similarly, phases of matter defined by classical codes are stable against symmetric perturbations. In the classical setting, we present improved locality bounds on the quasiadiabatic evolution operator between two nearby states in the same code phase.

quantum error correction↗
DOE OSTI2020

Controllable vortex shedding from dissipative exchange flows in ferromagnetic channels

Ferromagnetic channels subject to spin injection at one extremum sustain long-range coherent textures that carry spin currents known as dissipative exchange flows (DEFs). In the weak injection regime, spin currents carried by DEFs decay algebraically and extend through the length of the channel, a regime known as spin superfluidity. Similar to fluids, these structures are prone to phase slips that manifest as vortex-antivortex pairs. In this work, we numerically study vortex shedding from DEFs excited in a magnetic nanowire with a physical obstacle. Using micromagnetic simulations, we find regimes of laminar flow and vortex shedding as a function of obstacle position tunable by the spin injection sign and magnitude. Vortex-antivortex pairs translate forward (VF regime) or backward (VB regime) with respect to the detector's extremum, resulting in well-defined spectral features. Qualitatively similar results are obtained when temperature, anisotropy, and weak nonlocal dipole fields are included in the simulations. These results provide clear features associated with DEFs that may be detected experimentally in devices with nominally identical boundary conditions. Furthermore, our results suggest that obstacles can be considered as DEF control gates, opening an avenue to manipulate DEFs via physical defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Many-Body Scars as a Group Invariant Sector of Hilbert Space

We present a class of Hamiltonians H for which a sector of the Hilbert space invariant under a Lie group G, which is not a symmetry of H, possesses the essential properties of many-body scar states. These include the absence of thermalization and the “revivals” of special initial states in time evolution. Here, a particular class of examples concerns interacting spin-1/2 fermions on a lattice consisting of N sites (it includes deformations of the Fermi-Hubbard model as special cases), and we show that it contains two families of N + 1 scar states. One of these families, which was found in recent literature, comprises the well-known η-pairing states. We find another family of scar states that is U(N) invariant. Both families and most of the group-invariant scar states produced by our construction, in general, give rise to the off-diagonal long range order, which survives at high temperatures and is insensitive to the details of the dynamics. Such states could be used for reliable quantum information processing because the information is stored nonlocally and thus cannot be easily erased by local perturbations. In contrast, other scar states we find are product states, which could be easily prepared experimentally. The dimension of scar subspace is directly controlled by the choice of group G and can be made exponentially large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree–Fock exchange. Furthermore, these results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Comment on “Nonlocal nucleon matrix elements in the rest frame”

In a recent paper, “Nonlocal nucleon matrix elements in the rest frame” [1], it was observed that the next-to-leading order calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice quantum chromodynamics (QCD) rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. We argue that perturbative QCD breaks down at such a distance scale after resumming the associated large logarithms, while the Ansatz used in the analysis there is not justified in perturbation theory. Besides, we explain the observation in Ref. [1] and demonstrate that the Ansatz fails to describe the data for 𝑧 > 0.3 fm, showing an opposite trend. Finally, although Ref. [1] proposes multiplying the Ansatz by a Gaussian correction model, which is shown to reduce the discrepancy with the data, this does not legitimize the use of perturbative QCD at such distance scales.

Lattice QCD↗
DOE OSTI2022

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

Abstract We study the production of very light elements ( Z < 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10 −6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10 −6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r -process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗