Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pairing Hamiltonian”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Anomalous electrodynamics and quantum geometry in the Dirac-Harper model for a graphene bilayer

In this work, graphene bilayers with layer antisymmetric strains are studied using the Dirac-Harper model for a pair of single layer Dirac Hamiltonians coupled by a one-dimensional moiré-periodic interlayer tunneling amplitude. This model hosts low energy, nearly dispersionless bands near charge neutrality that support anomalous polarizations of its charge multipole distributions. These are analyzed introducing a generalized Berry curvature that encodes the field-induced dynamics of multipole fields allowed in a chiral medium with time reversal symmetry. The formulation identifies a reciprocity relation between responses to layer-symmetric and layer-antisymmetric in-plane electric fields and reveals momentum-space quantum oscillations produced by a spatial pattern of band inversions on the moiré scale.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ideal weak topological insulator and protected helical saddle points

Here we report the design and observation of an ideal weak topological insulator in a quasi-one-dimensional bismuth halide, Bi 4 I 1.2 Br 2.8 . Via angle-resolved photoemission spectroscopy, we identify that Bi 4 I 1.2 Br 2.8 hosts topological surface states on the (100)′ side surface in the form of two anisotropic π-offset Dirac cones while topologically dark on the (001) top surface. The results fully determine a unique side-surface Hamiltonian and thereby identify two pairs of nondegenerate helical saddle points. The fact that both the surface Dirac and helical saddle points are in the global bulk band gap of 195 meV suggests the potential of this system to be a fertile ground for topological many-body physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fermionic mean-field theory as a tool for studying spin Hamiltonians

The Jordan–Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Furthermore, Jordan–Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J 1 –J 2 Heisenberg models, as well as to the pairing or reduced Bardeen–Cooper–Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

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

The 1 S 0 Pairing Gap in Neutron Matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2026

Component-Level Inverse Design of Transmon Qubits Using Neural Networks

Designing a superconducting qubit to realize specific Hamiltonian parameters typically requires iterating through a time and compute-intensive forward loop in which the designer chooses a layout geometry, simulates it, extracts circuit parameters such as capacitances, and refines the geometry. We study the inverse version of this task using a neural-network workflow that maps target Hamiltonian parameters directly to component-level layout parameters, which we subsequently demonstrate on a planar transmon layout. During training, we pair the inverse model with a frozen forward surrogate model and evaluate the loss in Hamiltonian space rather than in layout-parameter space. In validation against a conventional EM solver, 97% of generated designs produce usable geometries, and the inverse-plus-surrogate pipeline reaches mean percent errors of 0.73% for qubit frequency and 1.58% for anharmonicity, comparable to or below the fabrication and simulation-to-measurement uncertainty expected for academic-process transmon devices of this type. A single pipeline query takes ~60 ms on CPU, versus ~2 min for a conventional EM capacitance extraction on the same hardware, a speedup of approximately 2,000x. Batching minimizes the AI model inference overhead, reducing the runtime to 3.1 microseconds per sample on CPU and 2.6 microseconds per sample on GPU at a batch size of 2048, resulting in speedups of 3.9 x 10^7 and 4.6 x 10^7, respectively, relative to a single conventional CPU EM extraction. Our results indicate that component-level inverse design usefully extends and complements conventional EM simulation, including for small datasets on the order of 1,000 samples.

Seidel, Olivia [Fermilab; Texas U., Arlington]↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Manipulating Growth and Propagation of Correlations in Dipolar Multilayers: From Pair Production to Bosonic Kitaev Models

Here, we study the nonequilibrium dynamics of dipoles confined in multiple stacked two-dimensional layers realizing a long-range interacting quantum spin 1/2 XXX model. We demonstrate that strong in-plane interactions can protect a manifold of collective layer dynamics. This then allows us to map the many-body spin dynamics to bosonic models. In a bilayer configuration we show how to engineer the paradigmatic two-mode squeezing Hamiltonian known from quantum optics, resulting in exponential production of entangled pairs and generation of metrologically useful entanglement from initially prepared product states. In multilayer configurations we engineer a bosonic variant of the Kitaev model displaying chiral propagation along the layer direction. Our study illustrates how the control over interactions, lattice geometry, and state preparation in interacting dipolar systems uniquely afforded by AMO platforms such as Rydberg and magnetic atoms, polar molecules, or trapped ions allows for the control over the temporal and spatial propagation of correlations for applications in quantum sensing and quantum simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stochastic exciton-scattering theory of optical line shapes: Renormalized many-body contributions

