Search NASA⌕ Search

SEARCH · Search NASA

Results for “Exact Diagonalization”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Hidden spin current in doped Mott antiferromagnets

We investigate the nature of doped Mott insulators using exact diagonalization and density matrix renormalization group methods. Persistent spin currents are revealed in the ground state, which are concomitant with a nonzero total momentum or angular momentum associated with the doped hole. The latter determines a nontrivial ground state degeneracy. By further making superpositions of the degenerate ground states with zero or unidirectional spin currents, we show that different patterns of spatial charge and spin modulations will emerge. Here, such anomaly persists for the odd numbers of holes, but the spin current, ground state degeneracy, and charge/spin modulations completely disappear for even numbers of holes, with the two-hole ground state exhibiting a d-wave symmetry. An understanding of the spin current due to a many-body Berry-like phase and its influence on the momentum distribution of the doped holes will be discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2020

Quench dynamics of two-component dipolar fermions subject to a quasiperiodic potential

Motivated by recent experiments in fermionic polar gases, we study the nonequilibrium dynamics of two-component dipolar fermions subject to a quasiperiodic potential. We investigate the localization of charge and spin degrees of freedom time evolving with a long-range spin-SU(2) symmetric fermionic Hamiltonian by calculating experimentally accessible dynamical observables. To study the nonequilibrium dynamics, we start the time evolution with two initial states at half-filling: (i) product state with doublons |↑↓0↑↓0↑↓0↑↓0↑↓$\rangle$ and (ii) product state with singlons |↑↓↑↓↑↓↑↓↑↓$\rangle$. We carried out the real-time evolution via the fermionic Hamiltonian using exact diagonalization (ED) and the time-dependent variational principle (TDVP) for finite matrix product states (MPS), within experimentally relevant timescales. For the product state with doublons, we observe a delocalized to localized phase transition varying disorder strengths by monitoring the decay of charge imbalance with time. For the long-range interacting Hamiltonian of our focus, and in the presence of strong enough disorder, starting the time evolution with singlons we find a strong reduction in the spin delocalization, contrary to results of previous studies using the disordered short-range (on-site) Hubbard model with SU(2) symmetry. In addition to the quench dynamics, we also demonstrate the localization of charge and spin in the full energy spectrum of the long-range spin-SU(2) symmetric Hamiltonian by monitoring spin and charge autocorrelation functions. Finally, our predictions for localization of both charge and spin should be observable in ultracold experiments with fermionic dipolar atoms subject to a quasiperiodic potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Valley Stoner instability of the composite Fermi sea

We study two-component electrons in the lowest Landau level at total filling factor ν T = 1/2 with anisotropic mass tensors and principal axes rotated by π/2 as realized in aluminum arsenide (AlAs) quantum wells. Combining exact diagonalization and the density matrix renormalization group we demonstrate that the system undergoes a quantum phase transition from a gapless state in which both flavors are equally populated to another gapless state in which all the electrons spontaneously polarize into a single flavor beyond a critical mass anisotropy of m x /m y ~7. We propose that this phase transition is a form of itinerant Stoner transition between a two-component and a single-component composite Fermi sea states and describe a set of trial wave functions which successfully capture the quantum numbers and shell filling effects in finite size systems as well as providing a physical picture for the energetics of these states. Our estimates indicate that the composite Fermi sea of AlAs is the analog of an itinerant Stoner magnet with a finite spontaneous valley polarization. Here, we pinpoint experimental evidence indicating the presence of Stoner magnetism in the Jain states surrounding ν = 1/2.

2-dimensional systems↗
DOE OSTI2021

Z 2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z 2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first order, in contrast to the conventional Z 2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic X X interaction while keeping the local gauge symmetry intact. We find that the Z 2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.

2-dimensional systems↗
DOE OSTI2022

Trions in twisted bilayer graphene

The strong-coupling phase diagram of magic-angle twisted bilayer graphene (TBG) predicts a series of exact one-particle charge ±1 gapped excitations on top of the integer-filled ferromagnetic ground states. Finite-size exact diagonalization studies showed that these are the lowest charge ±1 excitations in the system (for 10 nm screening length), with the exception of charge +1 at filling -1 in the chiral limit. In the current paper we show that this “trion bound state,” a 3-particle, charge 1 excitation of the insulating ferromagnetic ground state of the projected Hamiltonian of TBG, is the lowest charge +1 overall excitation at ν =-1, and also for some large (≈20 nm) screening lengths at ν =-2 in the chiral limit and with very small binding energy. At other fillings, we show that trion bound states do exist, but only for momentum ranges that do not cover the entire moiré Brillouin zone. The trion bound states (at different momenta) exist for finite parameter range w 0 /w 1 but they all disappear in the continuum far below the realistic values of w 0 /w 1 = 0.8. Thus we find the conditions for the existence of the trion bound state, a good variational wave function for it, and investigate its behavior for different screening lengths, at all integer fillings, on both the electron and hole sides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detection of chiral spin fluctuations driven by frustration in Mott insulators

Topologically ordered states, such as chiral spin liquids, have been proposed as candidates that host fraction alized excitations. However, detecting chiral character or proximity to these nontrivial states remains a challenge. Resonant Raman scattering can be a powerful tool for detecting chiral fluctuations, as the 𝐴 2⁢𝑔 channel probes excitations with broken time-reversal symmetry and local chiral order. Here, we use exact diagonalization to characterize the resonant 𝐴 2⁢𝑔 channel, alongside two-magnon scattering in 𝐵 1⁢𝑔 and 𝐸 𝑔 channels, for the Hubbard model on lattices with increasing levels of geometric spin frustration, where tuning the incident energy near the Mott gap reveals strong chiral spin excitation intensity. In conclusion, increased spin frustration in the Mott insulator results in an overall softening of the Raman 𝐴 2⁢𝑔 response, indicating a tendency toward low energy chiral-chiral fluctuations in Mott insulators with magnetic frustration and proximity to chiral spin liquid states that can potentially be tuned by external perturbations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quasiparticle excitations in a one-dimensional interacting topological insulator: Application for dopant-based quantum simulation

Here, we study the effects of electron-electron interactions on the charge excitation spectrum of the spinful Su-Schrieffer-Heeger (SSH) model, a prototype of a one-dimensional bulk obstructed topological insulator. In view of recent progress in the fabrication of dopant-based quantum simulators we focus on experimentally detectable signatures of interacting topology in finite lattices. To this end we use Lanczos-based exact diagonalization to calculate the single-particle spectral function in real space which generalizes the local density of states to interacting systems. Its spatial and spectral resolution allows for the direct investigation and identification of edge states. By studying the noninteracting limit, we demonstrate that the topological in-gap states on the boundary are robust against both finite-size effects as well as random bond and onsite disorder which suggests the feasibility of simulating the SSH model in engineered dopant arrays in silicon. While edge excitations become zero-energy spinlike for any finite interaction strength, our analysis of the spectral function shows that the single-particle charge excitations are gapped out on the boundary. Despite the loss of topological protection, we find that these edge excitations are quasiparticlelike as long as they remain within the bulk gap. Above a critical interaction strength of U c ≈5t these quasiparticles on the boundary lose their coherence which is explained by the merging of edge and bulk states. This is in contrast to the many-body edge excitations which survive the limit of strong coupling, as established in the literature. Our findings show that for moderate repulsive interactions the nontrivial phase of the interacting SSH model can be detected through remnant signatures of topological single-particle states using single-particle local measurement techniques such as scanning tunneling spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Estimation of biquadratic and bicubic Heisenberg effective couplings from multiorbital Hubbard models

Abstract We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U , from weak to strong coupling, and multiple ratios of the Hund coupling J H to U . The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U , we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S = 1 Heisenberg model where by symmetry both nearest-neighbor ( S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2 . For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with ( S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1 , J 2 , and J 3 . The strength of these coupling constants in the Heisenberg models depend on the U , J H , and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2 / J 1 and J 3 / J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2 / J 1 is less than 0.4 and J 3 / J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2 / J 1 and J 3 / J 1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Berry Phase and Model Wave Function in the Half-Filled Landau Level

We construct model wave functions for the half-filled Landau level parametrized by “composite fermion occupation-number configurations” in a two-dimensional momentum space, which correspond to a Fermi sea with particle-hole excitations. When these correspond to a weakly excited Fermi sea, they have a large overlap with wave functions obtained by the exact diagonalization of lowest-Landau-level electrons interacting with a Coulomb interaction, allowing exact states to be identified with quasiparticle configurations. We then formulate a many-body version of the single-particle Berry phase for adiabatic transport of a single quasiparticle around a path in momentum space, and evaluate it using a sequence of exact eigenstates in which a single quasiparticle moves incrementally. In this formulation the standard free-particle construction in terms of the overlap between “periodic parts of successive Bloch wave functions” is reinterpreted as the matrix element of a “momentum boost” operator between the full Bloch states, which becomes the matrix elements of a Girvin-MacDonald-Platzman density operator in the many-body context. Here, this allows the computation of the Berry phase for the transport of a single composite fermion around the Fermi surface. In addition to a phase contributed by the density operator, we find a phase of exactly π for this process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Beyond quantum cluster theories: multiscale approaches for strongly correlated systems

The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. Furthermore, we will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Localization properties and high-fidelity state transfer in hopping models with correlated disorder

We investigate a tight-binding chain featuring diagonal and off-diagonal disorder, these being modeled through the long-range-correlated fractional Brownian motion. Particularly, by employing exact diagonalization methods, we evaluate how the eigenstate spectrum of the system and how its related single-particle dynamics respond to both competing sources of disorder. Moreover, we report the possibility of carrying out efficient end-to-end quantum-state transfer protocols even in the presence of such generalized disorder due to the appearance of extended states around the middle of the band in the limit of strong correlations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2019

Phonon localization in binary alloys with diagonal and off-diagonal disorder: A cluster Green's function approach

We report the development and application of a method for carrying out computational investigations of the effects of mass and force-constant (FC) disorder on phonon spectra. The method is based on the recently developed typical medium dynamical cluster approach, which is a Green's function approach. Excellent quantitative agreement with previous exact diagonalization results establishes the veracity of the method. Application of the method to a model system of binary mass and an FC-disordered system leads to several findings. A narrow resonance, significantly below the van Hove singularity, that has been termed the boson peak, is seen to emerge for low soft particle concentrations. We show, using the typical phonon spectrum, that the states constituting the boson peak cross over from being completely localized to being extended as a function of increasing soft particle concentration. In general, an interplay of mass and FC disorder is found to be cooperative in nature, enhancing phonon localization over all frequencies. However, for a certain range of frequencies, and depending on material parameters, FC disorder can delocalize the states that were localized by mass disorder, and vice versa. Modeling vacancies as weakly bonded sites with vanishing mass, we find that vacancies, even at very low concentrations, are extremely effective in localizing phonons. Thus, inducing vacancies is proposed as a promising route for efficient thermoelectrics. Lastly, we use model parameters corresponding to the alloy system, Ni 1- x Pt x , and we show that mass disorder alone is insufficient to explain the pseudogap in the phonon spectrum; the concomitant presence of FC disorder is necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A j eff = 1/2 Kitaev material on the triangular lattice: the case of NaRuO 2

Motivated by recent reports of a quantum-disordered ground state in the triangular lattice compound NaRuO 2 , we derive a j eff = 1/2 magnetic model for this system by means of first-principles calculations. The pseudospin Hamiltonian is dominated by bond-dependent off-diagonal Γ interactions, complemented by a ferromagnetic Heisenberg exchange and a notably antiferromagnetic Kitaev term. In addition to bilinear interactions, we find a sizable four-spin ring exchange contribution with a strongly anisotropic character, which has been so far overlooked when modeling Kitaev materials. The analysis of the magnetic model, based on the minimization of the classical energy and exact diagonalization of the quantum Hamiltonian, points toward the existence of a rather robust easy-plane ferromagnetic order, which cannot be easily destabilized by physically relevant perturbations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2022

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately, and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence. Here, we adapt pseudodiagonalization to finite-temperature and metallic systems, where partially-occupied orbitals must be individually resolved with some accuracy. We apply pseudodiagonalization to the subspace eigenvalue problem that arises in Chebyshev-filtered subspace iteration. In tests on metallic and other systems for a range of temperatures, we show that pseudodiagonalization achieves similar rates of SCF convergence to exact diagonalization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