Search NASA⌕ Search

SEARCH · Search NASA

Results for “EIGENFUNCTION”

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 19 records

On Eigenfunctions of the Fourier Transform

A nontrivial example of an eigenfunction in the sense of the theory of distributions for the planar Fourier transform was described by the authors in their previous work. In this paper, a method for obtaining other eigenfunctions is proposed. Positive homogeneous distributions in ℝ{sup n} of order −n/2 are considered, and it is shown that F(ω)|x|{sup −n/2}, |ω| = 1, is an eigenfunction in the sense of the theory of distributions of the Fourier transform if and only if F(ω) is an eigenfunction of a certain singular integral operator on the unit sphere of ℝ{sup n}. Since Y{sub m,n}{sup (k)}(ω)|x|{sup −n/2}, where Y{sub m,n}{sup (k)} denote the spherical functions of order m in ℝ{sup n}, are eigenfunctions of the Fourier transform, it follows that Y{sub m,n}{sup (k)} are eigenfunctions of the above-mentioned singular integral operator. In the planar case, all eigenfunctions of the Fourier transform of the form F(ω)|x|{sup −1} are described by means of the Fourier coefficients of F(ω).

97 MATHEMATICS AND COMPUTING↗

Quantum Circuits for the Preparation of Spin Eigenfunctions on Quantum Computers

The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wave functions that are relevant in the behavior of the system under study. Hamiltonian symmetries are important instruments used to classify relevant many-particle wave functions and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and approximate preparation of total spin eigenfunctions on quantum computers are presented. Two different strategies are discussed and compared: exact recursive construction of total spin eigenfunctions based on the addition theorem of angular momentum, and heuristic approximation of total spin eigenfunctions based on the variational optimization of a suitable cost function. The construction of these quantum circuits is illustrated in detail, and the preparation of total spin eigenfunctions is demonstrated on IBM quantum devices, focusing on three- and five-spin systems on graphs with triangle connectivity.

97 MATHEMATICS AND COMPUTING↗

Torus geometry eigenfunctions of an interacting multi-Landau-level Hamiltonian

A short-ranged, rotationally symmetric multi-Landau-level model Hamiltonian for strongly interacting electrons in a magnetic field was proposed [A. Anand et al., Phys. Rev. Lett. 126, 136601 (2021)] with the key feature that it allows exact many-body eigenfunctions on the disk not just for quasiholes but for all charged and neutral excitations of the entire Jain sequence filling fractions. We extend this to geometries without full rotational symmetry, namely, the torus and cylinder geometries, and present their spectra. Exact diagonalization of the interaction on the torus produces the low-energy spectra at filling fraction v = n/(2⁢pn + 1) that is identical, up to a topological (2⁢pn + 1)-fold multiplicity, to that of the integer quantum Hall spectra at v = n, for the incompressible state as well as all excitations. While the ansatz eigenfunctions in the disk geometry cannot be generalized to closed geometries such as torus or sphere, we show how to extend them to cylinder geometry. Meanwhile, we show eigenfunctions for charged excitations at filling fractions between 1/3 and 2/5 can be written on the torus and the spherical geometries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On Justification of the Asymptotics of Eigenfunctions of the Absolutely Continuous Spectrum in the Problem of Three One-Dimensional Short-Range Quantum Particles with Repulsion

The present paper offers a new approach to the construction of the coordinate asymptotics of the kernel of the resolvent of the Schrödinger operator in the scattering problem of three onedimensional quantum particles with short-range pair potentials. Within the framework of this approach, the asymptotics of eigenfunctions of the absolutely continuous spectrum of the Schrödinger operator can be constructed. In the paper, the possibility of a generalization of the suggested approach to the case of the scattering problem of N particles with arbitrary masses is discussed.

97 MATHEMATICS AND COMPUTING↗

VLASOV EIGENFUNCTION ANALYSIS OF SPACE-CHARGE AND BEAM-BEAM EFFECTS

Space-charge and beam-beam interaction affect both incoherent and coherent motion of particles potentially leading to instabilities and deterioration of the beam parameters. An overview of these phenomena will be given with an emphasis on the observable spectral characteristics and the mitigation methods of their harmful effects.

43 PARTICLE ACCELERATORS↗

Multigroup Constant Calculation with Static α-Eigenvalue Monte Carlo for Time-Dependent Neutron Transport Simulations

