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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantum simulation of molecules without fermionic encoding of the wave function

Abstract Molecular simulations generally require fermionic encoding in which fermion statistics are encoded into the qubit representation of the wave function. Recent calculations suggest that fermionic encoding of the wave function can be bypassed, leading to more efficient quantum computations. Here we show that the two-electron reduced density matrix (2-RDM) can be expressed as a unique functional of the unencoded N -qubit-particle wave function without approximation, and hence, the energy can be expressed as a functional of the 2-RDM without fermionic encoding of the wave function. In contrast to current hardware-efficient methods, the derived functional has a unique, one-to-one (and onto) mapping between the qubit-particle wave functions and 2-RDMs, which avoids the over-parametrization that can lead to optimization difficulties such as barren plateaus. An application to computing the ground-state energy and 2-RDM of H 4 is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Methodology for Producing a Uniform Distribution of UO2 in a Tungsten Matrix

Current work at NASA's Marshall Space Flight Center (MSFC) is focused on the development CERMET fuel materials for Nuclear Thermal Propulsion (NTP). The CERMETs consist of uranium dioxide (UO2) fuel particles embedded in a tungsten (W) metal matrix. Initial testing of W-UO2 samples fabricated from fine angular powders performed reasonably well, but suffered from significant fuel loss during repeated thermal cycling due to agglomeration of the UO2 (1). The blended powder mixtures resulted in a non-uniform dispersion of the UO2 particles in the tungsten matrix, which allows rapid vaporization of the interconnected UO2 from the sample edges into the bulk material. Also, the angular powders create areas of stress concentrations due to thermal expansion mismatch, which eventually cracks the tungsten matrix. Evenly coating spherical UO2 particles with chemical vapor deposited (CVD) tungsten prior to consolidation was previously demonstrated to provide improved performance. However, the CVD processing technology is expensive and not currently available. In order to reduce cost and enhance performance, a powder coating process has been developed at MSFC to produce a uniform distribution of the spherical UO2 particles in a tungsten matrix. The method involves utilization of a polyethylene binder during mixing which leads to fine tungsten powders clinging to the larger UO2 spherical particles. This process was developed using HfO2 as a surrogate for UO2. Enough powder was mixed to make 8 discs (2cm diameter x 8mm thickness) using spark plasma sintering. A uniaxial pressure of 50 MPa was used at four different temperatures (2 samples at each temperature). The first two samples were heated to 1400C and 1500C respectively for 5 minutes. Densities for these samples were less than 85% of theoretical, so the time at temperature was increased to 20 minutes for the remaining samples. The highest densities were achieved for the two samples sintered at 1700C (approx. 92% of theoretical). Scanning electron microscopy (SEM) of the mixed powders and the sintered samples along with energy dispersive x-ray analysis was obtained. The SEM of the powders clearly show the fine W powder adhered to the larger HfO2 particles and a uniform distribution of HfO2 particles in a tungsten matrix upon densification. Vicker's Microhardness testing was also performed on all samples using 0.5, 1.0 and 2.0 kg loads. Five indents were made at each load level. All indents were placed in the tungsten matrix to assist as a proxy in measuring densification. The highest hardness value was obtained for the 1700C specimens. The hardness average for these samples was 312.14 MPa. This powder processing method has been applied to W/UO2 powders with the SEM of the powders appearing similar to the W/HfO2 powder images.

Tucker, Dennis S.↗
NASA NTRS1974

Elementary solutions of coupled model equations in the kinetic theory of gases

The method of elementary solutions is employed to solve two coupled integrodifferential equations sufficient for determining temperature-density effects in a linearized BGK model in the kinetic theory of gases. Full-range completeness and orthogonality theorems are proved for the developed normal modes and the infinite-medium Green's function is constructed as an illustration of the full-range formalism. The appropriate homogeneous matrix Riemann problem is discussed, and half-range completeness and orthogonality theorems are proved for a certain subset of the normal modes. The required existence and uniqueness theorems relevant to the H matrix, basic to the half-range analysis, are proved, and an accurate and efficient computational method is discussed. The half-space temperature-slip problem is solved analytically, and a highly accurate value of the temperature-slip coefficient is reported.

Kriese, J. T.↗
DOE OSTI2019

Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. Here, the accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics

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

Continuous microstructure evolution of ball rolled TWIP Ti15Mo alloy and its effect on phase precipitation behavior

