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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantifying the phase diagram and Hamiltonian of S = 1/2 kagome antiferromagnets: bridging theory and experiment

Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D, which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D. Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.

Jiang, Shengtao [SLAC National Accelerator Laborat↗
DOE OSTI2023

Mott insulators with boundary zeros

The topological classification of electronic band structures is based on symmetry properties of Bloch eigenstates of single-particle Hamiltonians. In parallel, topological field theory has opened the doors to the formulation and characterization of non-trivial phases of matter driven by strong electron-electron interaction. Even though important examples of topological Mott insulators have been constructed, the relevance of the underlying non-interacting band topology to the physics of the Mott phase has remained unexplored. Here, we show that the momentum structure of the Green’s function zeros defining the “Luttinger surface" provides a topological characterization of the Mott phase related, in the simplest description, to the one of the single-particle electronic dispersion. Considerations on the zeros lead to the prediction of new phenomena: a topological Mott insulator with an inverted gap for the bulk zeros must possess gapless zeros at the boundary, which behave as a form of “topological antimatter” annihilating conventional edge states. Placing band and Mott topological insulators in contact produces distinctive observable signatures at the interface, revealing the otherwise spectroscopically elusive Green’s function zeros.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Applications of reduced-basis methods to the nuclear single-particle spectrum

Reduced-basis methods provide a powerful framework for building efficient and accurate emulators. Although widely applied in many fields to simplify complex models, reduced-basis methods have only been recently introduced into nuclear physics. In this Letter we build an emulator to study the single-particle structure of atomic nuclei. By scaling a suitable mean-field Hamiltonian, a “universal” reduced basis is constructed capable of accurately and efficiently reproduce the entire single-particle spectrum of a variety of nuclei. Indeed, the reduced-basis model reproduces both ground- and excited-state energies as well as the associated wave functions with remarkable accuracy. Here our results bode well for more demanding applications that use Bayesian optimization to calibrate nuclear energy density functionals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Theoretical antiferromagnetism of ordered face-centered cubic Cr-Ni alloys

Contrary to prior calculations, the Ni-rich ordered structures of the Cr-Ni alloy system are found to be antiferromagnetic under semilocal density-functional theory. The optimization of local magnetic moments significantly increases the driving force for the formation of CrNi 2 , the only experimentally observed intermetallic phase. This structure's ab initio magnetism appears well described by a Heisenberg Hamiltonian with longitudinal spin fluctuations; itinerant Cr moments are induced only by the strength of exchange interactions. The role of magnetism at temperature is less clear and several scenarios are considered based on a review of experimental literature, specifically a failure of the theory, the existence of an overlooked magnetic phase transition, and the coupling of antiferromagnetism to chemical ordering. In conclusion, implications for related commercial and high-entropy alloys are discussed for each case.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multimodal Approach Reveals the Symmetry-Breaking Pathway to the Broken Helix in EuIn 2 ⁢As 2

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to EuIn 2 ⁢As 2 , a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of EuIn 2 ⁢As 2 . We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in EuIn 2 ⁢As 2 , but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in EuIn 2 ⁢As 2 . Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Light-induced electron dynamics of molecules in cavities: comparison of model Hamiltonians

The rapidly expanding field of polaritonic chemistry requires accurate theoretical simulations to understand new phenomena at the atomic scale. Computing the optoelectronic properties of molecules using established electronic structure methods is a careful balance of accuracy and computational expense, and expanding these methods to quantum electrodynamics to describe coupled cavity-molecule systems is an active topic of development. Key to these methods are the Hamiltonian operators representing the photon cavity modes. The recently introduced quantum electrodynamics time-dependent configuration interaction (QED-TDCI) method allows for the combination of electron dynamics simulations with quantum electrodynamics, enabling the simulation of time-dependent optoelectronic properties of cavity-molecule systems. Using this method, a comparison of two many-state QED Hamiltonians—the Pauli-Fierz and quantum Rabi model Hamiltonians—is presented, with a particular focus on time-dependent properties in applied electric fields.

Peyton, Benjamin G. (ORCID:0000000163190751)↗
DOE OSTI2022

Trotter Errors from Dynamical Structural Instabilities of Floquet Maps in Quantum Simulation

We study the behavior of errors in the quantum simulation of spin systems with long-range multibody interactions resulting from the Trotter-Suzuki decomposition of the time-evolution operator. We identify a regime where the Floquet operator underlying the Trotter decomposition undergoes sharp changes even for small variations in the simulation step size. This results in a time evolution operator that is very different from the dynamics generated by the targeted Hamiltonian, which leads to a proliferation of errors in the quantum simulation. These regions of sharp change in the Floquet operator, referred to as structural instability regions, appear typically at intermediate Trotter step sizes and in the weakly interacting regime, and are thus complementary to recently revealed quantum chaotic regimes of the Trotterized evolution [L. M. Sieberer et al. npj Quantum Inf. 5, 78 (2019); M. Heyl, P. Hauke, and P. Zoller, Sci. Adv. 5, eaau8342 (2019)]. We characterize these structural instability regimes in p-spin models, transverse-field Ising models with all-to-all p-body interactions, and analytically predict their occurrence based on unitary perturbation theory. We further show that the effective Hamiltonian associated with the Trotter decomposition of the unitary time-evolution operator, when the Trotter step size is chosen to be in the structural instability region, is very different from the target Hamiltonian, which explains the large errors that can occur in the simulation in the regions of instability. These results have implications for the reliability of near-term gate-based quantum simulators, and reveal an important interplay between errors and the physical properties of the system being simulated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Massless and massive photons within basis light-front quantization

We compute the structure function, transverse momentum dependent parton distributions (TMDs), and generalized parton distributions (GPDs) for the photon from the light-front quantum electrodynamics (QED) Hamiltonian, determined for its constituent bare photon and electron-positron Fock components within the basis light-front quantization framework. While the structure function and GPDs align well with expectations, some oscillatory behavior in the TMDs emerges as a finite-basis artifact. We also compute the GPD and the corresponding structure function for the longitudinal polarization mode of the massive photon. We further study the photon mass dependence of the calculated observables. After performing nonperturbative renormalization, we obtain a good quality description of the perturbative QED properties of the photon, as well as good agreement with the experimental data for the massless photon’s structure function.

Ab initio calculations↗
DOE OSTI2019

Efficient Multiscale Optoelectronic Prediction for Conjugated Polymers

Conjugated polymers represent a high-potential material class due to the tunability of their optoelectronic properties via straightforward processing protocols. However, to in silico tailor these optoelectronic properties, one must compute the conformationally dependent electronic structure over mesoscopic length and time scales. This task represents a challenging multiscale computational problem in which quantum chemistry and atomistic or coarse-grained molecular dynamics must be integrated. Recently, we introduced a methodology, called electronic coarse graining (ECG), that utilizes artificial neural networks to compute ab initio quality electronic structures using only the reduced degrees of freedom of coarse-grained models. Here, we adapt ECG to variable-molecular weight conjugated oligomers by casting the problem in the framework of sequence detection. A machine-learning (ML) architecture employing one-dimensional convolutional neural networks and long short-term memory networks is utilized to compute ground-state orbital energies, charge density distributions, and optical spectra solely from the coarse-grained model's configurational degrees of freedom. Robust molecular weight transferability for ECG is established via a A-ML approach that leverages model electronic Hamiltonians for ground and excited-state property determination. The accuracy and transferability of the ECG methodology opens the door for scalable optoelectronic property prediction in conjugated polymers directly from coarse-grained degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydrodynamics of quantum spin liquids

Quantum spin liquids are topological states of matter that arise in frustrated quantum magnets at low temperatures. At low energies, such states exhibit emergent gauge fields and fractionalized quasiparticles and can also possess enhanced global symmetries compared to their parent microscopic Hamiltonians. Here we study the consequences of this emergent gauge and symmetry structure for the hydrodynamics of quantum spin liquids. Specifically, we analyze two cases, the U(1) spin liquid with a Fermi surface and the SU(4)-symmetric “algebraic” spin liquid. We show that the emergent degrees of freedom in the spin liquid phase lead to a variety of additional hydrodynamic modes compared to the high-temperature paramagnetic phase. We identify a hydrodynamic regime for the internal U(1) gauge field common to both states, characterized by slow diffusion of the internal transverse photon.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Static and dynamical properties of the spin-$\frac{5}{2}$ nearly ideal triangular lattice antiferromagnet $\mathrm{Ba_3}$ $\mathrm{MnSb_2}$ $\mathrm{O_9}$

Here, we study the ground state and spin excitations in Ba 3 MnSb 2 O 9 , an easy-plane S = 5/2 triangular lattice antiferromagnet. By combining single-crystal neutron scattering, electric spin resonance (ESR), and spin wave calculations, we determine the frustrated quasi-two-dimensional spin Hamiltonian parameters describing the material. While the material has a slight monoclinic structural distortion, which could allow for isosceles-triangular exchanges and biaxial anisotropy by symmetry, we observe no deviation from the behavior expected for spin waves in the in-plane $120$° state. Even the easy-plane anisotropy is so small that it can only be detected by ESR in our study. In conjunction with the quasi-two-dimensionality, our study establishes that Ba 3 MnSb 2 O 9 is a nearly ideal triangular lattice antiferromagnet with the quasiclassical spin S = 5/2, which suggests that it has the potential for an experimental study of Z- or Z 2 -vortex excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Precision Mass Measurement of the Proton Dripline Halo Candidate 22 Al

Here, we report the first mass measurement of the proton-halo candidate 22 Al performed with the low energy beam ion trap facility’s 9.4 T Penning trap mass spectrometer at facility for rare isotope beams. This measurement completes the mass information for the lightest remaining proton-dripline nucleus achievable with Penning traps. 22 Al has been the subject of recent interest regarding a possible halo structure from the observation of an exceptionally large isospin asymmetry [J. Lee et al ., Large isospin asymmetry in Si22/O22 Mirror Gamow-Teller transitions reveals the halo structure of 22 Al , Phys. Rev. Lett. 125, 192503 (2020).]. The measured mass excess value of ME=18 092.5⁢(3) keV, corresponding to an exceptionally small proton separation energy of 𝑆𝑝=100.4⁢(8) keV, is compatible with the suggested halo structure. Our result agrees well with predictions from 𝑠⁢𝑑-shell USD Hamiltonians. While USD Hamiltonians predict deformation in the 22 Al ground state with minimal 1⁢𝑠 1/2 occupation in the proton shell, a particle-plus-rotor model in the continuum suggests that a proton halo could form at large quadrupole deformation. These results emphasize the need for a charge radius measurement to conclusively determine the halo nature.

binding energy↗
DOE OSTI2025

Quantum Mechanics of a Spherically Symmetric Causal Diamond in Minkowski Spacetime

We construct the phase space of a spherically symmetric causal diamond in ( d + 2 )-dimensional Minkowski spacetime. Utilizing the covariant phase space formalism, we identify the relevant degrees of freedom that localize to the d -dimensional bifurcate horizon and, upon canonical quantization, determine their commutators. On this phase space, we find two Iyer-Wald charges. The first of these charges, proportional to the area of the causal diamond, is responsible for shifting the null time along the horizon and has been well documented in the literature. The second charge is much less understood, being integrable for d ≥ 2 only if we allow for field-dependent diffeomorphisms and is responsible for changing the size of the causal diamond. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Dynamics of almost strong edge modes in spin chains away from integrability

Results are presented for the dynamics of an almost strong edge mode which is the quasistable Majorana edge mode occurring in nonintegrable spin chains. The dynamics of the edge mode is studied using exact diagonalization, and compared with time evolution with respect to an effective semi-infinite model in Krylov space obtained from the recursion method. The effective Krylov Hamiltonian is found to resemble a spatially inhomogeneous Su-Schrieffer-Heeger model where the hopping amplitude increases linearly with distance into the bulk, typical of thermalizing systems, but also has a staggered or dimerized structure superimposed on it. Here, the nonperturbatively long lifetime of the edge mode is shown to be due to this staggered structure which diminishes the effectiveness of the linearly growing hopping amplitude. On taking the continuum limit of the Krylov Hamiltonian, the edge mode is found to be equivalent to the quasistable mode of a Dirac Hamiltonian on a half line, with a mass which is nonzero over a finite distance, before terminating into a gapless metallic bulk. The analytic estimates are found to be in good agreement with the numerically obtained lifetimes of the edge mode.

1-dimensional systems↗