Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electron-mediated entanglement of two distant macroscopic ferromagnets within a nonequilibrium spintronic device

Using the nascent concept of quantum spin-transfer torque, we demonstrate that a current pulse can be harnessed to entangle quantum localized spins of two spatially separated ferromagnets (FMs) which are initially unentangled. The envisaged setup is composed of a spin-polarizer (FM p ) and a spin-analyzer (FM a ) FM layers separated by a normal metal (NM) spacer. The injection of a current pulse into the device leads to a time-dependent superposition of many-body states characterized by a high degree of entanglement between the spin degrees of freedom of the two distant FM layers. The nonequilibrium dynamics are due to the transfer of spin angular momentum from itinerant electrons to the localized spins via a quantum spin-torque mechanism that remains active even for collinear but antiparallel arrangements of the FM p and FM a magnetizations (a situation in which the conventional spin torque is absent). We quantify the mixed-state entanglement generated between the FM layers by tracking the time evolution of the full density matrix and analyzing the build-up of the mutual logarithmic negativity over time. The effect of decoherence and dissipation in the FM layers due to coupling to bosonic baths at finite temperature, the use of multielectron current pulses, and the dependence on the number of spins are also considered in an effort to ascertain the robustness of our predictions under realistic conditions. Lastly, we propose a “current-pump–x-ray-probe” scheme, utilizing ultrafast x-ray spectroscopy, that can witness nonequilibrium and transient entanglement of the FM layers by extracting its time-dependent quantum Fisher information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

Four-Spin Terms and the Origin of the Chiral Spin Liquid in Mott Insulators on the Triangular Lattice

At strong repulsion, the triangular-lattice Hubbard model is described by s=1/2 spins with nearest-neighbor antiferromagnetic Heisenberg interactions and exhibits conventional 120° order. Using the infinite density matrix renormalization group and exact diagonalization, we study the effect of the additional four-spin interactions naturally generated from the underlying Mott-insulator physics of electrons as the repulsion decreases. Although these interactions have historically been connected with a gapless ground state with emergent spinon Fermi surface, we find that, at physically relevant parameters, they stabilize a chiral spin liquid (CSL) of Kalmeyer-Laughlin (KL) type, clarifying observations in recent studies of the Hubbard model. Finally, we present a self-consistent solution based on a mean-field rewriting of the interaction to obtain a Hamiltonian with similarities to the parent Hamiltonian of the KL state, providing a physical understanding for the origin of the CSL.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Spin-1 Kitaev-Heisenberg model on a honeycomb lattice

We study the Kitaev-Heisenberg model with spin-1 local degree of freedom on a honeycomb lattice numerically by the infinite density matrix renormalization group method on a cylinder geometry. By tuning the relative value of the Kitaev and Heisenberg exchange couplings, we obtain the phase diagram with two spin liquid phases and four symmetry-broken phases. Here, we identify that the spin liquid phases are gapless by calculating the central charge at the pure Kitaev points without the Heisenberg interactions. Comparing to its spin-1/2 counterpart, the position and number of gapless modes of the spin-1 case are quite different. Due to the approximate Z 2 local conservations, the expectation value of the Wilson loop operator measuring the flux of each plaquette stays near to 1, and the static spin-spin correlations remain short range in the entire spin liquid phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Natural orbitals for the ab initio no-core configuration interaction approach

Ab initio no-core configuration interaction (NCCI) calculations for the nuclear many-body problem have traditionally relied upon an antisymmetrized product (Slater determinant) basis built from harmonic oscillator orbitals. The accuracy of such calculations is limited by the finite dimensions which are computationally feasible for the truncated many-body space. We therefore seek to improve the accuracy obtained for a given basis size by optimizing the choice of single-particle orbitals. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which maximizes the occupation of low-lying orbitals, thus accelerating convergence in a configuration-interaction basis, while also possibly providing physical insight into the single-particle structure of the many-body wave function. We describe the implementation of natural orbitals in the NCCI framework and examine the nature of the natural orbitals thus obtained, the properties of the resulting many-body wave functions, and the convergence of observables. After taking 3 He as an illustrative testbed, we explore aspects of NCCI calculations with natural orbitals for the ground state of the p-shell neutron halo nucleus 6 He .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Block orbital-selective Mott insulators: A spin excitation analysis

We present a comprehensive study of the spin excitations—as measured by the dynamical spin structure factor S(q,ω)—of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2×2 particular case was already found experimentally in the ladder material BaFe2Se3 that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of S(q,ω) using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard U strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Composite Fermi Liquid at Zero Magnetic Field in Twisted MoTe 2