Herein, multigroup constants for deterministic methods that preserve the time-dependent physics of the neutron transport equations are derived. Alternative multigroup constant weighting spectra are discussed: (1) the fundamental k-eigenfunction, (2) the fundamental α-eigenfunction, and (3) a composite of several α-modes. To generate the fundamental α-eigenfunction for calculating the multigroup constants, a static fundamental α-eigenvalue method is implemented into the open source Monte Carlo code OpenMC. Several static and kinetic problems are devised to verify the implementations and to investigate the relative performance of the alternative multigroup constant weighting spectra. Results emphasize that as a multigroup constant weighting spectrum, the fundamental α-eigenfunction offers physical characteristics that make it advantageous (in producing accurate solutions) over the typically used fundamental k-eigenfunction.

42 ENGINEERING↗

Higher-Order Bessel Equations Integrable in Elementary Functions

The eigenfunction problem for a scalar Euler operator leads to an ordinary differential equation, which is an analog of higher-order Bessel equations. Its solutions are expressed through elementary functions in the case where the corresponding Euler operator can be factorized in a certain appropriate way. We obtain a formula describing such solutions. We consider the problem on common eigenfunctions of two Euler operators and present commuting Euler operators of orders 4, 6, and 10 and a formula for their common eigenfunction and also commuting operators of orders 6 and 9.

97 MATHEMATICS AND COMPUTING↗

Non-Hermitian topological states in 2D line-graph lattices: evolving triple exceptional points on reciprocal line graphs

Non-Hermitian (NH) topological states, such as the doubly-degenerate nodes dubbed as exceptional points (EPs) in Bloch band structure of 2D lattices driven by gain and loss, have attracted much recent interest. We demonstrate theoretically that in the three-site edge-centered lattices, i.e. the so-called line-graph lattices, such as kagome lattice which is a line graph of hexagonal lattice, there exist three types of triply-degenerate EPs evolving intriguingly on another set of line graphs in the reciprocal space. A single TEP (STEP) with ±1/3 topological charge moves faithfully along the edges of reciprocal line graphs with varying gain and loss, while two STEPs merge distinctively into one unconventional orthogonal double TEP (DTEP) with ±2/3 charge at the vertices, which is characterized with two ordinary self-orthogonal eigenfunctions but one surprising 'orthogonal' eigenfunction. Differently, in a modified line-graph lattice with an off-edge-center site, the ordinary coalesced state of DTEPs emerges with three identical self-orthogonal eigenfunctions. Such NH states and their evolution can be generally realized in various artificial systems, such as photonic and sonic crystals, where light and sonic vortex beams with different fractional twisting can be found. Our findings shed new light on fundamental understanding of gapless topological states in NH systems in terms of creation and evolution of high-order EPs, and open up new research directions to further link line graph and flow network theory coupled with topological physics, especially under non-equilibrium gain/loss conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data-driven modeling and control of dynamical systems using Koopman and Perron-Frobenius operators

