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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Many-fermion simulation from the contracted quantum eigensolver without fermionic encoding of the wave function

Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding—a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schrödinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schrödinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen O 2 , and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Modeling the sequential dissociative double ionization of O 2 by ultrashort intense infrared laser pulses

A density matrix approach for sequential double ionization (DM-SDI) of molecules has been developed recently and was applied to the N 2 molecule. In this article, we extended the DM-SDI model to O 2 , which is a more complicated system to model than N 2 , due to its electronic structures and spin-orbit and laser couplings in the manifold of doubly charged states. We obtained a good agreement on the kinetic energy release spectrum of O + + O + from previous experiments. Thanks to the low computational cost of the model, we explored the mechanism behind the ionization and dissociation dynamics as well as the effects of lasers on the spectrum. Furthermore, this work will pave the way to model sequential dissociative double ionization of larger molecules and to probe molecular dynamics by measuring kinetic energy release spectra from this process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Coherence from multiorbital tunneling ionization of molecules

We present a simple and general coherence model for multiorbital tunnel ionization of molecules, which we incorporate into our previously developed density-matrix approach for sequential double ionization. Here, the influence of this coherence is investigated through simulations of single ionization and sequential double ionization of N 2 and O 2 using few-cycle near-infrared laser pulses. In the case of single ionization, our results reveal the crucial role played by this coherence in generating population inversion in N 2 + , suggesting a potential mechanism for air lasing. Regarding sequential double ionization, we observe only minor changes in the kinetic energy release spectra when the coherence is included, while noticeable differences in the angle-dependent dication yield for both N 2 and O 2 are found. Based on these findings, we recommend the inclusion of multiorbital tunnel ionization coherence in models for single ionization of general molecules, while suggesting that it can be safely neglected in the case of sequential double ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Theory of strong-field sequential double ionization of polyatomic molecules

Understanding of strong-field sequential double ionization (SDI) of molecules by a highly intense infrared laser pulse could be the key to observing electron motion in molecules in the attosecond to femtosecond timescale. Based on a novel density matrix approach for SDI (DM-SDI), a recent theoretical study has shown that SDI can be used as a probe to monitor the changes in vibronic coherence in homonuclear diatomic molecules. In this article, we extend the DM-SDI model to general molecules that could possess permanent dipole moments and arbitrary symmetry. We apply the model to SDI of a water molecule and identify the formation pathway of individual dication states. We further deduce that kinetic energy release spectra from two- and three-body fragments could carry the signature of vibronic coherence between the lowest two states of the water cation. Furthermore, our results suggest that observables from the SDI probe can be interpreted intuitively with only the knowledge of electronic structures of populated ionic states, making the SDI probe to be a highly desirable probing scheme for vibronic coherence in generic molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Probing vibronic coherence in charge migration in molecules using strong-field sequential double ionization

We propose a scheme for probing vibronic coherence in charge migration in molecules utilizing strong-field sequential double ionization. To demonstrate the feasibility of this approach, we perform full simulations of a pump-probe scheme employing few-cycle intense infrared pulses for N 2 and O 2 . We predict that the vibronic coherence between the pumped states will be directly imprinted in experimental observables such as kinetic energy release spectra and branching ratios of the dissociative dications. Our simulations are based on the recently developed density matrix approach for sequential double ionization model, which is capable of efficiently accounting for molecular orientations and enabling direct comparison with experimental results. Lastly, our findings strongly encourage the use of this probing scheme in future charge migration experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Nonlocal polarization interferometer for entanglement detection

We report a nonlocal interferometer capable of detecting entanglement and identifying Bell states statistically. This is possible due to the interferometer's unique correlation dependence on the antidiagonal elements of the density matrix, which have distinct bounds for separable states and unique values for the four Bell states. The interferometer consists of two spatially separated balanced Mach-Zehnder or Sagnac interferometers that share a polarization-entangled source. Correlations between these interferometers exhibit nonlocal interference, while single-photon interference is suppressed. This interferometer also allows for a unique version of the Clauser-Horne-Shimony-Holt Bell test where the local reality is the photon polarization. In conclusion, we present the relevant theory and experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Dirac open-quantum-system dynamics: Formulations and simulations

