Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigenvalue”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence

Abstract The onset of turbulent heat transport at a higher temperature gradient than the critical gradient of linear instability, known as the Dimits shift, is a recurring feature of nonlinear simulations for magnetically confined fusion plasmas. Resonance in the nonlinear coupling between the modes that dominate energy transfer can lead to suppression of turbulence and transport above the linear critical gradient. As an expression of this resonance, gyrokinetic simulations show a quasi-coherent interaction between streamers and sidebands coupled through the zonal flow within the Dimits regime. This mechanism is further confirmed by use of artificial complex frequencies which break the resonance. By incorporating corresponding saturation physics, the standard quasilinear model for rapid head flux prediction is improved, which can now predict reduced heat flux in the Dimits regime. In particular, the triplet correlation time, the lifetime of the nonlinear interaction, is shown to be well-approximated by combinations of linear eigenvalues, and yields good representations of the heat flux variation both in and above the Dimits regime. Thus, a reduced but predictive model for transport near the critical gradient of zonal-flow saturated turbulence now exists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Real time detection of multiple stable MHD eigenmode growth rates towards kink/tearing modes avoidance in DIII-D tokamak plasmas

Abstract Real time detection of time evolving growth rates of multiple stable magnetohydrodynamic (MHD) eigenmodes has been achieved in DIII-D tokamak experiments via multi-mode three-dimensional (3D) active MHD spectroscopy. The measured evolution of the multi-modes’ growth rates is in good accordance with the variation of the plasma β N . Using experimental equilibria, resistive MARS-F simulations found the two least stable modes to have comparable growth rates to those experimentally measured. Real time and offline calculations of the modes’ growth rates show comparable results and indicate that cleaner system input and output signals will improve the accuracy of the real time stability detection. Moreover, the shortest real time updating time window of multi-mode eigenvalues can be about 2 ms in DIII-D experiments. This real time monitoring of stable, macroscopic kink and tearing modes thus provides an effective tool for avoidance of the most common causes of tokamak disruption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Gaps labeling theorem for the bubble-diamond self-similar graphs

Abstract Motivated by the appearance of fractals in several areas of physics, especially in solid state physics and the physics of aperiodic order, and in other sciences, including the quantum information theory, we present a detailed spectral analysis for a new class of fractal-type diamond graphs, referred to as bubble-diamond graphs, and provide a gap-labeling theorem in the sense of Bellissard for the corresponding probabilistic graph Laplacians using the technique of spectral decimation. Labeling the gaps in the Cantor set by the normalized eigenvalue counting function, also known as the integrated density of states, we describe the gap labels as orbits of a second dynamical system that reflects the branching parameter of the bubble construction and the decimation structure. The spectrum of the natural Laplacian on limit graphs is shown generically to be pure point supported on a Cantor set, though one particular graph has a mixture of pure point and singularly continuous components.

Physics↗
DOE OSTI2024

Integrating subsystem embedding subalgebras and coupled cluster Green’s function: a theoretical foundation for quantum embedding in excitation manifold

Here, in this study, we introduce a novel approach to coupled-cluster Green's function (CCGF) embedding by seamlessly integrating conventional CCGF theory with the state-of-the-art sub-system embedding sub-algebras coupled cluster (SES-CC) formalism. This integration focuses primarily on delineating the characteristics of the sub-system and the corresponding segments of the Green's function, defined explicitly by active orbitals. Crucially, our work involves the adaptation of the SES-CC paradigm, addressing the left eigenvalue problem through a distinct form of Hamiltonian similarity transformation. This advancement not only facilitates a comprehensive representation of the interaction between the embedded sub-system and its surrounding environment but also paves the way for the quantum mechanical description of multiple embedded domains, particularly by employing the emergent quantum flow algorithms. Our theoretical underpinnings further set the stage for a generalization to multiple embedded sub-systems. This expansion holds significant promise for the exploration and application of non-equilibrium quantum systems, enhancing the understanding of system–environment interactions. In doing so, the research underscores the potential of SES-CC embedding within the realm of quantum computations and multi-scale simulations, promising a good balance between accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Exploiting a derivative discontinuity estimate for accurate G0W0 ionization potentials and electron affinities