Here, the pursuit of exotic phases of matter outside of the extreme conditions of a quantizing magnetic field is a long-standing quest of solid state physics. Recent experiments have observed spontaneous valley polarization and fractional Chern insulators in zero magnetic field in twisted bilayers of MoTe 2 , at partial filling of the topological valence band (ν =-2/3 and -3/5). We study the topological valence band at half filling, using exact diagonalization and density matrix renormalization group calculations. We discover a composite Fermi liquid (CFL) phase even at zero magnetic field that covers a large portion of the phase diagram near twist angle ~3.6°. The CFL is a non-Fermi liquid phase with metallic behavior despite the absence of Landau quasiparticles. We discuss experimental implications including the competition between the CFL and a Fermi liquid, which can be tuned with a displacement field. The topological valence band has excellent quantum geometry over a wide range of twist angles and a small bandwidth that is, remarkably, reduced by interactions. These key properties stabilize the exotic zero field quantum Hall phases. Finally, we present an optical signature involving “extinguished” optical responses that detects Chern bands with ideal quantum geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Finite-size subthermal regime in disordered SU ( N ) -symmetric Heisenberg chains

SU⁡(N) symmetry is incompatible with the many-body localized (MBL) phase, even when strong disorder is present. However, recent studies have shown that finite-size SU⁡(2) systems exhibit nonergodic, subthermal behavior, characterized by the breakdown of the eigenstate thermalization hypothesis, and by the excited eigenstates entanglement entropy that is intermediate between area and volume law. In this paper, we extend previous studies of the SU⁡(2)-symmetric disordered Heisenberg model to larger systems, using the time-dependent density matrix renormalization group (tDMRG) method. We simulate quench dynamics from weakly entangled initial states up to long times, finding robust subthermal behavior at stronger disorder. Although we find an increased tendency towards thermalization at larger system sizes, the subthermal regime persists at intermediate time scales, nevertheless, and therefore should be accessible experimentally. At weaker disorder, we observe signatures of thermalization; however, entanglement entropy exhibits slow sublinear growth, in contrast to conventional thermalizing systems. Furthermore, we study dynamics of the SU⁡(3)-symmetric disordered Heisenberg model. Similarly, strong disorder drives the system into subthermal regime, albeit thermalizing phase is broader compared to the SU⁡(2) case. Finally, our findings demonstrate the robustness of the subthermal regime in spin chains with non-Abelian continuous symmetry, and are consistent with eventual thermalization at large system sizes and long time scales, suggested by previous studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Spin dynamics of the generalized quantum spin compass chain

Here we calculate the dynamical spin structure factor of the generalized spin-1/2 compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb 2 O 6 . It features bond-dependent interactions and interpolates between an Ising chain and a one-dimensional variant of Kitaev's honeycomb spin model. The structure factor, in turn, is found to interpolate from gapped and nondispersive in the Ising limit to gapless with nontrivial continua in the Kitaev limit. In particular, the component of the structure factor perpendicular to the Ising directions changes abruptly at the Kitaev point into a dispersionless continuum due to the emergence of an extensive ground-state degeneracy. We show this continuum is consistent with analytical Jordan-Wigner results. We also discuss implications for future inelastic scattering experiments and applications to materials, particularly CoNb 2 O 6 .

1-dimensional spin chains↗
DOE OSTI2023

Quantum Phase Diagram and Spontaneously Emergent Topological Chiral Superconductivity in Doped Triangular-Lattice Mott Insulators

The topological superconducting state is a highly sought-after quantum state hosting topological order and Majorana excitations. Here, in this Letter, we explore the mechanism to realize the topological superconductivity (TSC) in the doped Mott insulators with time-reversal symmetry (TRS). Through large-scale density matrix renormalization group study of an extended triangular-lattice t–J model on the six- and eight-leg cylinders, we identify a d + id-wave chiral TSC with spontaneous TRS breaking, which is characterized by a Chern number C=2 and quasi-long-range superconducting order. We map out the quantum phase diagram with by tuning the next-nearest-neighbor (NNN) electron hopping and spin interaction. In the weaker NNN-coupling regime, we identify a pseudogaplike phase with a charge stripe order coexisting with fluctuating superconductivity, which can be tuned into d -wave superconductivity by increasing the doping level and system width. The TSC emerges in the intermediate-coupling regime, which has a transition to a d-wave superconducting phase with larger NNN couplings. The emergence of the TSC is driven by geometrical frustrations and hole dynamics which suppress spin correlation and charge order, leading to a topological quantum phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Tracing the Mott-Hubbard transition in one-dimensional Hubbard models without Umklapp scattering

