Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hamiltonian structure”

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
NASA NTRS2003

Microscopic Modeling of Intersubband Optical Processes in Type II Semiconductor Quantum Wells: Linear Absorption

Linear absorption spectra arising from intersubband transitions in semiconductor quantum well heterostructures are analyzed using quantum kinetic theory by treating correlations to the first order within Hartree-Fock approximation. The resulting intersubband semiconductor Bloch equations take into account extrinsic dephasing contributions, carrier-longitudinal optical phonon interaction and carrier-interface roughness interaction which is considered with Ando s theory. As input for resonance lineshape calculation, a spurious-states-free 8-band kp Hamiltonian is used, in conjunction with the envelop function approximation, to compute self-consistently the energy subband structure of electrons in type II InAs/AlSb single quantum well structures. We demonstrate the interplay of nonparabolicity and many-body effects in the mid-infrared frequency range for such heterostructures.

Li, Jian-Zhong↗
NASA NTRSDate not supplied

Simulating a Low-Weight Encoding of the Fermi-Hubbard Model Feedback and Support with QAOA

Recent work [1] has proposed a compact fermion-to-qubit mapping for the Fermi-Hubbard model (FHM), potentially lowering the overhead needed to simulate such a system on quantum hardware. In prior work [2], we numerically investigated how well a quantum annealer could prepare ground states of the FHM encoded using the resource-efficient mapping. In order to compare these results to gate-model systems, we turn to a QAOA-style approach and replace the annealing schedule with both Trotterizations thereof and also more general Hamiltonian-variational ansatze based on the terms of the mapped FH Hamiltonian. We investigate the achievable ground state fidelities of these gate sequences, consider various optimizations of the ansatze structure, and offer resource estimations of the optimized QAOA protocols for near-term hardware.

quantum computing↗
DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2020

First-principles coupled cluster theory of the electronic spectrum of transition metal dichalcogenides

The electronic properties of two-dimensional transition metal dichalcogenides (2D TMDs) have attracted much attention during the last decade. Here we show how a diagrammatic ab initio coupled cluster singles and doubles (CCSD) treatment paired with a careful thermodynamic limit extrapolation in two dimensions can be used to obtain converged band gaps for monolayer materials in the MoS 2 family. We find CCSD gaps to lie in the upper range of the spread of $\textit{GW}$ approximation based on density functional theory (DFT) simulations, and also find slightly higher effective hole masses compared to previous reports. We also investigate the ability of CCSD to describe trion states, finding a reasonable qualitative structure, but poor excitation energies due to the lack of screening of three-particle excitations in the effective Hamiltonian. Our study provides an independent high-level benchmark of the role of many-body effects in 2D TMDs and showcases the potential strengths and weaknesses of diagrammatic coupled cluster approaches for realistic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

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

Reducing Qubit Requirements while Maintaining Numerical Precision for the Variational Quantum Eigensolver: A Basis-Set-Free Approach

We present a basis-set-free approach to the variational quantum eigensolver using an adaptive representation of the spatial part of molecular wave functions.Our approach directly determines system-specific representations of qubit Hamiltonians while fully omitting globally defined basis sets. In this work, we use directly determined pair-natural orbitals on the level of second-order perturbation theory. This results in compact qubit Hamiltonians with high numerical accuracy. We demonstrate initial applications with compact Hamiltonians on up to 22 qubits where conventional representation would for the same systems require 40–100 or more qubits. Here, we further demonstrate reductions in the quantum circuits through the structure ofthe pair-natural orbitals.

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

Structure of rotational bands in 109 Rh

Rotational bands in 109 Rh are investigated in a simple model and in the interacting boson-fermion model. We have developed a solvable extended transitional Hamiltonian by adding a two-configuration mixing term. Results suggest that 109 Rh is a good candidate for triaxiality and shape coexistence. Mixing between 3/2 + states with K = 1/2 and 3/2 is found to be weak, as evidenced by the E2 strengths.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Crystallography of hyperbolic lattices

Hyperbolic lattices are a revolutionary platform for tabletop simulations of holography and quantum physics in curved space and facilitate efficient quantum error correcting codes. Their underlying geometry is non-Euclidean, and the absence of Bloch's theorem precludes the straightforward application of the often indispensable energy band theory to study model Hamiltonians on hyperbolic lattices. Motivated by recent insights into hyperbolic band theory, we initiate a crystallography of hyperbolic lattices. Here we show that many hyperbolic lattices feature a hidden crystal structure characterized by unit cells, hyperbolic Bravais lattices, and associated symmetry groups. Using the mathematical framework of higher-genus Riemann surfaces and Fuchsian groups, we derive a list of example hyperbolic {p,q} lattices and their hyperbolic Bravais lattices, including five infinite families and several graphs relevant for experiments in circuit quantum electrodynamics and topolectrical circuits. This dramatically simplifies the computation of energy spectra of tight-binding Hamiltonians on hyperbolic lattices, from exact diagonalization on the graph to solving a finite set of equations in terms of irreducible representations. The significance of this achievement needs to be compared to the all-important role played by conventional Euclidean crystallography in the study of solids. We exemplify the high potential of this approach by constructing and diagonalizing finite-dimensional Bloch wave Hamiltonians.

2-dimensional systems↗
NASA NTRS1989

Distributed structural control using multilayered piezoelectric actuators. I - Analysis

