Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

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

Multidimensional coherent spectroscopy of light-driven states and their collective modes in multiband superconductors

We present a comprehensive theory of light-controlled multiband superconductivity, and apply it to predict distinctive signatures of light-driven superconducting (SC) states in terahertz multidimensional coherent spectroscopy (THz-MDCS) experiments. We first derive gauge-invariant Maxwell-Bloch equations for multiband BCS superconductors with spatial fluctuations. We consider driving electromagnetic fields determined self-consistently by Maxwell's equations. By calculating the THz-MDCS spectra measured experimentally in the clean SC limit, we identify unique signatures of finite-momentum Cooper-pairing states that live longer than the laser pulse. They are controlled by a pair of THz laser pulses with well-defined relative phase (pulse pair). The pseudospin oscillators that describe the properties of these SC states are parametrically driven by both finite-momentum Cooper pairing and by time oscillations of the order-parameter relative phase. We show that such strong parametric driving leads to drastic changes in the THz-MDCS spectral shape from the predictions of third-order nonlinear susceptibility calculations. These spectral changes strongly depend on the interband-to-intraband interaction ratio and on the collective modes of the light-driven state. For negligible interband interaction, the spectra show a transition with increasing field, from traditional pump-probe, four-wave-mixing, and third-harmonic generation peaks determined by the laser frequency to sidebands determined by the excitations of the driven system. These sidebands emerge from difference-frequency Raman processes in the nonequilibrium SC state. For interband couplings weaker than the intraband pairing, we show that the Leggett phase collective mode leads to harmonic sidebands around the traditional pump-probe peaks. Additional Higgs collective mode peaks result from light-induced inversion-symmetry breaking in a thin-film geometry. For strong interband coupling, we find a transition from a nonequilibrium finite Cooper-pair momentum state characterized by hybrid-Higgs amplitude mode peaks in THz-MDCS spectra to a driven state identified experimentally by the emergence of Floquet-type sidebands at bi-Higgs frequencies. Those dominant bi-Higgs-frequency satellites are manifestations of a new order parameter relative phase collective mode that characterizes the nonequilibrium SC state. The predicted interaction- and field-dependent transitions in the spectral profile allow us to propose THz-MDCS experiments for quantum tomography of light-driven superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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

Recovery of Damaged Information and the Out-of-Time-Ordered Correlators

The evolution with a complex Hamiltonian generally leads to information scrambling. A time-reversed dynamics unwinds this scrambling and thus leads to the original information recovery. In this work, we show that if the scrambled information is, in addition, partially damaged by a local measurement, then such a damage can still be treated by application of the time-reversed protocol. This information recovery is described by the long-time saturation value of a certain out-of-time-ordered correlator of local variables. We also propose a simple test that distinguishes between quantum and reversible classical chaotic information scrambling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effective one-band models for the one-dimensional cuprate Ba 2-x Sr x CuO 3+δ

In this work, we consider a multiband Hubbard model H m for Cu and O orbitals in Ba 2-x Sr x CuO 3+δ similar to the three-band model for two-dimensional cuprates. The hopping parameters are obtained from maximally localized Wannier functions derived from ab initio calculations. Using the cell perturbation method, we derive both a generalized t–J model H tJ and a one-band Hubbard model H H to describe the low-energy physics of the system. H tJ has the advantage of having a smaller relevant Hilbert space, facilitating numerical calculations, while additional terms should be included in H H to accurately describe the multiband physics of H m . Using H tJ and the density matrix renormalization group method, we calculate the wave-vector-resolved photoemission and discuss the relevant features in comparison with recent experiments. In agreement with previous calculations, we find that the addition of an attractive nearest-neighbor interaction of the order of the nearest-neighbor hopping shifts the weight from the 3k F to the holon-folding branch. Kinetic effects also contribute to this process.

1-dimensional systems↗
DOE OSTI2020

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

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↗
DOE OSTI2025

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2021

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↗
DOE OSTI2025

A multidimensional approach to quantum state tomography of photoelectron wavepackets

There is a growing interest in reconstructing the density matrix of photoelectron wavepackets, in particular in complex systems where decoherence can be introduced either by a partial measurement of the system or through coupling with a stochastic environment. To this end, several methods to reconstruct the density matrix, quantum state tomography protocols, have been developed and tested on photoelectrons ejected from noble gases following absorption of extreme ultraviolet (XUV) photons from attosecond pulses. It remains a challenge to obtain model-free, single scan protocols that can reconstruct the density matrix with high fidelities. Current methods require extensive measurements or involve complex fitting of the signal. Efficient single-scan reconstructions would be of great help to increase the number of systems that can be studied. We propose a new and more efficient protocol that is able to reconstruct the continuous variable density matrix of a photoelectron in a single time delay scan. It is based on measuring the coherences of a photoelectron created by absorption of an XUV pulse using a broadband infrared (IR) probe that is scanned in time and a narrowband IR reference that is temporally fixed to the XUV pulse. We illustrate its performance for a Fano resonance in He as well as mixed states in Ar arising from spin-orbit splitting. We show that the protocol results in excellent fidelities and near-perfect estimation of the purity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Root-N Krylov-space correction vectors for spectral functions with the density matrix renormalization group

In this work, we propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov-space decomposition to compute the correction vectors. Our approach entails the calculation of the root-N (N=2 is the standard square root) of the Hamiltonian propagator using Krylov-space decomposition and repeating this procedure N times to obtain the actual correction vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, whereas achieving better computational performance at large N. We apply our method to spin and charge spectral functions of t-J and Hubbard models in the challenging two-leg ladder geometry and provide evidence that the root-N approach reaches a much improved spectral resolution compared to the conventional correction vector.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

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↗
DOE OSTI2021

Nonequilibrium optical response of a one-dimensional Mott insulator

We define, compute, and analyze the nonequilibrium differential optical conductivity of the one-dimensional extended Hubbard model at half-filling after applying a pump pulse, using the time-dependent density matrix renormalization group method. The melting of the Mott insulator is accompanied by a suppression of the local magnetic moment and ensuing photogeneration of doublon-holon pairs. The differential optical conductivity reveals (i) mid-gap states related to parity-forbidden optical states and (ii) strong renormalization and hybridization of the excitonic resonance and the absorption band, yielding a Fano resonance. We offer evidence and interpret such a resonance as a signature of nonequilibrium optical excitations resembling excitonic strings, (bi)excitons, and unbound doublon-holon pairs, depending on the magnitude of the intersite Coulomb repulsion. We discuss our results in the context of pump and probe spectroscopy experiments on organic Mott insulators.

1-dimensional systems↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy-particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). In this work, we apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

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

Intertwined charge, spin, and pairing orders in doped iron ladders

Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe2S3. The model is constructed by using ab initio hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the (π,0) magnetic ordering at half filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard U coupling, other magnetic patterns, such as (0,π), are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intraorbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge density wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such a scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