In this work, we present an open-system interaction formalism for the Dirac equation. Overcoming a complexity bottleneck of alternative formulations, our framework enables efficient numerical simulations (utilizing a typical desktop) of relativistic dynamics within the von Neumann density matrix and Wigner phase-space descriptions. Employing these instruments, we gain important insights into the effect of quantum dephasing for relativistic systems in many branches of physics. In particular, the conditions for robustness of Majorana spinors against dephasing are established. Using the Klein paradox and tunneling as examples, we show that quantum dephasing does not suppress negative energy particle generation. Hence, the Klein dynamics is also robust to dephasing.

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

Unconventional fermionic pairing states in a monochromatically tilted optical lattice

We study the one-dimensional attractive fermionic Hubbard model under the influence of periodic driving with the time-dependent density matrix renormalization group method. We show that the system can be driven into an unconventional pairing state characterized by a condensate made of Cooper pairs with a finite center-of-mass momentum similar to a Fulde-Ferrell state. We obtain results both in the laboratory and the rotating reference frames demonstrating that the momentum of the condensate can be finely tuned by changing the ratio between the amplitude and the frequency of the driving. In particular, by quenching this ratio to the value corresponding to suppression of the tunneling and the Coulomb interaction strength to zero, we are able to “freeze” the condensate. We finally study the effects of different initial conditions and compare our numerical results to those obtained from a time-independent Floquet theory in the large frequency regime. Lastly, our work offers the possibility of engineering and controlling unconventional pairing states in fermionic condensates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Path-integral Monte Carlo study of particles obeying quantum mechanics and classical statistics

Ultracold atomic systems have been of great research interest in the past, with more recent attention being paid to systems of mixed species. In this work, we carry out nonperturbative path-integral Monte Carlo (PIMC) simulations of $N$ distinguishable particles at finite temperature, which can be thought of as an ultracold atomic system containing $N$ distinct species. We use the PIMC approach to calculate thermodynamic properties of particles interacting via hard-sphere and hard-cavity potentials. The first problem we study is a two-particle system interacting via a hard-sphere and hard-cavity interaction in order to test the effectiveness of two approximations for the thermal density matrix corresponding to these potentials. We then apply the PIMC method to a system of many hard-sphere particles under periodic boundary conditions at varying temperature in order to calculate the energy per particle, pressure, and specific heat of the system. We examine how finite-size effects impact the results of PIMC simulations of hard-sphere particles and when the thermodynamic limit has been reached. Finally, our results provide microscopic benchmarks for a system containing distinguishable particles, which can be thought of as a limiting case for ultracold atomic systems of mixed species.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Highly dispersive magnons with spin-gap-like features in the frustrated ferromagnetic S = 1 2 chain compound Ca 2 Y 2 Cu 5 O 10 detected by inelastic neutron scattering

