Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Lower bounds on entanglement entropy without twin copy

We discuss the possibility of estimating experimentally the von Neumann entanglement entropy S A v N of a symmetric bipartite quantum system A B by using the basic measurement counts (bitstrings) for a single copy of a prepared state. Using exact diagonalization and analog simulations performed with the publicly available QuEra facilities for chains and ladders of Rydberg atoms, we calculate the Shannon entropy S A B X associated with the bitstrings of adiabatically prepared ground states and the reduced entropies S A X and S B X obtained from the marginal probabilities in A and B . We then calculate the classical mutual information I A B X = S A X + S B X − S A B X , which is a lower bound on S A v N . We show that for a broad range of lattice spacing and detuning, I A B X is typically 20% below S A v N in regions where S A v N is large and a less close bound in regions where S A v N is low. We argue that this use of the easily available bitstrings provides a robust and efficient way to explore empirically the phase diagram of qubit-based quantum simulators and identify critical regions. Published by the American Physical Society 2025

Meurice, Yannick (ORCID:0000000209959694)↗
DOE OSTI2025

Superspin renormalization and slow relaxation in random spin systems

We develop an excited-state real-space renormalization group (RSRG-X) formalism to describe the dynamics of conserved densities in randomly interacting spin-12 systems. Our formalism is suitable for systems with U(1) and Z2 symmetries, and we apply it to chains of randomly positioned spins with dipolar XX+YY interactions, as arise in Rydberg quantum simulators and other platforms. The formalism generates a sequence of effective Hamiltonians that provide approximate descriptions for dynamics on successively smaller energy scales. These effective Hamiltonians involve “superspins”: two-level collective degrees of freedom constructed from (anti)aligned microscopic spins. Conserved densities can then be understood as relaxing via coherent collective spin flips. For the well-studied simpler case of randomly interacting nearest-neighbor XX+YY chains, the superspins reduce to single spins. Our formalism also leads to a numerical method capable of simulating the dynamics up to an otherwise inaccessible combination of large system size and late time. Focusing on disorder-averaged infinite-temperature autocorrelation functions, in particular the spin survival probability Sp¯(t), we demonstrate quantitative agreement between our algorithm and exact diagonalization (ED) at low but nonzero frequencies. Such agreement holds for chains with nearest-neighbor, next-nearest-neighbor, and long-range dipolar interactions. Our results indicate decay of Sp¯(t) slower than any power law and feature no significant deviation from the ∼1/ln2(t) asymptote expected from the infinite-randomness fixed-point of the nearest-neighbor model. We also apply the RSRG-X formalism to two-dimensional long-range systems of moderate size and find slow late-time decay of Sp¯(t).

Zhao, Yi J↗
DOE OSTI2026

Spectroscopic Demarcation of Emergent Photons and Spinons in a Dipolar-Octupolar Quantum Spin Liquid

The identification of fractionalized excitations in quantum spin liquids (QSLs) remains a central challenge in condensed matter physics. In dipolar-octupolar (DO) pyrochlores, such as Ce 2⁢ Zr 2 ⁢O 7 , the candidate 𝜋-flux quantum spin ice (QSI) state is predicted to host both gapless emergent photons and a continuum of spinons. However, resolving these modes at zero field is complicated by their spectral overlap and the presence of nonmagnetic scattering near zero energy. Here, in this study, we report neutron scattering experiments on Ce 2 ⁢Zr 2 ⁢O 7 under a magnetic field along the [1, 1, 1] direction. In contrast to previous unpolarized studies at zero field that relied on high-temperature subtraction, we use a same-temperature high-field subtraction protocol to isolate the photon and spinon modes. Leveraging the selective coupling of the magnetic field to the dipolar degrees of freedom, we demonstrate the spectroscopic demarcation of these excitations. We observe that weak fields (≈ 0.15 T) suppress the low-energy photon weight while leaving the high-energy spinon continuum robust, albeit hardened. Our results, supported by gauge mean-field theory and exact diagonalization calculations, provide strong evidence for the 𝜋-flux QSI state and introduce a powerful field-tuning protocol for investigating DO-QSLs.