Spectral line shapes provide a window into the local environment coupled to a quantum transition in the condensed phase. In this paper, we build upon a stochastic model to account for non-stationary background processes produced by broad-band pulsed laser stimulation, as distinguished from those for stationary phonon bath. In particular, we consider the contribution of pair-fluctuations arising from the full bosonic many-body Hamiltonian within a mean-field approximation, treating the coupling to the system as a stochastic noise term. Herein, using the Itô transformation, we consider two limiting cases for our model, which lead to a connection between the observed spectral fluctuations and the spectral density of the environment. In the first case, we consider a Brownian environment and show that this produces spectral dynamics that relax to form dressed excitonic states and recover an Anderson–Kubo-like form for the spectral correlations. In the second case, we assume that the spectrum is Anderson–Kubo like and invert to determine the corresponding background. Using the Jensen inequality, we obtain an upper limit for the spectral density for the background. The results presented here provide the technical tools for applying the stochastic model to a broad range of problems.

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

Quantum simulations for strong-field QED

Quantum field theory in the presence of strong background fields contains interesting problems where quantum computers may someday provide a valuable computational resource. In the noisy intermediate-scale quantum era it is useful to consider simpler benchmark problems in order to develop feasible approaches, identify critical limitations of current hardware, and build new simulation tools. Here we perform quantum simulations of strong-field QED (SFQED) in 3 + 1 dimensions, using real-time nonlinear Breit-Wheeler pair production as a prototypical process. The strong-field QED Hamiltonian is derived and truncated in the Furry-Volkov mode expansion, and the interactions relevant for Breit-Wheeler are transformed into a quantum circuit. Quantum simulations of a “null double slit” experiment are found to agree well with classical simulations following the application of various error mitigation strategies, including an asymmetric depolarization algorithm which we develop and adapt to the case of Trotterization with a time-dependent Hamiltonian. We also discuss longer-term goals for the quantum simulation of SFQED. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Coulomb confinement in the Hamiltonian limit

The Gribov-Zwanziger scenario attributes the phenomenon of confinement to the instantaneous interaction term in the QCD Hamiltonian in the Coulomb gauge. For a static quark-antiquark pair, it leads to a potential energy that increases linearly with the distance between them. Lattice studies of the SU(2) Yang-Mills theory determined the corresponding (Coulomb) string tension for sources in the fundamental representation, 𝜎 𝐶 , to be about 3 times larger than the Wilson loop string tension, 𝜎 𝐹 . It is far above the Zwanziger variational bound, 𝜎 𝐶 ≥ 𝜎 𝐹 . We argue that the value often reported in the literature is artificially inflated. We examine the lattice definition of the instantaneous potential, find the source of the string tension’s enhancement, and perform its improved determination in SU(2) lattice gauge theory. We report our conservative estimate for the value of the Coulomb string tension as 𝜎 𝐶 /𝜎 𝐹 = 2.0 ± 0.4 and discuss its phenomenological implications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Coupling of Nondegenerate Topological Modes in Nitrogen Core-Doped Graphene Nanoribbons

Nitrogen core-doping of graphene nanoribbons (GNRs) allows trigonal planar carbon atoms along the backbone of GNRs to be substituted by higher-valency nitrogen atoms. The excess valence electrons are injected into the π-orbital system of the GNR, thereby changing not only its electronic occupation but also its topological properties. We have observed this topological change by synthesizing dilute nitrogen core-doped armchair GNRs with a width of five atoms (N 2 -5-AGNRs). The incorporation of pairs of trigonal planar nitrogen atoms results in the emergence of topological boundary states at the interface between doped and undoped segments of the GNR. These topological boundary states are offset in energy by approximately ΔE = 300 meV relative to the topological end states at the termini of finite 5-AGNRs. Scanning tunneling microscopy (STM) and spectroscopy (STS) reveal that for finite GNRs the two types of topological states can interact through a linear combination of orbitals, resulting in a pair of asymmetric hybridized states. This behavior is captured by an effective Hamiltonian of nondegenerate diatomic molecules, where the analogous interatomic hybridization interaction strength is tuned by the distance between GNR topological modes.

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

Schematic model for induced fission in a configuration-interaction approach