In this paper, we report inelastic neutron scattering experiments in Ca 2 Y 2 Cu 5 O 10 and map out the full one-magnon dispersion which extends up to a record value of 53 meV for frustrated ferromagnetic (FM) edge-sharing CuO 2 chain (FFESC) cuprates. A homogeneous spin-1/2 chain model with a FM nearest-neighbor (NN), an antiferromagnetic (AFM) next-nearest-neighbor (NNN) inchain, and two diagonal AFM interchain couplings (ICs) analyzed within linear spin-wave theory (LSWT) reproduces well the observed strong dispersion along the chains and a weak one perpendicularly. The ratio α=| J a 2 / J a 1 | of the FM NN and the AFM NNN couplings is found as ~0.23, close to the critical point α c =1/4 which separates ferromagnetically and antiferromagnetically correlated spiral magnetic ground states in single chains, whereas α c >0.25 for coupled chains is considerably upshifted even for relatively weak IC. Although the measured dispersion can be described by homogeneous LSWT, the scattering intensity appears to be considerably reduced at ~11.5 and ~28 meV. The gaplike feature at 11.5 meV is attributed to magnon-phonon coupling whereas based on density matrix renormalization group simulations of the dynamical structure factor the gap at 28 meV is considered to stem partly from quantum effects due to the AFM IC. Another contribution to that gap is ascribed to the intrinsic superstructure from the distorting incommensurate pattern of CaY cationic chains adjacent to the CuO 2 ones. It gives rise to nonequivalent CuO 4 units and Cu-O-Cu bond angles Φ and a resulting distribution of all exchange integrals. The J 's fitted by homogeneous LSWT are regarded as average values. The record value of the FM NN integral J 1 =24 meV among FFESC cuprates can be explained by a nonuniversal Φ(≠90°) and Cu-O bond length dependent anisotropic mean direct FM Cu-O exchange $\overline{K}_{pd}$~120 meV, similar to a value of 105 meV for Li 2 CuO 2 , in accord with larger values for La 2 CuO 4 and CuGeO 3 (~110 meV) reported by Braden et al . [Phys. Rev. B 54, 1105 (1996)] phenomenologically. Lastly, enhanced $K_{pd}$ values are also needed to compensate a significant AFM J dd ≥ 6 meV from the dd channel, generic for FFESC cuprates but ignored so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Dynamical mean field theory simulations with the adaptive sampling configuration interaction method

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero-temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here, we present the ASCI-DMFT method and example calculations on the one- and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground-state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Finally, our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Phase transitions in the Rényi entropies of a 2+1D large- N interacting vector model

Here we consider the reduced density matrix for a disk in the ground state of the interacting fixed points of the large- N O(N) vector model in 2 + 1 dimensions. Using the map to the free energy on $\mathbb H_2 × S^1$, we show that there is an instability in the Rényi entropies at n = 1 in the N → ∞ limit, indicating a phase transition at or near this Rényi parameter. We close with a discussion of the finite- N case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Nematic spin liquid phase in a frustrated spin-1 system on the square lattice

Frustration in quantum spin systems promote a variety of novel quantum phases. An important example is the frustrated spin-1 model on the square lattice with the nearest-neighbor bilinear ($J_1$) and biquadratic ($K_1$) interactions. Here we provide strong evidence for a nematic spin liquid phase in a range of $K_1/J_1$ near the SU(3)-symmetric point ($J_1 = K_1$), based on the linear flavor-wave theory and extensive density matrix renormalization group calculation. This phase displays no spin dipolar or quadrupolar order, preserves translational symmetry but spontaneously breaks $C_4$ lattice rotational symmetry, and possesses fluctuations peaked at the wave vector (π, 2π/3) . The spin excitation gap drops rapidly with system size and appears to be gapless, and the nematic order is attributed to the dominant (π, 2π/3) fluctuations. Our results provide a novel mechanism for electronic nematic order and, more generally, open up a new avenue to explore frustration-induced exotic ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Chiral spin liquid with spinon Fermi surfaces in the spin-$\frac{1}{2}$ triangular Heisenberg model

Here, we study the interplay of competing interactions in spin-$\frac{1}{2}$ triangular Heisenberg model through tuning the first- ($J_1$), second- ($J_2$), and third-neighbor ($J_ 3$) couplings. Based on a large-scale density-matrix renormalization group calculation, we identify a quantum phase diagram of the system and discover a gapless chiral spin-liquid (CSL) phase in the intermediate $J_2$ and $J_ 3$ regime. This CSL state spontaneously breaks time-reversal symmetry with finite scalar chiral order, and it has gapless excitations implied by a vanishing spin triplet gap and a finite central charge on the cylinder. Moreover, the central charge grows rapidly with the cylinder circumference, indicating emergent spinon Fermi surfaces. To understand the numerical results we propose a parton mean-field spin-liquid state, the U (1) staggered flux state, which breaks time-reversal symmetry with chiral edge modes by adding a Chern insulator mass to Dirac spinons in the U (1) Dirac spin liquid. This state also breaks lattice rotational symmetries and possesses two spinon Fermi surfaces driven by nonzero $J_2$ and $J_3$, which naturally explains the numerical results. This realizes an example of a gapless CSL state with coexisting spinon Fermi surfaces and chiral edge states, demonstrating the rich family of interesting quantum phases emergent from competing interactions in triangular-lattice magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Continuous phase transition between bosonic integer quantum Hall liquid and a trivial insulator: Evidence for deconfined quantum criticality