Here we apply the density-matrix renormalization group (DMRG) method to a one-dimensional Hubbard model that lacks Umklapp scattering and thus provides an ideal case to study the Mott-Hubbard transition analytically and numerically. The model has a linear dispersion and displays a metal-to-insulator transition when the Hubbard interaction U equals the bandwidth U c = W, where the single-particle gap opens linearly, Δ(U ≥ W) = U - W. The simple nature of the elementary excitations permits us to determine numerically with high accuracy the critical interaction strength and the gap function in the thermodynamic limit. The jump discontinuity of the momentum distribution $n_k$ at the Fermi wave number $k_F$ = 0 cannot be used to locate accurately U c from finite-size systems. However, the slope of $n_k$ at the band edges $k_B$ = ±π, reveals the formation of a single-particle bound state which can be used to determine U c reliably from $n_k$ using accurate finite-size data.

1-dimensional systems↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗
DOE OSTI2023

Characterizing random-singlet state in two-dimensional frustrated quantum magnets and implications for the double perovskite Sr 2 CuTe 1 – x W x O 6

Motivated by the experimental observation of a nonmagnetic phase in compounds with frustration and disorder, we study the ground state of a spin-1/2 square-lattice Heisenberg model with randomly distributed nearest-neighbor J 1 and next-nearest-neighbor J 2 couplings. By using the density matrix renormalization group (DMRG) calculation on a cylinder system with a circumference of up to ten lattice sites, we identify a disordered phase between the Néel and stripe magnetic phase with growing J 2 /J 1 in the presence of strong bond randomness. The vanished spin-freezing parameter indicates the absence of spin-glass order. The large-scale DMRG results unveil the size-scaling behaviors of the spin-freezing parameter, the power-law decay of the average spin correlation, and the exponential decay of the typical spin correlation, which all agree with the corresponding behavior in the one-dimensional random-singlet (RS) state and characterize the RS nature of this disordered phase. The DMRG simulation also provides insights and opportunities for characterizing a class of nonmagnetic states in two-dimensional frustrated magnets with disorder. Here, we also compare with existing experiments and suggest more measurements for understanding the spin-liquid-like behaviors in the double perovskite Sr 2 CuTe 1–x W x O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Non-Hermitian approach for quantum plasmonics

Herein, we examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. Furthermore, while such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful. Published under license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Phase Diagram of the Easy-Axis Triangular-Lattice 𝐽 1 −𝐽 2 Model

Our work studies the largely unexplored phase diagram of the S = 1/2 easy-axis triangular-lattice J 1 -J 2 model, motivated by recent interest in rare-earth and transition-metal compounds that exhibit strong quantum fluctuations and broad spin excitation continua—hallmarks of spin-liquid behavior. Combining density-matrix renormalization group simulations with analytical insights, we present a comprehensive study of the model that describes these materials. Our results provide compelling evidence for a robust spin-liquid phase stabilized between the exotic supersolid Y phase and the collinear stripe phase, remarkably persistent even in the presence of symmetry-breaking anisotropy. Additionally, we analyze the supersolid Y phase, offering a quantitative characterization of its order parameters and clarifying the surprising absence of a ferromagnetic moment; both features relevant to recent experiments on Ising-like triangular-lattice magnets. Altogether, our work provides the necessary framework and important theoretical guidance to the ongoing searches of the spin liquid and other exotic states in the rare-earth and transition-metal triangular-lattice compounds, as well as for advancing the understanding of anisotropic-exchange magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Polariton spectra under the collective coupling regime. II. 2D non-linear spectra

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton–polariton dynamics described by the Holstein–Tavis–Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (⁠$\mathscr{L}$-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton–polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from $\mathscr{O}$($\mathscr{K}$ 2 )$\mathscr{O}$(T 3 ) to $\mathscr{O}$($\mathscr{K}$)$\mathscr{O}$(T 0 ) for the 2DES simulation, where $\mathscr{K}$ is the number of states and T is the number of time steps of propagation. As a result, we further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways.

2D non-linear spectra↗
DOE OSTI2025

Topological phase diagram of generalized Su-Schrieffer-Heeger models with interactions

We investigate interacting Su-Schrieffer-Heeger (SSH) chains with two- and three-site unit cells using density matrix renormalization group (DMRG) simulations. By selecting appropriate filling fractions and sweeping across interaction strength 𝐽 𝑧 and dimerization 𝛿, we map out their phase diagrams and identify transition lines via entanglement entropy and magnetization measurements. In the two-site model, we observe the emergence of an interaction-induced antiferromagnetic intermediate phase between the topologically trivial and nontrivial regimes, as well as a critical region at negative 𝐽 𝑧 with suppressed magnetization and finite-size scaling of entanglement entropy. In contrast, the three-site model lacks an intermediate phase and exhibits asymmetric edge localization and antiferromagnetic ordering in both positive and negative 𝐽 𝑧 regimes. We further examine the response of edge states to Ising perturbations. In the two-site model, zero-energy edge modes are topologically protected and remain robust up to a finite interaction strength. However, in the three-site model, where the edge states reside at finite energy, this protection breaks down. Despite this, the edge-localized nature of these states survives in the form of polarized modes whose spatial profiles reflect the noninteracting limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