Fermionic Reduced Density Low-Rank Matrix Completion, Noise Filtering, and Measurement Reduction in Quantum Simulations
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “DENSITY MATRIX”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A multisite decomposition of the tensor network path integrals
Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH
We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Accelerating the density-functional tight-binding method using graphical processing units
Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Multiscale characterization and representation of variability in ceramic matrix composites
Low density, high strength, and high creep and oxidation resistance properties of ceramic matrix composites (CMCs) make them an ideal choice for use in extreme environments in space and military applications. This paper presents a detailed characterization study of structural and manufacturing flaws in Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs at different length-scales. Energy-dispersive spectroscopy (EDS) is used for the chemical characterization of the material’s elemental constituents. High-resolution multiscale graphs obtained from scanning electron microscope (SEM) and confocal laser scanning microscope (LSM) are used to characterize the distribution and morphology of defects at different length scales. This is followed by the classification and quantification of the common manufacturing defects. An image processing algorithm based on the image segmentation process is developed to quantify the variability of various scale-dependent architectural parameters. Finally, a three-dimensional stochastic representative volume element (SRVE) generation algorithm is developed to provide precise representations of material textures at multiple length scales. The developed algorithm accurately accounts for material features and flaws based on a range of multiscale structural and defects characterization results.
36 MATERIALS SCIENCE↗
237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes
A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Oxide dispersoid coherency of a ferritic-martensitic 12Cr oxide-dispersion-strengthened alloy under self-ion irradiation
Ferritic-martensitic oxide-dispersion-strengthened (ODS) alloy has shown excellent mechanical property and high radiation tolerance. However, the stability of dispersoids during displacive irradiation in the individual ferritic and martensitic phases is unclear. In this work, the correlation among dispersoid coherency, size, density, and matrix phase are studied in dual-phase 12Cr ODS after 100 peak displacements-per-atom (dpa) irradiation at 475°C, using 3.5 MeV Fe 2+ self-ions. The size and density changes of coherent and incoherent dispersoids were analyzed as a function of irradiation depth. The average dispersoid size decreased after irradiation in both phases, and the large incoherent dispersoids in the tempered martensite phase underwent a more dramatic change than those in the ferrite phase. The dispersoid density significantly increased in the ferrite phase within the irradiated region, mostly resulting in an increase in coherent dispersoid density. On the other hand, only a small change in density was observed in the tempered martensite phase for both coherent and incoherent dispersoids. This study shows that dispersoid evolution is dramatically different in ferrite and tempered martensite phases.
36 MATERIALS SCIENCE↗
Elucidating the impact of fiber source on polypropylene/hemp composite performance for the automotive industry
Given their high strength-to-weight ratio, there is an ever-increasing volume of plastics being used in the automotive industry as plastics aid in the charge to lightweight vehicles for improved fuel efficiency. However, these plastics are often landfilled at their end-of-life, which has given rise to the demand for sustainable materials and waste management alternatives compared to purely synthetic systems. Natural fiber composites have been explored as a viable material option to reduce the environmental impact of plastic use in automobiles while simultaneously ensuring the part performance is not sacrificed. Herein, we explored the use of hemp/polypropylene (PP) composites in which US-sourced hemp is compared to internationally sourced and industrially available hemp. There appears to be a minimal impact on the composite properties regardless of fiber sourcing, and the addition of a natural filler to the PP matrix results in up to a 367% increase in Young's modulus, 126% increase in heat deflection temperature, and comparable water uptake performance. It should be noted that the natural filler addition does increase density by up to 13% due to the higher density of natural fibers compared to the low-density PP matrix. A modified rule of mixtures calculation revealed that the composite materials produced in this study demonstrated good agreement with analytical modeling. Finally, a screening analysis was performed exploring the transportation of hemp fibers, and the results build a strong case for regionalized manufacturing of automotive parts.
Hubbard, Amber M. [Oak Ridge National Laboratory (↗
High-performance strategies for the recent MRSF-TDDFT in GAMESS
Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy
We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems
Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.
Calculus of variations↗
Quasi-Fermi liquid behavior in a one-dimensional system of interacting spinless fermions
We present numerical evidence for a paradigm in one-dimensional interacting fermion systems, whose phenomenology has traits of both Luttinger liquids and Fermi liquids. This state, dubbed a quasi-Fermi liquid, possesses a discontinuity in its fermion occupation number at the Fermi momentum. The excitation spectrum presents particlelike quasiparticles and absence of holelike quasiparticles, giving rise instead to edge singularities. Such a state is realized in a one-dimensional spinless fermion lattice Hamiltonian by fine-tuning the interactions to a regime where they become irrelevant in the renormalization group sense. We show, using uniform infinite matrix products states and finite-entanglement scaling analysis, that the system ground state is characterized by a Luttinger parameter K = 1 and a discontinuous jump in the fermion occupation number. We support the characterization with calculations of the spectral function that show a particle-hole asymmetry reflected in the existence of well-defined Landau quasiparticles above the Fermi level and edge singularities without the associated quasiparticles below. Furthermore, these results indicate that the quasi-Fermi liquid paradigm can be realized beyond the low-energy perturbative realm.
1-dimensional systems↗
Efficient simulation of low-temperature physics in one-dimensional gapless systems
Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O(cN 0 logβ) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Hubbard-corrected oxide formation enthalpies without adjustable parameters
A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.
transition metal oxides↗
Mana and thermalization: Probing the feasibility of near-Clifford Hamiltonian simulation
Quantum hydrodynamics is the emergent classical dynamics governing transport of conserved quantities in generic strongly interacting quantum systems. Here, recent matrix product operator methods [1,2] have made simulations of quantum hydrodynamics in 1+1D tractable, but they do not naturally generalize to 2+1D or higher, and they offer limited guidance as to the difficulty of simulations on quantum computers. Near-Clifford simulation algorithms are not limited to one dimension, and future error-corrected quantum computers will likely be bottlenecked by non-Clifford operations. We therefore investigate the non-Clifford resource requirements for simulation of quantum hydrodynamics using mana, a resource theory of non-Clifford operations. For infinite-temperature starting states, we find that the mana of subsystems quickly approaches zero, while for starting states with energy above some threshold the mana approaches a nonzero value. Surprisingly, in each case the finite-time mana is governed by the subsystem entropy, not the thermal state mana; we argue that this is because mana is a sensitive diagnostic of finite-time deviations from canonical typicality.
1-dimensional spin chains↗
Spectral function of the $J_1$ – $J_2$ Heisenberg model on the triangular lattice
Spectral probes, such as neutron scattering, are crucial for characterizing excitations in quantum many-body systems and the properties of quantum materials. Among the most elusive phases of matter are quantum spin liquids, which have no long-range order even at zero temperature and host exotic fractionalized excitations with nontrivial statistics. These phases can occur in frustrated quantum magnets, such as the paradigmatic Heisenberg model with nearest- and next-nearest-neighbor exchange interactions on the triangular lattice, the so-called $J_1$ – $J_2$ model. In this work, we compute the spectral function using large-scale matrix product state simulations across the three different phases of this model's phase diagram, including a quantum spin-liquid phase at intermediate $J_2$/$J_1$. Despite a plethora of theoretical and experimental studies, the exact nature of this phase is still contested, with the dominant candidates being a gapped $\mathbb{Z}_2$, a gapless U(1) Dirac, and a spinon Fermi surface quantum spin-liquid state. We find a V-shaped spectrum near the center of the Brillouin zone ($Γ$ point), a key signature of a spinon Fermi surface, observed in prior neutron scattering experiments. However, we find a small gap near the $Γ$ point, ruling out such a phase. Furthermore, we find localized gapless excitations at the corner of the Brillouin zone boundary (K point) and the middle of the edge of the Brillouin zone boundary (M point), ruling out the gapped $\mathbb{Z}_2$ spin-liquid phase. Here, our results imply that the intermediate spin-liquid phase is a gapless $U(1)$ Dirac spin liquid, and provide clear signatures to detect this phase in future neutron scattering experiments.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Distributed out-of-memory NMF on CPU/GPU architectures
We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .
97 MATHEMATICS AND COMPUTING↗
Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction
Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