Highlights: • Gradient microstructure was induced by ball rolling method in a Ti15Mo alloy. • Twin density increased and then decreased as increasing distance from surface. • Ultrafine (α + β) duplex structure was formed in the heavily deformed region. • α-phase in plate was precipitated in twined region in slightly deformed region. • Strengthening by formation of sub/nano-grains was notable than that by twins. - Abstract: A β-type Ti15Mo alloy with {332}<113> twinning-induced plasticity (TWIP) effect was utilized to investigate the continuous microstructure evolution of β-matrix by ball rolling method. Precipitation behavior of α-phase and strengthening mechanism were examined in regions with different strain. Amounts of intersected twins were present in the slightly deformed region of β-matrix, and their area fraction decreased with an increase of distance from rolled surface. The β-matrix close to surface was heavily deformed to result in the formation of abundant subgrains and nanograins, which was caused by cooperative occurrence of twinning refinement, intragranular dislocation activities and subgrain rotation. Mechanical twins provided the preferential nucleation sites for plate-like α-phase with strong variant selection. Equiaxed α-phase preferentially nucleated at the boundaries of subgrains and nanograins with weak variant selection to form the ultrafine (α + β) duplex structure. Due to the formation of subgrains and nanograins as well as the strain hardening, strengthening became more significant than that by twinning refinement solely. Similarly, the ultrafine (α + β) duplex structure exhibited a better strengthening effect than that by precipitation of plate-like α-phase. The diversities of deformation microstructure combined with different precipitation behaviors of α-phase are desirable for controlling the mechanical properties of β-type TWIP titanium alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

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

Variational mixed quantum/semiclassical simulation of dihalogen guest and rare-gas solid host dynamics

A variational mixed quantum-semiclassical theory for the internal nuclear dynamics of a small molecule and the induced small-amplitude coherent motion of a low-temperature host medium is developed, tested, and used to simulate the temporal evolution of nonstationary states of the internal molecular and surrounding medium degrees of freedom. In this theory, termed the Fixed Vibrational Basis/Gaussian Bath (FVB/GB) method, the system is treated fully quantum mechanically while Gaussian wave packets are used for the bath degrees of freedom. An approximate time-dependent wave function of the entire model is obtained instead of just a reduced system density matrix, so the theory enables the analysis of the entangled system and bath dynamics that ensues following initial displacement of the internal-molecular (system) coordinate from its equilibrium position. The norm- and energy-conserving properties of the propagation of our trial wave function are natural consequences of the Dirac-Frenkel-McLachlan variational principle. The variational approach also stabilizes the time evolution in comparison to the same ansatz propagated under a previously employed locally quadratic approximation to the bath potential and system-bath interaction terms in the bath-parameter equations of motion. Dynamics calculations are carried out for molecular iodine in a 2D krypton lattice that reveal both the time-course of vibrational decoherence and the details of host-atom motion accompanying energy dissipation and dephasing. This work sets the stage for the comprehensive simulation of ultrafast time-resolved optical experiments on small molecules in low-temperature solids.

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

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Generalization of the Kohn-Sham system that can represent arbitrary one-electron density matrices

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of noninteracting particles, is the workhorse of the theory. The particular form of the Kohn-Sham wave function admits only idempotent one-electron density matrices whereas wave functions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept and introducing a suitable dot product as well as a probability density, a noninteracting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans' theorem are available. Lastly, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single-determinant wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Toward quantum networking with frequency-bin qudits

Quantum networking holds tremendous promise in transforming computation and communication. Entangled-photon sources are critical for quantum repeaters and networking, while photonic integrated circuits are vital for miniaturization and scalability. In this talk, we focus on generating and manipulating frequency-bin entangled states within integrated platforms. We encode quantum information as a coherent superposition of multiple optical frequencies; this approach is favorable due to its amenability to high-dimensional entanglement and compatibility with fiber transmission. We successfully generate and measure the density matrix of biphoton frequency combs from integrated silicon nitride microrings, fully reconstructing the state in an 8 × 8 two-qudit Hilbert space, the highest so far for frequency bins. Moreover, we employ Vernier electro-optic phase modulation methods to perform time-resolved measurements of biphoton correlation functions. Currently, we are exploring bidirectional pumping of microrings to generate indistinguishable entangled pairs in both directions, aiming to demonstrate key networking operations such as entanglement swapping and Greenberger–Horne–Zeilinger state generation in the frequency domain.

Myilswamy, Karthik V.↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

On minimizing the probability of misclassification for linear feature selection

The use of techniques for feature selection permits treatment of classification problems in spaces of reduced dimensions. A method is considered of linear feature selection for n-dimensional observation vectors which belong to one of two populations, where each population is described by a known multivariate normal density function. More specifically, the problem of finding a 1xn transformation matrix B for which the probability of misclassification with respect to the one-dimensional transformed density functions was minimized was considered. Theoretical results are presented which give rise to a numerically tractable expression for the variation in the probability of misclassification with respect to B. Using this expression a computational procedure is discussed for obtaining a B which minimizes the probability of misclassification. Preliminary numerical results are discussed.

Guseman, L. F., Jr.↗
NASA NTRS1991

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗
NASA NTRS2006

High-Performance Scanning Acousto-Ultrasonic System

