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 307 records · Page 17

TeraChem: Accelerating electronic structure and ab initio molecular dynamics with graphical processing units

Developed over the past decade, TeraChem is an electronic structure and ab initio molecular dynamics software package designed from the ground up to leverage graphics processing units (GPUs) to perform large-scale ground and excited state quantum chemistry calculations in the gas and the condensed phase. TeraChem’s speed stems from the reformulation of conventional electronic structure theories in terms of a set of individually optimized high-performance electronic structure operations (e.g., Coulomb and exchange matrix builds, one- and two-particle density matrix builds) and rank-reduction techniques (e.g., tensor hypercontraction). Recent efforts have encapsulated these core operations and provided language-agnostic interfaces. Finally, this greatly increases the accessibility and flexibility of TeraChem as a platform to develop new electronic structure methods on GPUs and provides clear optimization targets for emerging parallel computing architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BCS-BEC crossover in a (t 2g ) 4 excitonic magnet

The condensation of spin-orbit-induced excitons in t 2 g 4 electronic systems is attracting considerable attention. At large Hubbard U , antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ( J eff = 1 ). Furthermore, we show that even at intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ( λ ) in one dimension using the density matrix renormalization group, while in two dimensions we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full λ versus U phase diagrams for both one- and two-dimensional lattices. Our main result is that at intermediate Hubbard U , increasing λ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O ( a ) and O ( 10 a ) in 1 d and 2 d , respectively, with a being the lattice spacing. Further increasing λ , the system eventually crosses over to the BEC limit (with r coh ≪ a ).

1-dimensional systems↗

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗

Optical gain for the interband optical transition in InAsP/InP quantum well wire in the influence of laser field intensity

Intense high frequency laser field induced electronic and optical properties of heavy hole exciton in the InAs{sub 0.8}P{sub 0.2}/InP quantum wire is studied taking into account the geometrical confinement effect. Laser field related exciton binding energies and the optical band gap in the InAs{sub 0.8}P{sub 0.2}/InP quantum well wire are investigated. The optical gain, for the interband optical transition, as a function of photon energy, in the InAs{sub 0.8}P{sub 0.2}/InP quantum wire, is obtained in the presence of intense laser field. The compact density matrix method is employed to obtain the optical gain. The obtained optical gain in group III-V narrow quantum wire can be applied for achieving the preferred telecommunication wavelength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Time-dependent approach to inelastic scattering spectroscopies in and away from equilibrium: Beyond perturbation theory

We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schrödinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, as an example.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Linear Light Amplifier and Amplification of N-Photon States

We consider a linear quantum amplifier consisting of NA two-level atoms and study the problem of amplification of N-photon states. The N-photon states are associated with N-quantum states of the harmonic oscillator. We show that the process of interaction of the electromagnetic field with atoms can be associated with some transformation of the phase space and functions defined on this phase space. We consider the Husimi functions QN(q, p) of N-quantum states of the harmonic oscillator, which are defined on the phase space, investigate transformation of these functions, and find an explicit form of the density matrix of the amplified N-photon state. © 2019, Springer Science+Business Media, LLC, part of Springer Nature.

60 APPLIED LIFE SCIENCES↗

Lorentzian geometry for detecting qubit entanglement

We describe a new approach based on Lorentzian geometry to detect qubit entanglement. The treatment is based physically, on the causal structure of Minkowski spacetime, and mathematically, on a Lorentzian Singular Value Decomposition. A surprising feature is the natural emergence of “Energy conditions” used in Relativity. All states satisfy a “Dominant Energy Condition” (DEC) and separable states satisfy the Strong Energy Condition(SEC), while entangled states violate the SEC. We thus propose a test for two qubit entanglement which is an alternative to the positive partial transpose (PPT) test. This test is based on the partial Lorentz transformation (PLT) on individual qubits. Apart from testing for entanglement, our approach also enables us to construct a separable form for the density matrix in those cases where it exists. Our approach leads to a simple graphical three dimensional representation of the state space which shows the entangled states within the set of all states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-pressure inelastic neutron scattering study of the anisotropic S = 1 spin chain [ Ni ( HF 2 ) ( 3 - Clpyradine ) 4 ] B F 4

[Ni(HF 2 )(3-Clpyradine) 4 ]BF 4 (NBCT) is a one-dimensional, S=1 spin-chain material that shows no long-range magnetic order down to thermometer temperatures of 0.1 K. Previous ambient pressure inelastic neutron scattering experiments identified NBCT to be in the large-D quantum paramagnetic phase of the D/J phase diagram, where D is the axial single-ion anisotropy and J is the intrachain superexchange. In this work, we extend the previous experiments to a hydrostatic pressure of 0.9 GPa. By comparing to density matrix renormalization group calculations, we find D/J increases from 1.5 to 3.2 as pressure increases from 0 GPa to 0.9 GPa, which pushes the system further into the large-D phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multisliced gausslet basis sets for electronic structure

