Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “operator 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 91 records · Page 5

Symmetry-specific orientational order parameters for complex structures

A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linear stability analysis via simulated annealing and accelerated relaxation

Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectral localizer for line-gapped non-Hermitian systems

Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Para-particle oscillator simulations on a trapped-ion quantum computer

Deformed oscillators allow for a generalization of the standard fermions and bosons, namely, for the description of para-particles. Such particles remain hypothetical and unobserved in nature; yet, they can model physical phenomena, such as topological phases of matter. Here, we report the digital quantum simulation of para-particle oscillators by mapping para-particle states to the state of a qubit register, which allows us to identify the para-particle oscillator Hamiltonian as an XY model and further digitize the system onto a universal set of gates. In both instances, the gate depth grows polynomially with the number of qubits used. To establish the validity of our results, we experimentally simulate the dynamics of para-fermions and para-bosons, demonstrating full control of para-particle oscillators on a quantum computer. Furthermore, we compare the overall performance of the digital simulation of dynamics of the driven para-Fermi oscillator to a recent analog quantum simulation result.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixed quantum-classical methods for polaron spectral functions

In this work, using two distinct semiclassical approaches—namely, the mean-field Ehrenfest method and the mapping approach to surface hopping—we investigate the spectral function of a single charge interacting with phonons on a lattice. This quantity is relevant for the description of angle-resolved photoemission experiments. Focusing on the one-dimensional Holstein model, we compare the performance of these approaches across a range of coupling strengths and lattice sizes, exposing the relative strengths and weaknesses of each. We demonstrate that these approaches can be efficiently applied with reasonable accuracy to ab initio polaron models. Furthermore, our work provides a route to the calculation of spectral properties in realistic electron–phonon-coupled systems in a computationally inexpensive manner with encouraging accuracy.

Atomic and molecular spectra↗

A simple fourth order propagator based on the Magnus expansion in the Liouville space: Application to a Λ-system and assessment of the rotating wave approximation

A simple fourth-order propagator [Ture and Jang, J. Phys. Chem. A 128, 2871 (2024)] based on the Magnus expansion is extended to the Liouville space for both closed-system and Lindbladian open-system quantum dynamics. For both dynamics, commutator free versions of fourth-order propagators are provided as well. These propagators are then applied to the dynamics of a driven Λ-system, where Lindblad terms represent the effect of a photonic bath. For both dynamics, the accuracy of the rotating wave approximation (RWA) for the matter–radiation interaction is assessed. We confirmed reasonable performance of RWA for weak and resonant fields. However, small errors appear for moderate fields and substantial errors can be found for strong fields where coherent population trapping can still be expected. We also found that the presence of bath for open-system quantum dynamics consistently reduces the errors of the RWA. These results provide quantitative information on how the RWA breaks down beyond weak field or for non-resonant cases. Major results are benchmarked against results of our sixth-order ME-based propagator. Finally, we also provide numerical comparison of our algorithms with other fourth-order algorithms for the Λ-system. These confirm reasonable performance of our simple propagators and the improvement gained through commutator-free expressions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Critical Review of Heat Pipe Experiments in Nuclear Energy Applications

This critical review provides heat pipe (HP) experimental data sets that contain pertinent information regarding nuclear technology that may be beneficial to researchers. Heat pipes have been shown to have a tremendously positive impact on nuclear technologies and will continue to become a more prevalent technology as more nuclear reactor concepts embrace this robust technology. Most previous reviews may focus on only a specific HP design or application, and some are backdated. Herein, this critical review extends previous efforts; integrates and summarizes previously reported HP experimental efforts; and provides updates with recently reported results in the literature for HPs in all nuclear-related applications, including space power (thermal radiators, core cooling, and electricity production), microreactors (emergency core cooling, hybrid control rods, and reactor core cooling), and HP involvement in other nuclear-related technologies (spent fuel pool cooling). The two main objectives of this critical review are (1) to facilitate the development of HP codes by outlining some of the existing experimental data sets to validate their codes and directing developers to these efforts and (2) to provide comprehensive information regarding the vast applicability of HPs used in the nuclear industry, including the theory of operation and limitations to supplement researchers in the development of new ideas for potential applications in nuclear-related technologies. The review clearly shows extensive and diverse experimental data sets for HPs developed under diverse testing conditions depending on the available nuclear application for validation purposes. Thus, this critical review is oriented to providing attention to the existing efforts rather than determining gaps in HP research.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lepton flavor violation: From muon decays to muon colliders

