Search NASA⌕ Search

SEARCH · Search NASA

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 775 records · Page 43

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗

Magnetic field-induced intermediate quantum spin liquid with a spinon Fermi surface

The Kitaev model with an applied magnetic field in the H ∥ [ 111 ] direction shows two transitions: from a nonabelian gapped quantum spin liquid (QSL) to a gapless QSL at H c 1 ≃ 0.2 K and a second transition at a higher field H c 2 ≃ 0.35 K to a gapped partially polarized phase, where K is the strength of the Kitaev exchange interaction. We identify the intermediate phase to be a gapless U(1) QSL and determine the spin structure function S ( k ) and the Fermi surface ϵ F S ( k ) of the gapless spinons using the density matrix renormalization group (DMRG) method for large honeycomb clusters. Further calculations of static spin-spin correlations, magnetization, spin susceptibility, and finite temperature-specific heat and entropy corroborate the gapped and gapless nature of the different field-dependent phases. In the intermediate phase, the spin-spin correlations decay as a power law with distance, indicative of a gapless phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stripe order enhanced superconductivity in the Hubbard model

Unidirectional (“stripe”) charge density wave order has now been established as a ubiquitous feature in the phase diagram of the cuprate high-temperature superconductors, where it generally competes with superconductivity. Nonetheless, on theoretical grounds it has been conjectured that stripe order (or other forms of “optimal” inhomogeneity) may play an essential positive role in the mechanism of high-temperature superconductivity. Here, we report density matrix renormalization group studies of the Hubbard model on long four- and six-leg cylinders, where the hopping matrix elements transverse to the long direction are periodically modulated—mimicking the effect of putative period 2 stripe order. We find that even modest amplitude modulations can enhance the long-distance superconducting correlations by many orders of magnitude and drive the system into a phase with a substantial spin gap and superconducting quasi–long-range order with a Luttinger exponent, K sc ~ 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ground-state phase diagram of the t-t ′ -J model

Significance Theoretical efforts to understand the superconducting cuprates have focused on a few simplified but challenging models: notably the Hubbard and the t - t ′ - J models. Using density matrix renormalization group (DMRG) simulations on width-8 cylinders, we study the ground-state phase diagram of the t - t ′- J model. The next nearest-neighbor hopping t ′ allows one to distinguish between hole ( t ′ < 0 ) and electron doping ( t ′ > 0 ). We find for the magnetic, charge, and single-particle properties that this model roughly captures the differences between the electron- and hole-doped cuprates. However, the superconducting properties differ. For hole doping ( t ′ < 0 ), there is an absence of superconductivity, while for electron doping ( t ′ > 0), we find d -wave superconductivity coexisting with antiferromagnetism and a secondary triplet p -wave pairing component.

Jiang, Shengtao↗

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗

Orbital optimisation in the perfect pairing hierarchy: applications to full-valence calculations on linear polyacenes

Here, we describe the implementation of orbital optimisation for the models in the perfect pairing hierarchy. Orbital optimisation, which is generally necessary to obtain reliable results, is pursued at perfect pairing (PP) and perfect quadruples (PQ) levels of theory for applications on linear polyacenes, which are believed to exhibit strong correlation in the π space. While local minima and σ-π symmetry breaking solutions were found for PP orbitals, no such problems were encountered for PQ orbitals. The PQ orbitals are used for single-point calculations at PP, PQ and perfect hextuples (PH) levels of theory, both only in the π subspace, as well as in the full σπ valence space. It is numerically demonstrated that the inclusion of single excitations is necessary also when optimised orbitals are used. PH is found to yield good agreement with previously published density matrix renormalisation group data in the π space, capturing over 95% of the correlation energy. Full-valence calculations made possible by our novel, efficient code reveal that strong correlations are weaker when larger basis sets or active spaces are employed than in previous calculations. The largest full-valence PH calculations presented correspond to a (192e,192o) problem.

74 ATOMIC AND MOLECULAR PHYSICS↗

Ab initio no-core solutions for 6 Li