Abstract The GW approximation has become an important tool for predicting charged excitations of isolated molecules and condensed systems. Its popularity can be attributed to many factors, including a favorable scaling and relatively good accuracy. In practical applications, the GW is often performed as a one-shot perturbation known as G 0 W 0 . Unfortunately, G 0 W 0 suffers from a strong starting point dependence and is often not as accurate as one would need. Self-consistent GW methodologies alleviate these problems but come with a marked increase in computational cost. In this manuscript, we propose the use of an estimate of the exchange-correlation derivative discontinuity to provide a remarkably good starting point for G 0 W 0 calculations, yielding ionization potentials and electron affinities with eigenvalue self-consistent GW quality at no additional cost. We assess the quality of the resulting methodology with the GW 100 benchmark set and compare its advantages over other similar methods.

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

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2022

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2022

Clustering with general photo- z uncertainties: application to Baryon Acoustic Oscillations

ABSTRACT Photometric data can be analysed using the 3D correlation function ξp to extract cosmological information via e.g. measurement of the Baryon Acoustic Oscillations (BAO). Previous studies modeled ξp assuming a Gaussian photo-z approximation. In this work we improve the modeling by incorporating realistic photo-z distribution. We show that the position of the BAO scale in ξp is determined by the photo-z distribution and the Jacobian of the transformation. The latter diverges at the transverse scale of the separation s⊥, and it explains why ξp traces the underlying correlation function at s⊥, rather than s, when the photo-z uncertainty σz/(1+ z) ≳ 0.02. We also obtain the Gaussian covariance for ξp. Due to photo-z mixing, the covariance of ξp shows strong off-diagonal elements. The high correlation of the data causes some issues to the data fitting. None the less, we find that either it can be solved by suppressing the largest eigenvalues of the covariance or it is not directly related to the BAO. We test our BAO fitting pipeline using a set of mock catalogs. The data set is dedicated for Dark Energy Survey Year 3 (DES Y3) BAO analyses and includes realistic photo-z distributions. The theory template is in good agreement with mock measurement. Based on the DES Y3 mocks, ξp statistic is forecast to constrain the BAO shift parameter α to be 1.001 ± 0.023, which is well consistent with the corresponding constraint derived from the angular correlation function measurements. Thus, ξp offers a competitive alternative for the photometric data analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

AlgoSCR: an algorithm for solar contamination removal from radio interferometric data

Hydrogen intensity mapping is a new field in astronomy that promises to make three-dimensional maps of the matter distribution of the Universe using the redshifted 21cm line of neutral hydrogen gas (HI). Several ongoing and upcoming radio interferometers, such as Tianlai, CHIME, HERA, HIRAX, etc., are using this technique. These instruments are designed to map large swaths of the sky by drift scanning over periods of many months. One of the challenges of the observations is that the daytime data are contaminated by strong radio signals from the Sun. In the case of Tianlai, this results in almost half of the measured data being unusable. We try to address this issue by developing an algorithm for solar contamination removal (AlgoSCR) from the radio data. The algorithm is based on an eigenvalue analysis of the visibility matrix and hence is applicable only to interferometers. We apply AlgoSCR to simulated visibilities, as well as real daytime data from the Tianlai dish array. The algorithm can reduce strong solar contamination by about 95 per cent without seriously affecting other weaker sky signals and thus makes the data usable for certain applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Polaritons and excitons: Hamiltonian design for enhanced coherence

The primary questions motivating this report are: Are there ways to increase coherence and delocalization of excitation among many molecules at moderate electronic coupling strength? Coherent delocalization of excitation in disordered molecular systems is studied using numerical calculations. The results are relevant to molecular excitons, polaritons, and make connections to classical phase oscillator synchronization. In particular, it is hypothesized that it is not only the magnitude of electronic coupling relative to the standard deviation of energetic disorder that decides the limits of coherence, but that the structure of the Hamiltonian—connections between sites (or molecules) made by electronic coupling—is a significant design parameter. Inspired by synchronization phenomena in analogous systems of phase oscillators, some properties of graphs that define the structure of different Hamiltonian matrices are explored. The report focuses on eigenvalues and ensemble density matrices of various structured, random matrices. Some reasons for the special delocalization properties and robustness of polaritons in the single-excitation subspace (the star graph) are discussed. The key result of this report is that, for some classes of Hamiltonian matrix structure, coherent delocalization is not easily defeated by energy disorder, even when the electronic coupling is small compared to disorder.

Science & Technology - Other Topics↗
DOE OSTI2023

Kinematics and dynamics of disclination lines in three-dimensional nematics