A high-performance scanning acousto-ultrasonic system, now undergoing development, is designed to afford enhanced capabilities for imaging microstructural features, including flaws, inside plate specimens of materials. The system is expected to be especially helpful in analyzing defects that contribute to failures in polymer- and ceramic-matrix composite materials, which are difficult to characterize by conventional scanning ultrasonic techniques and other conventional nondestructive testing techniques. Selected aspects of the acousto-ultrasonic method have been described in several NASA Tech Briefs articles in recent years. Summarizing briefly: The acousto-ultrasonic method involves the use of an apparatus like the one depicted in the figure (or an apparatus of similar functionality). Pulses are excited at one location on a surface of a plate specimen by use of a broadband transmitting ultrasonic transducer. The stress waves associated with these pulses propagate along the specimen to a receiving transducer at a different location on the same surface. Along the way, the stress waves interact with the microstructure and flaws present between the transducers. The received signal is analyzed to evaluate the microstructure and flaws. The specific variant of the acousto-ultrasonic method implemented in the present developmental system goes beyond the basic principle described above to include the following major additional features: Computer-controlled motorized translation stages are used to automatically position the transducers at specified locations. Scanning is performed in the sense that the measurement, data-acquisition, and data-analysis processes are repeated at different specified transducer locations in an array that spans the specimen surface (or a specified portion of the surface). A pneumatic actuator with a load cell is used to apply a controlled contact force. In analyzing the measurement data for each pair of transducer locations in the scan, the total (multimode) acousto-ultrasonic response of the specimen is utilized. The analysis is performed by custom software that extracts parameters of signals in the time and frequency domains. The computer hardware and software provide both real-time and postscan processing and display options. For example, oscilloscope displays of waveforms and power spectral densities are available in real time. Images can be computed while scanning continues. Signals can be digitally preprocessed and/or post-processed by filtering, windowing, time-segmenting, and running-waveform-averaging algorithms. In addition, the software affords options for off-line simulation of the waveform-data-acquisition and scanning processes. In tests, the system has been shown to be capable of characterizing microstructural changes and defects in SiC/SiC and C/SiC ceramic-matrix composites. Delaminations, variations in density, microstructural changes attributable to infiltration by silicon, and crack-space indications (defined in the next sentence) have been revealed in images formed from several time- and frequency-domain parameters of scanning acousto-ultrasonic signals. The crack-space indications were image features that were not revealed by other nondestructive testing methods and are so named because they turned out to mark locations where cracking eventually occurred.

Roth, Don↗
DOE OSTI2014

SU-E-I-41: Non-Cartesian MR Image Reconstruction with Integrated Gradient Non-Linearity Correction

Purpose: Nonlinearities in the spatial encoding gradients of MRI systems cause geometric distortion in images. Typically, this is retrospectively corrected via image-domain interpolation (a.k.a., “gradwarp”) albeit with a loss of spatial resolution. For non-Cartesian MRI, the latter problem is exaggerated by noise and undersampling artifact. In this study, we describe a novel correction strategy that accounts for gradient nonlinearities during — rather than after — non-Cartesian MRI reconstruction, and demonstrate that this approach mitigates the resolution loss that can occur with standard methods. Methods: To test the proposed method, the American College of Radiology (ACR) quality control phantom was scanned on at 1.5 T (General Electric, v16.0, “zoom” gradient) using a 1.6x undersampled 3D non- Cartesian Shells trajectory (GRE, FOV=24 cm3, 120 shells, 16552 shots, 512 readout, matrix=2403). Image reconstruction was first performed via standard k-space density-compensated gridding and retrospectively corrected via cubic spline interpolation. Image reconstruction was then separately performed using a k-space and image-domain densitycompensated type-3 non-uniform fast Fourier transform (NUFFT), which provides a direct mapping between non-Cartesian k-space samples and warped image space voxel locations. Thus, no separate distortion correction procedure is needed for the proposed approach. The gradient distortion field was determined using vendor provided calibration data. Results: Phantom scan results show that both processing approaches successfully correct geometric distortion. However, visual inspection of the ACR phantom spatial resolution inserts shows that the proposed strategy preserves the resolution of the nominal (uncorrected) reconstruction while “gradwarp” imparts marked spatial blurring (especially for the 1.0 and 1.1 mm inserts) and thus resolution loss. Conclusion: We've presented a novel reconstruction strategy for non-Cartesian MRI that corrects for gradient nonlinearities during — rather than after — reconstruction, and thus better preserves image resolution than traditional interpolation-based methods. This approach is expected to be especially advantageous when imaging with non-standard magnet geometries. NIH RR018898; NIH EB10065.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Comment on “Additional Insights Between Fermi-Löwdin Orbital SIC and the Localization Equation Constraints in SIC-DFT”

In a recent Letter, Aquino and Wong (AW) analyze the connection between the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) and the traditional Perdew– Zunger self-interaction correction (PZ-SIC) methods. In this work, the authors present a derivation of the conditions required to minimize the SIC density functional theory (DFT) total energy in FLO-SIC. These involve the gradient of the energy with respect to parameters known as Fermi-orbital descriptors (FODs) (called orbital “centroids” in ref 1) that are used together with the total charge density to determine the Fermi- Löwdin orbitals (FLOs). AW state that the conditions, the vanishing of the FOD gradient components, “resemble” the localization equations (LE), the minimum energy conditions in traditional PZ-SIC. The LE involve the equality of certain off-diagonal matrix elements of different orbital-dependent Hamiltonians. While AW are careful to say that their equations do not prove the mathematical equivalence of the two methods, the discussion in the text and the graphical table of contents figure suggest that vanishing FOD gradients imply satisfaction of the LE and vice versa, and hence that the methods are equivalent. We believe this ambiguity clouds an important issue, rather than clarifying it.

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

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↗