Here, we solve for properties of 6 Li in the ab initio no-core full configuration (NCFC) approach and we separately solve for its ground state and ${J}^{\pi }={2}_{2}^{+}$ resonance with the Gamow shell model (GSM) in the Berggren basis. We employ both the JISP16 and chiral ${\mathrm{NNLO}}_{\mathrm{opt}}$ realistic nucleon–nucleon interactions and investigate the ground state energy, excitation energies, point proton root mean square (rms) radius and a suite of electroweak observables. We also extend and test methods to extrapolate the ground state energy, point proton rms radius, and electric quadrupole moment. We attain improved estimates of these observables in the NCFC approach by using basis spaces up through ${N}_{\max }=18$ that enable more definitive comparisons with experiment. Using the density matrix renormalization group approach with the JISP16 interaction, we find that we can significantly improve the convergence of the GSM treatment of the 6Li ground state and ${J}^{\pi }={2}_{2}^{+}$ resonance by adopting a natural orbital single-particle basis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Observation of quantum entanglement in top quark pair production in proton–proton collisions at $\sqrt{s}$ = 13 TeV

Entanglement is an intrinsic property of quantum mechanics and is predicted to be exhibited in the particles produced at the Large Hadron Collider. A measurement of the extent of entanglement in top quark-antiquark ($t\bar{t}$) events produced in proton–proton collisions at a center-of-mass energy of 13 TeV is performed with the data recorded by the CMS experiment at the CERN LHC in 2016, and corresponding to an integrated luminosity of 36.3 fb −1 . The events are selected based on the presence of two leptons with opposite charges and high transverse momentum. An entanglement-sensitive observable D is derived from the top quark spin-dependent parts of the $t\bar{t}$ production density matrix and measured in the region of the $t\bar{t}$ production threshold. Values of $D \lt -1/3$ are evidence of entanglement and D is observed (expected) to be $-0.480^{+0.026}_{-0.029}$ ($-0.467^{+0.026}_{-0.029}$) at the parton level. With an observed significance of 5.1 standard deviations with respect to the non-entangled hypothesis, this provides observation of quantum mechanical entanglement within $t\bar{t}$ pairs in this phase space. This measurement provides a new probe of quantum mechanics at the highest energies ever produced.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superconductor to exciton condensate transition in a model copper-oxide material

Abstract Superconductivity and exciton condensation are fundamental phenomena in condensed matter physics, associated with the condensation of electron–electron and electron–hole pairs, respectively, into coherent quantum states. In this study, we present evidence of a superconductor to exciton condensate transition within the context of the three-band Hubbard model of copper-oxide-like materials. As the electron–electron repulsion increases, the superconducting phase is superseded by exciton condensation. In support of theoretical predictions—not yet realized experimentally—we observe the coexistence of the two condensates in the vicinity of the transition where the quantum states become a superposition of electron–electron and electron–hole condensates. Coexistence is rigorously computed from large eigenvalues and their eigenvectors in both the two-electron reduced density matrix (2-RDM) and the particle-hole RDM, which we obtain from a direct variational ground-state energy minimization with respect to the 2-RDM by semidefinite programming. We further discern that adjacent d orbitals and intervening p orbitals facilitate electron–electron pairing between copper orbitals, thereby supporting the superexchange mechanism for superconductivity. These observations suggest the feasibility of witnessing a superconductor to exciton condensate transition in copper-oxide analogs, bearing significant implications for identifying materials conducive to efficient transport processes.

Schouten, Anna O. (ORCID:0000000264255219)↗

A limit on the total lepton number in the Universe from BBN and the CMB

At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Eigenstate thermalization hypothesis in conformal field theory

In this work, we investigate the eigenstate thermalization hypothesis (ETH) in d + 1 dimensional conformal field theories by studying the reduced density matrices in energy eigenstates. We show that if the local probes of the finitely excited primary eigenstates satisfy ETH, then any finite energy observable with support on a subsystem of finite size satisfies ETH. In two dimensions, we discover that if ETH holds locally, the finite size reduced density matrix of states created by heavy primary operators is well-approximated by a projection to the Virasoro identity block.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photoproduction of the Λ(1520) hyperon with a 9 GeV photon beam at GlueX