An exact kinematic law for the motion of disclination lines in nematic liquid crystals as a function of the tensor order parameter Q is derived. Unlike other order parameter fields that become singular at their respective defect cores, the tensor order parameter remains regular. Following earlier experimental and theoretical work, the disclination core is defined to be the line where the uniaxial and biaxial order parameters are equal, or equivalently, where the two largest eigenvalues of Q cross. This allows an exact expression relating the velocity of the line to spatial and temporal derivatives of Q on the line, to be specified by a dynamical model for the evolution of the nematic. By introducing a linear core approximation for Q, analytical results are given for several prototypical configurations, including line interactions and motion, loop annihilation, and the response to external fields and shear flows. Behaviour that follows from topological constraints or defect geometry is highlighted. Finally, the analytic results are shown to be in agreement with three-dimensional numerical calculations based on a singular Maier–Saupe free energy that allows for anisotropic elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Co-designing Spectral Transformation Oracles with Hybrid Oscillator-Qubit Quantum Processors: From Algorithms to Compilation

We co-design a family of quantum eigenvalue transformation oracles that can be efficiently implemented on hybrid discrete- or continuous-variable (qubit or qumode) hardware. To illustrate the oracle’s representation-theoretic power and near-term experimental accessibility, we encode a Gaussian imaginary time-evolution spectral filter. As a result, we define a continuous linear combination of unitaries block encoding. Due to the ancillary qumode’s infinite-dimensional nature, continuous-variable qumodes constitute a powerful compilation tool for encoding continuous spectral functions without discretization errors while minimizing resource requirements. We then focus on the ubiquitous task of preparing eigenstates in quantum spin models. For completeness, we provide an end-to-end compilation which expresses high-level oracles in terms of an experimentally realizable instruction set architecture in both 1D and 2D. Finally, we examine the leading-order effects of physical errors and highlight open research directions. Our algorithms scale linearly with the spatial extent of the target system and are applicable to both near-term and large-scale quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Double copy root of Hawking thermality

The Hawking radiation spectrum from a collapsing null shell can be derived via the double copy of a simpler gauge theory calculation. Analyzing the non-abelian Yang-Mills root of this process, we demonstrate that the radiation spectrum is thermal in the color charge eigenvalue 𝜆, not energy. Considering the 𝑆⁢𝑈⁡(𝑁 𝑐 ) gauge theory in the large 𝑁 𝑐 limit, we find the differential spectrum ⅆ⁢𝑁/ⅆ⁢𝜆 is a product of the gravitationally familiar Planck-like factor and the color phase space density, modeled here as the Wigner semicircle from random matrix theory. This reveals that apparent energy thermality in gravity is the direct dual of charge thermality in its underlying non-abelian gauge theory.

Carrasco, John Joseph M. [Northwestern University,↗
DOE OSTI2025

Universality of Rényi Entropy in Conformal Field Theory

We use the thermal effective theory to prove that, for the vacuum state in any conformal field theory in 𝑑 dimensions, the 𝑛th Rényi entropy 𝑆$^{(𝑛)}_{𝐴}$ behaves as 𝑆$^{(𝑛)}_{𝐴}$ = [𝑓⁡/(2⁢𝜋⁢𝑛) 𝑑−1 ]⁢[Area⁡(∂𝐴)/(𝑑−2)⁢𝜀 𝑑−2 ]⁢(1+𝑂⁡(𝑛)) in the 𝑛 → 0 limit when the boundary of the entanglement domain 𝐴 is spherical with the UV cutoff 𝜀. The theory dependence is encapsulated in the cosmological constant 𝑓 in the thermal effective action. Using this result, we estimate the density of states for large eigenvalues of the modular Hamiltonian for the domain 𝐴. In two dimensions, we can use the hot spot idea, which describes the effective action in the high-temperature limit when the temperature is position-dependent, to derive more powerful formulas valid for arbitrary positive 𝑛. We discuss the difference between two and higher dimensions and clarify the applicability of the hot spot idea. We also use the thermal effective theory to derive an analog of the Cardy formula for boundary operators in higher dimensions.

Conformal field theory↗
DOE OSTI2025

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗
DOE OSTI2025

Lanczos Algorithm, the Transfer Matrix, and the Signal-to-Noise Problem

This Letter introduces a method for determining the energy spectrum of lattice quantum chromodynamics by applying the Lanczos algorithm to the transfer matrix and using a bootstrap generalization of the Cullum-Willoughby method to filter out spurious eigenvalues. Proof-of-principle analyses of the simple harmonic oscillator and the lattice quantum chromodynamics proton mass demonstrate that this method provides faster ground-state convergence than the “effective mass,” which is related to the power-iteration algorithm. Lanczos provides more accurate energy estimates than multistate fits to correlation functions with small imaginary times while achieving comparable statistical precision. Two-sided error bounds are computed for Lanczos results and guarantee that excited-state effects cannot shift Lanczos results far outside their statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