Timoshenko beam theory is applied to beams with multiple layers of piezoelectric material attached. The model is developed using a Hamiltonian approach, and includes a model of the external electrical circuits to which the piezoelectric layers are attached, as well as a complete set of boundary conditions. The resulting model is then formulated in state space, and compared to a state space model of the standard Timoshenko beam.

Cudney, H. H.↗
DOE OSTI2022

Variational quantum eigensolver with reduced circuit complexity

Abstract The variational quantum eigensolver (VQE) is one of the most promising algorithms to find eigenstates of a given Hamiltonian on noisy intermediate-scale quantum devices (NISQ). The practical realization is limited by the complexity of quantum circuits. Here we present an approach to reduce quantum circuit complexity in VQE for electronic structure calculations. Our ClusterVQE algorithm splits the initial qubit space into clusters which are further distributed on individual (shallower) quantum circuits. The clusters are obtained based on mutual information reflecting maximal entanglement between qubits, whereas inter-cluster correlation is taken into account via a new “dressed” Hamiltonian. ClusterVQE therefore allows exact simulation of the problem by using fewer qubits and shallower circuit depth at the cost of additional classical resources, making it a potential leader for quantum chemistry simulations on NISQ devices. Proof-of-principle demonstrations are presented for several molecular systems based on quantum simulators as well as IBM quantum devices.

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

Beam dynamics in independent phased cavities

Linear accelerators containing the sequence of independently phased cavities with constant geometrical velocity along each structure are widely used in practice. The chain of cavities with identical cell lengths is utilized within a certain beam velocity range, with subsequent transformation to the next chain with higher cavity velocity. Design and analysis of beam dynamics in this type of accelerator are usually performed using numerical simulations. A full theoretical description of particle acceleration in an array of independent phased cavities has not been developed. In the present paper, we provide an analytical treatment of beam dynamics in such linacs employing Hamiltonian formalism. We begin our analysis with an examination of beam dynamics in an equivalent traveling wave of a single cavity, propagating within accelerating section with constant phase velocity. We then consider beam dynamics in arrays of cavities, utilizing an effective traveling wave propagating along with the whole accelerator with the velocity of synchronous (reference) particle. The analysis concluded with the determination of the matched beam conditions. Finally, we present a beam dynamics study in 805 MHz Coupled Cavity Linac of the LANSCE accelerator facility.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Approximation of nearly-periodic symplectic maps via structure-preserving neural networks

A continuous-time dynamical system with parameter ε is nearly-periodic if all its trajectories are periodic with nowhere-vanishing angular frequency as ε approaches 0. Nearly-periodic maps are discrete-time analogues of nearly-periodic systems, defined as parameter-dependent diffeomorphisms that limit to rotations along a circle action, and they admit formal U(1) symmetries to all orders when the limiting rotation is non-resonant. For Hamiltonian nearly-periodic maps on exact presymplectic manifolds, the formal U(1) symmetry gives rise to a discrete-time adiabatic invariant. In this paper, we construct a novel structure-preserving neural network to approximate nearly-periodic symplectic maps. This neural network architecture, which we call symplectic gyroceptron, ensures that the resulting surrogate map is nearly-periodic and symplectic, and that it gives rise to a discrete-time adiabatic invariant and a long-time stability. This new structure-preserving neural network provides a promising architecture for surrogate modeling of non-dissipative dynamical systems that automatically steps over short timescales without introducing spurious instabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Six-dimensional matching of intense beam with linear accelerating structure

Beam matching is a common technique that is routinely employed in accelerator design with the aim of minimizing beam losses and preservation of beam brightness. Despite being widely used, a full theoretical understanding of beam matching in 6D remains elusive. In this work, we present an analytical treatment of 6D beam matching of a high-intensity beam onto an RF structure. We begin our analysis within the framework of a linear model, and apply the averaging method to a set of 3D beam envelope equations. Accordingly, we obtain a matched solution that is comprised of smoothed envelopes and periodic terms, describing envelope oscillations with the period of the focusing structure. We then consider the nonlinear regime, where the beam size is comparable with the separatrix size. Stating with a Hamiltonian analysis in 6D phase space, we attain a self-consistent beam profile and show that it is significantly different from the commonly used ellipsoidal shape. Subsequently, we analyze the special case of an equilibrium with equal space charge depression between all degrees of freedom. Comparison of beam dynamics for equipartitioned, equal space charge depression, and equal emittances beams is given. Finally, we present experimental results on beam matching in the LANSCE linac.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

Initial value and two point boundary value solutions to the Clohessy-Wiltshire equations

The nonhomogeneous Clohessy-Wiltshire (C-W) equations are formulated and solved as an initial value problem in the form structure of linear systems theory. The state transition matrix (STM) and its inverse are obtained explicitly in both Newtonian and Hamiltonian form. It is shown that the STM for the C-2 equations possesses a special property making its inverse easily obtainable. Since solutions to the C-W equations are needed in two-point boundary value form to construct a good mission design tool for orbit transfer, the Lambert problem is solved in the context of the C-W equations.

Mullins, Larry D.↗
NASA NTRS2004

Spin lifetime tuning in zincblende heterostructures and applications to spin devices

We present analytical expressions for the D'yakonov-Perel' spin relaxation rates under the combined action of bulk and structural inversion asymmetry for zincblende heterostructures when terms up to linear and third order in k are included in the Hamiltonian. We see for heterostructures that, under the right conditions, the lowest-order-in-k component of the spin relaxation tensor can be made to vanish for all spin components at the same time. We study how the inclusion of terms of higher order in k affects these results.

spin filters↗