We introduce highly local basis sets for electronic structure which are very efficient for correlation calculations near the complete basis set limit. Our approach is based on gausslets, recently introduced waveletlike smooth orthogonal functions. We adapt the gausslets to particular systems using one-dimensional (1D) coordinate transformations, putting more basis functions near the nuclei, while maintaining orthogonality. Furthermore, three-dimensional basis functions are composed out of products of the 1D functions in an efficient way called multislicing. We demonstrate these bases with both Hartree-Fock and density matrix renormalization group calculations on hydrogen chain systems. With both methods, we can go to higher accuracy in the complete basis set limit than is practical for conventional Gaussian basis sets, with errors near 0.1 mH per atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulation of Nonequilibrium Dynamics on a Quantum Computer

We present a hybrid quantum-classical algorithm for the time evolution of out-of-equilibrium thermal states. Here, the method depends on classically computing a sparse approximation to the density matrix and, then, time-evolving each matrix element via the quantum computer. For this exploratory study, we investigate a time-dependent Ising model with five spins on the Rigetti Forest quantum virtual machine and a one spin system on the Rigetti 8Q-Agave quantum processor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multireference linearized coupled cluster theory for strongly correlated systems using matrix product states

We propose a multireference linearized coupled cluster theory using matrix product states (MPSs-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles, for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual Density Matrix Renormalization Group (DMRG) algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here, MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including two- and three-dimensional solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How quantum phases on cylinders approach the two-dimensional limit

We consider the properties of $T$ = $0$ quantum phases of matter - especially superconducting and analogous spin-liquid phases - on infinite cylinders of width $L$ ⊥ and analyze the ways in which the $L$ ⊥ → ∞ (two-dimensional, 2d) limit is approached. This problem is interesting in its own right but is particularly important in the context of extrapolating accessible density matrix renormalization group (DMRG) results on model strongly interacting problems to the desired 2d limit. In conclusion, various methods for drawing firm conclusions about the quantum phases in 2d from relatively small $L$ ⊥ results are illustrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange

We report on the use of the adaptively compressed exchange (ACE) operator to accelerate many-body perturbation theory (MBPT) calculations, including G 0 W 0 and the Bethe–Salpeter equation (BSE), for hybrid density functional theory starting points. We show that by approximating the exact exchange operator with the low-rank ACE operator, substantial computational savings can be achieved with systematically controllable errors in the quasiparticle energies computed with full-frequency G 0 W 0 and the optical absorption spectra and vertical excitation energies computed by solving the BSE within density matrix perturbation theory. Our implementation makes use of the ACE-accelerated electronic Hamiltonian to carry out both G 0 W 0 and BSE without explicitly computing empty states. We show the robustness of the approach and present the computational gains obtained on both the central processing unit and graphics processing unit nodes. In conclusion, our work will facilitate the exploration and evaluation of fine-tuned hybrid starting points aimed at enhancing the accuracy of MBPT calculations without involving computationally demanding self-consistency in Hedin’s equations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model dynamics for quantum computing

A model master equation suitable for quantum computing dynamics is presented. In an ideal quantum computer (QC), a system of qubits evolves in time unitarily and, by virtue of their entanglement, interfere quantum mechanically to solve otherwise intractable problems. In the real situation, a QC is subject to decoherence and attenuation effects due to interaction with an environment and with possible short-term random disturbances and gate deficiencies. The stability of a QC under such attacks is a key issue for the development of realistic devices. We assume that the influence of the environment can be incorporated by a master equation that includes unitary evolution with gates, supplemented by a Lindblad term. Lindblad operators of various types are explored; namely, steady, pulsed, gate friction, and measurement operators. In the master equation, we use the Lindblad term to describe short time intrusions by random Lindblad pulses. The phenomenological master equation is then extended to include a nonlinear Beretta term that describes the evolution of a closed system with increasing entropy. An external Bath environment is stipulated by a fixed temperature in two different ways. Here we explore the case of a simple one-qubit system in preparation for generalization to multi-qubit, qutrit and hybrid qubit–qutrit systems. This model master equation can be used to test the stability of memory and the efficacy of quantum gates. The properties of such hybrid master equations are explored, with emphasis on the role of thermal equilibrium and entropy constraints. Several significant properties of time-dependent qubit evolution are revealed by this simple study. - Highlights: • General realistic density matrix evolution model shown for one qubit. • Unitary pulses & bias yield quantum gates without harmful phase accumulations. • Gate friction, dissipation included by random Lindblad noise pulses. • Bath & entropy terms yield equilibrium; test algorithms, gates act on Larmor grid. • Measurement modeled by strong, fast Lindblad pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Natural ionization orbitals for interpreting electron detachment processes

A compact orbital representation of ionization processes is described utilizing the difference of calculated one-particle density matrices. Natural orbital analysis involving this difference density matrix simplifies interpretation of electronic detachment processes and allows differentiation between one-electron transitions and shake-up/shake-off transitions, in which one-electron processes are accompanied by excitation of a second electron into the virtual orbital space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lambda-scheme spectroscopy in the cat-state field

We investigate the spectroscopical scheme of a three-level system in a Λ-configuration utilizing resonant optical field being a superposition of two Glauber coherent states. The source of the field is explicitly included into the theoretical model. In the limit of strong non-classical field we propose an ansatz for solution of the master equation and calculate the steady-state density matrix. The work of both classical and nonclassical field in the steady-state regime was calculated. It was found that the presence of non-classical field creates strong atom-field correlations that results in absence of dark resonances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