This dissertation studies the data-driven modeling and control problem of nonlinear systems by exploiting the linear operator theoretic framework involving Koopman and Perro-Frobenius operator. A systematic linear-operator based controller design procedure has been established, which can be used to solve a variety of nonlinear control problems, including feedback stabilization using control Lyapunov functions, optimal quadratic regulation using Koopman eigenfunctions and convex optimization formulation of optimal control problem using P-F and Koopman operator approximation. As the core of data-driven modeling, we first propose a new algorithm for the finite-dimensional approximation of the linear transfer Koopman and Perron-Frobenius operator from time-series data. We argue that the existing approach for the finite-dimensional approximation of these transfer operators such as Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD) do not capture two important properties of these operators, namely positivity and Markov property. The algorithm we propose preserves these two properties. We call the proposed algorithm as naturally structured DMD (NSDMD) since it retains the inherent properties of these operators. Naturally structured DMD algorithm leads to a better approximation of the steady-state dynamics of the system regarding computing Koopman and Perron- Frobenius operator eigenfunctions and eigenvalues. However, preserving positivity property is critical for capturing the real transient dynamics of the system. This positivity property of the transfer operators and it's finite-dimensional approximation play an important role for controller and estimator design of nonlinear systems. To solve the feedback stabilization problem for nonlinear control systems, we tried to take advantage of the Koopman operator framework. The Koopman operator approach provides a linear representation for a nonlinear dynamical system and a bilinear representation for a nonlinear control system. The problem of feedback stabilization of a nonlinear control system is then transformed to the stabilization of a bilinear control system. We propose a control Lyapunov function (CLF)-based approach for the design of stabilizing feedback controllers for the bilinear system. The search for finding a CLF for the bilinear control system is formulated as a convex optimization problem. This leads to a schematic procedure for designing CLF-based stabilizing feedback controllers for the bilinear system and hence the original nonlinear system. Another advantage of the proposed controller design approach outlined in this dissertation is that it does not require explicit knowledge of system dynamics. In particular, the bilinear representation of a nonlinear control system in the Koopman eigenfunction space can be obtained from time-series data. Next, we study the optimal quadratic regulation problem for nonlinear systems. The linear operator theoretic framework involving the Koopman operator is used to lift the dynamics of nonlinear control system to an infinite-dimensional bilinear system. The optimal quadratic regulation problem for nonlinear system is formulated in terms of the finite-dimensional approximation of the bilinear system. A convex optimization-based approach is proposed for solving the quadratic regulator problem for bilinear system. We applied a variety of examples and compared the simulation results between our framework and conventional LQR control using linearized model. For more general optimal control problems, we provide a density-function based convex formulation for the optimal control problem of the nonlinear system. The convex formulation relies on the duality result in the stability theory of a dynamical system involving density function and Perron-Frobenius operator. The optimal control problem is formulated as an infinite-dimensional convex optimization program. The finite-dimensional approximation of the optimization problem relies on the recent advances made in the data-driven computation of the Koopman operator, which is dual to the Perron-Frobenius operator. Simulation results are presented to demonstrate the application of the developed framework.

Huang, Bowen↗

Properties of a potential energy matrix in oscillator basis

Matrix elements of potential energy are examined in detail. We consider a model problem — a particle in a central potential. The most popular forms of central potential are taken up, namely, square-well potential, Gaussian, Yukawa and exponential potentials. We study eigenvalues and eigenfunctions of the potential energy matrix constructed with oscillator functions. It is demonstrated that eigenvalues coincide with the potential energy in coordinate space at some specific discrete points. We establish approximate values for these points. It is also shown that the eigenfunctions of the potential energy matrix are the expansion coefficients of the spherical Bessel functions in a harmonic oscillator basis. We also demonstrate a close relation between the separable approximation and L{sup 2} basis (J-matrix) method for the quantum theory of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unitarization from geometry

We study the perturbative unitarity of scattering amplitudes in general dimensional reductions of Yang–Mills theories and general relativity on closed internal manifolds. For the tree amplitudes of the dimensionally reduced theory to have the expected high-energy behavior of the higher-dimensional theory, the masses and cubic couplings of the Kaluza–Klein states must satisfy certain sum rules that ensure there are nontrivial cancellations between Feynman diagrams. These sum rules give constraints on the spectra and triple overlap integrals of eigenfunctions of Laplacian operators on the internal manifold and can be proven directly using Hodge and eigenfunction decompositions. One consequence of these constraints is that there is an upper bound on the ratio of consecutive eigenvalues of the scalar Laplacian on closed Ricci-flat manifolds with special holonomy. This gives a sharp bound on the allowed gaps between Kaluza–Klein excitations of the graviton that also applies to Calabi–Yau compactifications of string theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Weighted Composition Operators for Learning Nonlinear Dynamics

Operator theoretic methods in dynamical system have been dominated by the use of Koopman operators and their continuous time counterparts, such as Koopman Generators and Liouville Operators. The advantage gained from their use primarily stems from the ability to extract subspaces and eigenfunctions within a space of observables that are invariant with respect to the Koopman operator over that space. When this occurs, a dynamic mode decomposition of the systems state provides a linear model for the dynamical system. Not all Koopman operators have eigenfunctions that may be exploited in this manner. However, the framework can still be leveraged for approximations using other operators. In this setting, we present a different operator for the study of dynamical systems, the weighted composition operator. These operators are compact for a wide range of dynamics and spaces, and through their interactions with occupation kernels and vector valued kernels, they admit an estimation of the underlying dynamics. Here, this manuscript presents a new algorithm for the data driven study of dynamical systems from data, and also provides two numerical experiments where convergence is achieved as a proof of concept.

97 MATHEMATICS AND COMPUTING↗

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