We investigate the unique potential of a high-energy muon collider to probe lepton-flavor-violating signals arising from physics beyond the Standard Model (SM). Low-energy, precision searches for charged lepton flavor violation (LFV) are projected to dramatically improve their sensitivity in the coming years and could provide the first evidence of new physics. We interpret the sensitivity of these searches in terms of a set of LFV operators in the SM effective field theory. The same operators are then probed at the TeV scale via new, high-energy processes only available at a high-energy muon collider, such as 𝜇⁢𝜇 → 𝜇⁢𝜏 or the scattering of a muon of an electroweak gauge boson into LFV final states. We find that, for most operators, a muon collider could confirm signals if they are seen at future low-energy experiments, whereas for certain flavor combinations it extends the reach to scales well beyond those accessible at lower energies. We also project the sensitivity of a muon collider to lepton-flavor-violating decays of the SM Higgs boson and demonstrate improved sensitivity to ℎ → 𝑒⁢𝜏 and ℎ → 𝜇⁢𝜏 by an order of magnitude compared to the High-Luminosity LHC. The importance of having multiple, complementary probes is illustrated by considering both various combinations of operators and relative sizes of flavor-violating transitions between generations under various assumptions for the flavor structure of new physics.

Asadi, Pouya [University of California, Santa Cruz↗

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗

Thermal WIMPs and the scale of new physics: global fits of Dirac dark matter effective field theories

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework GAMBIT. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from Planck, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

State-operator correspondence in celestial conformal field theory

The bulk-to-boundary dictionary for 4D celestial holography is given a new entry defining 2D boundary states living on oriented circles on the celestial sphere. The states are constructed using the 2D CFT state-operator correspondence from operator insertions corresponding to either incoming or outgoing particles which cross the celestial sphere inside the circle. The BPZ construction is applied to give an inner product on such states whose associated bulk adjoints are shown to involve a shadow transform. Scattering amplitudes are then given by BPZ inner products between states living on the same circle but with opposite orientations. 2D boundary states are found to encode the same information as their 4D bulk counterparts, but organized in a radically different manner.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Second-order spectral line shift comparisons

The second-order spectral line width formulae from the projection operator and kinetic theory methods were recently compared. It was shown that a systematic expansion of the projection operator width expression including initial correlations formally agrees with the second-order kinetic theory result. It is now shown that the second-order dynamic shifts are also formally the same. The static shifts, however, differ due to an ad hoc treatment of electron-electron correlations in the projection operator method. The approximation is necessary in order to screen the radiator-electron interactions. The differences, however, are expected to be small. Finally, the results suggest using the rigorous and more compact second-order width and shift expressions from the kinetic theory method as the starting point for spectral line shape calculations. At line center, however, the projection operator second-order expression for the width and shift simplifies and reduces to the kinetic theory result.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of second-order spectral line widths formulae

Frequently used second-order spectral line width formulae from the projector operator and kinetic theory methods have not been formally compared previously. Here it is shown that a systematic second-order expansion of the projection operator expression including initial correlations agrees with the second-order kinetic theory result. The agreement assumes a common approximation in the projector operator method that introduces a screened radiator–perturber interaction to account for neglected electron–electron correlations. On the other hand, it is shown that the usual width expression from the projection operator approach neglecting initial correlations differs from kinetic theory. The differences, however, are at least third order in the radiator–perturber interaction. The comparisons suggest using the more compact width expression from kinetic theory, which includes initial correlations and a systematic screening of the radiator–perturber interactions, as the starting point for second-order width calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗

Relational bulk reconstruction from modular flow

Abstract The entanglement wedge reconstruction paradigm in AdS/CFT states that for a bulk qudit within the entanglement wedge of a boundary subregion$$ \overline{A} $$ A ¯ , operators acting on the bulk qudit can be reconstructed as CFT operators on$$ \overline{A} $$ A ¯ . This naturally fits within the framework of quantum error correction, with the CFT states containing the bulk qudit forming a code protected against the erasure of the boundary subregionA. In this paper, we set up and study a framework for relational bulk reconstruction in holography: given two code subspaces both protected against erasure of the boundary regionA, the goal is to relate the operator reconstructions between the two spaces. To accomplish this, we assume that the two code subspaces are smoothly connected by a one-parameter family of codes all protected against the erasure ofA, and that the maximally-entangled states on these codes are all full-rank. We argue that such code subspaces can naturally be constructed in holography in a “measurement-based” setting. In this setting, we derive a flow equation for the operator reconstruction of a fixed code subspace operator using modular theory which can, in principle, be integrated to relate the reconstructed operators all along the flow. We observe a striking resemblance between our formulas for relational bulk reconstruction and the infinite-time limit of Connes cocycle flow, and take some steps towards making this connection more rigorous. We also provide alternative derivations of our reconstruction formulas in terms of a canonical reconstruction map we call the modular reflection operator.

Physics↗