The deconfined quantum critical point, a prototype Landau-forbidden transition, could exist, in principle, in the phase transitions involving a symmetry-protected topological phase; however, examples of such kinds of transitions in physical systems are rare beyond one-dimensional systems. Here, using a density-matrix renormalization-group calculation, we unveil a bosonic integer quantum Hall phase in a two-dimensional correlated honeycomb lattice, by full identification of its internal structure from the topological K matrix. Furthermore, we demonstrate that imbalanced periodic chemical potentials can destroy the bosonic integer quantum Hall state and drive it into a featureless trivial (Mott) insulator, where all physical observables evolve smoothly across the critical point. At the critical point, the entanglement entropy reveals a characteristic scaling behavior, which is consistent with the critical field theory as an emergent QED3 with two flavors of Dirac fermions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-order harmonic generation in solid C 60

High-order harmonic generation (HHG) has unleashed the power of strong laser physics in solids. Here we investigate HHG from a large system, solid C 60 , with 240 valence electrons engaging harmonic generation at each crystal momentum. We employ the density functional theory and the time-dependent Liouville equation of the density matrix to compute HHG signals. We find that under a moderately strong laser pulse, HHG signals reach 15th order, consistent with the experimental results from C 60 plasma. The helicity dependence in solid C 60 is weak, due to the high symmetry. In contrast to the general belief, HHG is unsuitable for band structure mapping in C 60 . However, we find a window of opportunity using a long wavelength, where harmonics are generated through multiple-photon excitation. In particular, the fifth-order harmonic energies closely follow the transition energy dispersion between the valence and conduction bands. Finally, this finding is expected to motivate future experimental investigations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetic phase transitions in quantum spin-orbital liquids

Here, we investigate the spin and orbital correlations of a superexchange model with spin $\textit{S}$ = 1 and orbital $\textit{L}$ = 1 relevant for $5d^4$ transition-metal Mott insulators, using exact diagonalization and density matrix renormalization group (DMRG). For spin-orbit coupling λ = 0, the orbitals are in an entangled state that is decoupled from the spins. We find two phases with increasing λ: (I) the $\textit{S}$2 phase with two peaks in the structure factor for λ ≤ λ$_{c1}$ ≈ 0.34$\textit{J}$ where $\textit{J}$ is the ferromagnetic exchange; and, (II) the $\textit{S}$1 phase for λ$_{c1}$ < λ ≤ λ$_{c2}$ ≈ 1.2$\textit{J}$ with emergent antiferromagnetic correlations. Both $\textit{S}$1 and $\textit{S}$2 phases are shown to exhibit power-law correlations, indicative of a gapless spectrum. Upon increasing λ > λ$_{c2}$ leads to a product state of local spin-orbital singlets that exhibit exponential decay of correlations, indicative of a gapped phase. We obtain insights into the phases from the well-known Uimin-Lai-Sutherland model in an external field that provides an approximate description of our model within mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Two-magnon bound states in the Kitaev model in a [111] field

It is now well established that the Kitaev honeycomb model in a magnetic field along the [111] direction harbors an intermediate gapless quantum spin liquid (QSL) phase sandwiched between a gapped non-Abelian QSL at low fields $H < H_{c1}$ and a partially polarized phase at high fields $H > H_{c2}$. Here, we analyze the low-field and high-field phases and phase transitions in terms of single- and two-magnon excitations using exact diagonalization and density matrix renormalization group methods. In this work, we find that the energy to create a bound state of two-magnon $Δ_p$ becomes lower than the energy to create a single spin flip $Δ_s$ near $H_{c2}$. In the entire Kitaev spin liquid, $Δ_p < Δ_s$ and both gaps vanish at $H_{c2}$. We make testable predictions for magnon pairing that could be observable in Raman scattering measurements on Kitaev QSL candidate materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