The GlueX experiment is located at the Thomas Jefferson National Accelerator Facility (JLab) in Newport News, VA, USA. It features a hermetic 4π detector with excellent tracking and calorimetry capabilities. Its 9 GeV linearly polarized photon beam is produced from the 12 GeV electron beam, delivered by JLab’s Continuous Electron Beam Accelerator Facility (CEBAF), via bremsstrahlung on a thin diamond and is incident on a LH2 target. GlueX recently finished its first data taking period and published first results. The main goal of GlueX is to measure gluonic excitations of mesons. These so-called hybrid or exotic mesons are predicted by Quantum Chromodynamics (QCD) but haven’t been experimentally confirmed yet. They can have quantum numbers not accessible by ordinary quark-antiquark pairs which helps in identifying them using partial wave analysis techniques. The search for exotic mesons requires a very good understanding of photoproduction processes in a wide range of final states, one of them being pK + K - which contains many meson and baryon reactions. The Λ(1520) is a prominent hyperon resonance in this final state and is the subject of this presentation. This talk will give an introduction to the GlueX experiment and show preliminary results for the photoproduction of the Λ(1520) hyperon. The measurement of important observables like the photon beam asymmetry and spin-density matrix elements will be discussed and an outlook to possible measurements of further hyperon states in the pK + K - final state will be given.

Pauli, Peter↗

Ultrafast charge dynamics and photoluminescence in bilayer MoS 2

Our examination of the interplay of ultrafast charge dynamics and electron–phonon interaction in the AA' stacked bilayer MoS 2 provides a microscopic basis for understanding the features (two peaks) in the emission spectrum. We demonstrate that while the initial accumulation of excited charge occurs at and near the Q point of the two-dimensional Brillioun zone, emission takes place predominantly through two pathways: direct charge recombination at the K point and indirect phonon-assisted recombination of electrons at the K valley and holes at the Γ hill of the Brillouin zone. Analysis of the wave vector dependencies of the electron–phonon interaction traces the higher energy peak to phonon-assisted relaxation of the excited electrons from the Q to the K valley in the conduction band. Our results thus reveal the importance of ultrafast charge dynamics in understanding photoemissive properties of a few-layer transition-metal dichalcogenide. These calculations are based on time dependent density functional theory in the density matrix formulation.

Materials Science↗

Mutual information-assisted adaptive variational quantum eigensolver

Adaptive construction of ansatz circuits offers a promising route towards applicable variational quantum eigensolvers on near-term quantum hardware. Those algorithms aim to build up optimal circuits for a certain problem and ansatz circuits are adaptively constructed by selecting and adding entanglers from a predefined pool. In this work, we propose a way to construct entangler pools with reduced size by leveraging classical algorithms. Our method uses mutual information between the qubits in classically approximated ground state to rank and screen the entanglers. The density matrix renormalization group method is employed for classical precomputation in this work. We corroborate our method numerically on small molecules. Our numerical experiments show that a reduced entangler pool with a small portion of the original entangler pool can achieve same numerical accuracy. Here, we believe that our method paves a new way for adaptive construction of ansatz circuits for variational quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating noise on a quantum processor: interactions between a qubit and resonant two-level system bath

Material defects fundamentally limit the coherence times of superconducting qubits, and manufacturing completely defect-free devices is not yet possible. Therefore, understanding the interactions between defects and a qubit in a real quantum processor design is essential. We build a model that incorporates the standard tunneling model, the electric field distributions in the qubit, and open quantum system dynamics, and draws from the current understanding of two-level system (TLS) theory. Specifically, we start with one million TLSs distributed on the surface of a qubit and pick the 200 systems that are most strongly coupled to the qubit. We then perform a full Lindbladian simulation that explicitly includes the coherent coupling between the qubit and the TLS bath to model the time dependent density matrix of resonant TLS defects and the qubit. We find that the 200 most strongly coupled TLSs can accurately describe the qubit energy relaxation time. This work confirms that resonant TLSs located in areas where the electric field is strong can significantly affect the qubit relaxation time, even if they are located far from the Josephson junction (JJ). Similarly, a strongly-coupled resonant TLS located in the JJ does not guarantee a reduced qubit relaxation time if a more strongly coupled TLS is far from the JJ. In addition to the coupling strengths between TLSs and the qubit, the model predicts that the geometry of the device and the TLS relaxation time play a significant role in qubit dynamics. Our work can provide guidance for future quantum processor designs with improved qubit coherence times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