We model fission at barrier-top energies in a simplified model space that permits comparison of different components of the residual nucleon-nucleon interaction. The model space is built on particle-hole excitations of reference configurations. These are Slater determinants of uniformly spaced orbitals characterized only by their quantum numbers and orbital energies. The residual interaction in the Hamiltonian includes the diabatic interaction connecting similar orbitals at different deformations, the pairing interaction between like nucleons, and a schematic off-diagonal neutron-proton interaction. We find that the fission reaction probability is sensitive to the off-diagonal neutron-proton interaction much more than to the pairing and the diabatic interactions. In particular, the transmission coefficients become insensitive to the strength of the pairing interaction when the neutron-proton interaction is large. Here, we also find that the branching ratio is insensitive to the final-state scission dynamics, as is assumed in the well-known Bohr-Wheeler theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Trapped-ion quantum simulation of collective neutrino oscillations

It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of N neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. In conclusion, we finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evidence for valence-bond pairing in a one-dimensional two-orbital system

Valence bond (VB) states as the formation mechanism of Cooper pairs, eventually leading to high-temperature superconductivity, remain a controversial topic. Although various VB-like states find variational relevance in the description of specific spin models and quantum spin liquids, in the realm of many-body fermionic Hamiltonians, the evidence for such states as ground-state wave functions remains elusive, challenging the valence-bond pairing mechanism. Here, we present evidence of a VB ground state with pairing tendencies, particularly at finite doping. Here, we achieved this for the generic two-orbital Hubbard model in low dimension, where the VB states can be associated with the presence of the topological order manifested by edge states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Doping topological Dirac semimetal with magnetic impurities: Electronic structure of Mn-doped Cd 3 ⁢As 2

The prospect of transforming a Dirac topological semimetal (TSM) into a Weyl TSM phase, following doping by magnetic impurities, is central to TSM applications. The magnetic field from polarized 𝑑 levels of magnetic impurities produces a field with a sharp local structure. To what extent magnetic impurities act in the same manner as an applied field and what are the effects of such a field on the electronic structure of a Dirac TSM are the subject of this paper. We present electronic structure calculations of bulk Cd 3 ⁢As 2 with substitutional doping of Mn impurities in the dilute alloy range. Quasiparticle 𝐺⁡𝑊 (QS⁢𝐺⁡𝑊) ab initio electronic structure calculations are used in conjunction with 𝑘 · 𝑝 model Hamiltonian calculations. As expected, we observe the splitting of the Dirac points into pairs of Weyl points following the doping with Mn. We also show that the electronic structure of Mn-doped Cd 3⁢ As 2 can be emulated by the electronic structure of pristine Cd 3⁢ As 2 with an appropriate external magnetic field. Some properties of the conductivity of bulk Cd 3 ⁢As 2 for different magnetic field orientations are also investigated. Our results inform future opportunities for unique device functionality based on band structure tuning not found in conventional magnetic Weyl TSM.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Long-range three-dimensional magnetic structures of the spin S =1 hexamer cluster fedotovite-like A 2 Cu 3 O ( SO 4 ) 3 ( A 2 = K 2 , NaK, Na 2 ): A neutron diffraction study

The crystal and magnetic structures of the spin S = 1 hexamer cluster fedotovite-like A 2 Cu 3 O(SO 4 ) 3 (A 2 = K 2 , NaK, Na 2 ) were studied by neutron powder diffraction at temperatures 1.6–290 K. The crystal structures in all compounds are well refined in the monoclinic space group C2/c. The basic magnetic units of the compounds are copper hexamers, which are coupled by weak superexchange interactions giving rise to three-dimensional long-range magnetic order below 3.0 < T N < 4.7 K. We have found that for A 2 = K 2 and NaK the propagation vector of the magnetic structure is k = [0,0,0], and the coupling of the Cu hexamers is ferromagnetic (FM) along the ab diagonal and antiferromagnetic along the bc diagonal. In contrast, for A = Na the propagation vector is k = [0,1,0], and the Cu hexamers are coupled antiferromagnetically (AFM) along the ab diagonal. Further, the hexamers are formed by three Cu pairs arranged along the b axis. The calculated spin expectation values $\langle s \rangle$ for the simplest symmetric spin Hamiltonian (obtained from inelastic neutron spectroscopy) of the isolated hexamers in the mean field amounted to $\langle s \rangle$ = 3/8 for the side spins Cu1 and Cu2 and $\langle s \rangle$ = 1/4 for Cu3 in the middle. The Cu spins are FM coupled in pairs and AFM between neighboring pairs. The experimental magnetic moments of the Cu 2+ ions turn out to be not completely collinear due to spin frustrations within the weak interhexamer interactions. The sizes of magnetic moments of Cu in the hexamers determined from the diffraction data are in fair agreement with the calculated values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