Gao, Bin [Rice University, Houston, TX (United Sta↗
DOE OSTI2023

A path integral ground state Monte Carlo algorithm for entanglement of lattice bosons

A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dimer description of the SU(4) antiferromagnet on the triangular lattice

In systems with many local degrees of freedom, high-symmetry points in the phase diagram can provide an important starting point for the investigation of their properties throughout the phase diagram. In systems with both spin and orbital (or valley) degrees of freedom such a starting point gives rise to SU(4)-symmetric models. Here we consider SU(4)-symmetric "spin" models, corresponding to Mott phases at half-filling, i.e. the six-dimensional representation of SU(4). This may be relevant to twisted multilayer graphene. In particular, we study the SU(4) antiferromagnetic "Heisenberg" model on the triangular lattice, both in the classical limit and in the quantum regime. Carrying out a numerical study using the density matrix renormalization group (DMRG), we argue that the ground state is non-magnetic. We then derive a dimer expansion of the SU(4) spin model. An exact diagonalization (ED) study of the effective dimer model suggests that the ground state breaks translation invariance, forming a valence bond solid (VBS) with a 12-site unit cell. Finally, we consider the effect of SU(4)-symmetry breaking interactions due to Hund's coupling, and argue for a possible phase transition between a VBS and a magnetically ordered state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Subsystem Rényi entropy of thermal ensembles for SYK-like models

The Sachdev-Ye-Kitaev model is an N N -modes fermionic model with infinite range random interactions. In this work, we study the thermal Rényi entropy for a subsystem of the SYK model using the path-integral formalism in the large- N N limit. The results are consistent with exact diagonalization and can be well approximated by thermal entropy with an effective temperature when subsystem size M\leq N/2 M ≤ N / 2 . We also consider generalizations of the SYK model with quadratic random hopping term or U(1) U ( 1 ) charge conservation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Evolution of Charge-Lattice Dynamics across the Kuramoto Synchronization Phase Diagram of Quantum Tunneling Polarons in Cuprate Superconductors

Because of its sensitivity to the instantaneous structure factor, S(Q,t = 0), Extended X-ray Absorption Fine Structure (EXAFS) is a powerful tool for probing the dynamic structure of condensed matter systems in which the charge and lattice dynamics are coupled. When applied to hole-doped cuprate superconductors, EXAFS has revealed the presence of internal quantum tunneling polarons (IQTPs). An IQTP arises in EXAFS as a two-site distribution for certain Cu–O pairs, which is also duplicated in inelastic scattering but not observed in standard diffraction measurements. The Cu–Sr pair distribution has been found to be highly anharmonic and strongly correlated to both the IQTPs and to superconductivity, as, for example, in YSr 2 Cu 2.75 Mo 0.25 O 7.54 (T c =84 K). In order to describe such nontrivial, anharmonic charge-lattice dynamics, we have proposed a model Hamiltonian for a prototype six-atom cluster, in which two Cu-apical-O IQTPs are charge-transfer bridged through Cu atoms by an O atom in the CuO 2 plane and are anharmonically coupled via a Sr atom. By applying an exact diagonalization procedure to this cluster, we have verified that our model indeed produces an intricate interplay between charge and lattice dynamics. Then, by using the Kuramoto model for the synchronization of coupled quantum oscillators, we have found a first-order phase transition for the IQTPs into a synchronized, phase-locked phase. Most importantly, we have shown that this transition results specifically from the anharmonicity. Finally, we have provided a phase diagram showing the onset of the phase-locking of IQTPs as a function of the charge-lattice and anharmonic couplings in our model. We have found that the charge, initially confined to the apical oxygens, is partially pumped into the CuO 2 plane in the synchronized phase, which suggests a possible connection between the synchronized dynamic structure and high-temperature superconductivity (HTSC) in doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Neural Network Solver for Small Quantum Clusters

Machine learning approaches have recently been applied to the study of various problems in physics. Most of these studies are focused on interpreting the data generated by conventional numerical methods or the data on an existing experimental database. An interesting question is whether it is possible to use a machine learning approach, in particular a neural network, for solving the many-body problem. In this paper, we present a neural network solver for the single impurity Anderson model, the paradigm of an interacting quantum problem in small clusters. We demonstrate that the neural-network-based solver provides quantitative accurate results for the spectral function as compared to the exact diagonalization method. This opens the possibility of utilizing the neural network approach as an impurity solver for other many-body numerical approaches, such as the dynamical mean field theory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determining Ground-State Phase Diagrams on Quantum Computers via a Generalized Application of Adiabatic State Preparation

Quantum phase transitions materialize as level crossings in the ground-state energy when the parameters of the Hamiltonian are varied. The resulting ground-state phase diagrams are straightforward to determine by exact diagonalization on classical computers, but are challenging on quantum computers because of the accuracy needed and the near degeneracy of the competing states close to the level crossings. On the other hand, classical computers are limited to small system sizes, which quantum computers may help overcome. In this work, we use a local adiabatic ramp for state preparation to allow us to directly compute ground-state phase diagrams on a quantum computer via time evolution. This methodology is illustrated by examining the ground states of the XY model with a magnetic field in the z-direction in one dimension. We are able to calculate an accurate phase diagram on both two- and three-site systems using IBM quantum machines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Non-LTE radiative transfer with lambda-acceleration - Convergence properties using exact full and diagonal lambda-operators

We investigate the convergence properties of Lambda-acceleration methods for non-LTE radiative transfer problems in planar and spherical geometry. Matrix elements of the 'exact' A-operator are used to accelerate convergence to a solution in which both the radiative transfer and atomic rate equations are simultaneously satisfied. Convergence properties of two-level and multilevel atomic systems are investigated for methods using: (1) the complete Lambda-operator, and (2) the diagonal of the Lambda-operator. We find that the convergence properties for the method utilizing the complete Lambda-operator are significantly better than those of the diagonal Lambda-operator method, often reducing the number of iterations needed for convergence by a factor of between two and seven. However, the overall computational time required for large scale calculations - that is, those with many atomic levels and spatial zones - is typically a factor of a few larger for the complete Lambda-operator method, suggesting that the approach should be best applied to problems in which convergence is especially difficult.

Macfarlane, J. J.↗
DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2014

Line Interference Effects Using a Refined Robert-Bonamy Formalism: the Test Case of the Isotropic Raman Spectra of Autoperturbed N2

A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N2 for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N2, opening the way to the analysis of more complex molecular systems.

nitrogen↗
DOE OSTI2024

Fingerprints of composite fermion Lambda levels in scanning tunneling microscopy

A composite fermion (CF) is a topological quasiparticle that emerges from a nonperturbative attachment of vortices to electrons in strongly correlated two-dimensional materials. Similar to noninteracting fermions that form Landau levels in a magnetic field, CFs can fill analogous “Lambda” levels, giving rise to the fractional quantum Hall (FQH) effect of electrons. Here, we show that Lambda levels can be directly visualized through the characteristic peak structure in the signal obtained via spectroscopy with scanning tunneling microscopy (STM) on a FQH state. Complementary to transport, which probes the low-energy properties of CFs, we show that high-energy features in STM spectra can be interpreted in terms of Lambda levels. We numerically demonstrate that STM spectra can be accurately modeled using Jain's CF theory. Our results show that STM provides a powerful tool for revealing the anatomy of FQH states and identifying physics beyond the noninteracting CF paradigm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Exact solution of the A{sub 2}{sup (2)} model with non-diagonal boundary reflections

The exact solution of the A{sub 2}{sup (2)} model (or the Izergin–Korepin model) with non-diagonal boundary reflections of the type I is studied by using the off-diagonal Bethe ansatz method. The complete operator product identities of the transfer matrices are obtained with the help of fusion. Based on them and the asymptotic behaviors as well as the values at special points, the inhomogeneous T−Q relations of the system are constructed. The method and the results in this paper can be generalized to the high rank A{sub n}{sup (2)} case directly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Superconductivity in a system of interacting spinful semions

Noninteracting particles obeying certain fractional statistics have been predicted to exhibit superconductivity. We discuss the issue in an attractively interacting system of spinful semions on a lattice by numerically investigating the presence of off-diagonal long-range order at zero temperature. For this purpose we construct a Hubbard model wherein two semions with opposite spin can virtually coincide while maintaining consistency with the fractional braiding statistics. Clear off-diagonal long-range order is seen in the strong coupling limit, consistent with the expectation that a pair of semions obeys Bose statistics. We find that the semion system behaves similarly to a system of fermions with the same attractive Hubbard U interaction for a wide range of U, suggesting that semions also undergo a BCS to BEC crossover as a function of U.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