Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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

Dark Energy Survey Year 3 results: galaxy–halo connection from galaxy–galaxy lensing

ABSTRACT Galaxy–galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter haloes, which is important both for galaxy evolution and cosmology. We extend the measurement and modelling of the galaxy–galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly non-linear scales (∼100 kpc). This extension enables us to study the galaxy–halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redmagic) and a magnitude-limited galaxy sample (maglim). We find that redmagic (maglim) galaxies typically live in dark matter haloes of mass log10(Mh/M⊙) ≈ 13.7 which is roughly constant over redshift (13.3−13.5 depending on redshift). We constrain these masses to ${\sim}15{{\ \rm per\ cent}}$, approximately 1.5 times improvement over the previous work. We also constrain the linear galaxy bias more than five times better than what is inferred by the cosmological scales only. We find the satellite fraction for redmagic (maglim) to be ∼0.1−0.2 (0.1−0.3) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints, and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Scalable Approach to Minimize Charging Costs for Electric Bus Fleets

Incorporating battery electric buses into bus fleets faces three primary challenges: a BEB’s extended refuel time, the cost of charging, both by the consumer and the power provider, and large compute demands for planning methods. When BEBs charge, the additional demands on the grid may exceed hardware limitations, so power providers divide a consumer’s energy needs into separate meters even though doing so is expensive for both power providers and consumers. Prior work has developed a number of strategies for computing charge schedules for bus fleets; however, prior work has not worked to reduce costs by aggregating meters. Additionally, because many works use mixed integer linear programs, their compute needs make planning for commercial-sized bus fleets intractable. This work presents a multi-program approach to computing charge plans for electric bus fleets. The proposed method solves a series of subproblems where the solution to the charge problem becomes more refined with each problem, moving closer to the optimal schedule. The results demonstrate how runtimes are reduced by using intermediate subproblems to refine the bus charge solution so that the proposed method can be applied to large bus fleets of 100+ buses. Not only will we demonstrate that runtimes scale linearly with the number of buses but we will also show how the proposed method scales to large bus fleets of over 100 buses while managing the monthly cost of energy.

Mortensen, Daniel (ORCID:0000000276494452)↗
DOE OSTI2022

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the growth rate of structure from the small-scale clustering of the luminous red galaxy sample

ABSTRACT We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\, h^{-1}{\rm Mpc}$ with a model based on the aemulus cosmological emulator to measure the growth rate of cosmic structure, parametrized by fσ8. We obtain a measurement of fσ8(z = 0.737) = 0.408 ± 0.038, which is 1.4σ lower than the value expected from 2018 Planck data for a flat ΛCDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1\text{--}60\, h^{-1}{\rm Mpc}$ modelled by the aemulus cosmological emulator and find a 4.5σ tension in the amplitude of the halo velocity field with the Planck + ΛCDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

Using Graph Coloring to Compute Total Derivatives More Efficiently in OpenMDAO

When they are applicable, gradient based optimization algorithms are the most efficient way to solve design optimization problems. Although gradient based methods are generally efficient, they can be made significantly more so through the usage of analytic techniques to compute the necessary total derivatives. The traditional forward (direct) and reverse (adjoint) analytic techniques have computational costs that scale linearly with the number of design variables and the number of constraints, respectively. In this work, we present an application of a graph coloring algorithm to the analytic techniques for computing total derivative Jacobians in order to achieve much better computational scaling than the pure analytic methods can provide alone. A detailed theoretical explanation of how coloring algorithms interact with analytic derivative methods is presented that illustrates specific types of sparsity patterns that must be present in total derivative Jacobians in order for this coloring technique to be effective. The new technique has been implemented as a feature in the OpenMDAO framework and the implementation is demonstrated on two example problems. The performance on the example problems up to 50% reduction in compute cost for optimizations with bi-directional coloring compared to traditional constraint aggregation. Additionally, the results show how coloring technique alleviates some of the numerical difficulties that constraint aggregation can cause, leading to the ability to solve larger problems. It is expected that the new method will have wide applicability to multidisciplinary optimization problems, and that its availability in OpenMDAO will offer significant computational savings for users without the need for them to implement the coloring algorithm themselves.

Trajectory↗
DOE OSTI2023

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Spatiotemporal Downscaling Model for Solar Irradiance Forecast Using Nearest-Neighbor Random Forest and Gaussian Process

Accurate solar photovoltaic (PV) capacity estimation requires high-resolution, site-specific solar irradiance data to account for localized variability. However, global datasets, such as the National Solar Radiation Database (NSRDB), provide regional averages that fail to capture the fine-scale fluctuations critical for large-scale grid integration. This limitation is particularly relevant in the context of increasing distributed energy resources (DERs) penetration, such as rooftop PV. Additionally, it is critical to the implementation of the U.S. Federal Energy Regulatory Commission (FERC) Order 2222, which facilitates DER participation in U.S. bulk power markets. To address this challenge, this study evaluates Nearest-Neighbor Random Forest (NNRF) and Nearest-Neighbor Gaussian Process (NNGP) models for spatiotemporal downscaling of global solar irradiance data. By leveraging historical irradiance and meteorological data, these models incorporate spatial, temporal, and feature-based correlations to enhance local irradiance predictions. The NNRF model, a machine-learning approach, prioritizes computational efficiency and predictive accuracy, while the NNGP model offers a level of interpretability and prediction uncertainty by numerically quantifying correlations and dependencies in the data. Model validation was conducted using day-ahead predictions. The results showed that the average Goodness of Fit (GoF) of the NNRF model of 90.61% across all eight sites outperformed the GoF of the NNGP of 85.88%. Additionally, the computational speed of NNRF was 2.5 times faster than the NNGP. Finally, the NNGP displayed polynomial scaling while the NNRF scaled linearly with increasing number of nearest neighbors. Additional validation of the model on five sites in Puerto Rico further confirmed the superiority of the NNRF model over the NNGP model. These findings highlight the robustness and computational efficiency of NNRF for large-scale solar irradiance downscaling, making it a strong candidate for improving PV capacity estimation and real-time electricity market integration for DERs.

Asiedu, Shadrack (ORCID:0009000646004826)↗
NASA NTRS1999

Effect of a Nonplanar Melt-Solid Interface on Lateral Compositional Distribution during Unidirectional Solidification of a Binary Alloy with a Constant Growth Velocity V: Theory - Part 1

Infrared detector materials, such as Hg(1-x)Cd(x)Te, Hg(1-x)Zn(x)Te have energy gaps almost linearly proportional to its composition. Due to the wide separation of liquidus and solidus curves of their phase diagram, there are compositional segregations in both of axial and radial directions of these crystals grown in the Bridgman system unidirectionally with constant growth rate. It is important to understand the mechanisms which affect lateral segregation such that large uniform radial composition crystal is possible. Following Coriell, etc's treatment, we have developed a theory to study the effect of a curved melt-solid interface shape on the lateral composition distribution. The system is considered to be cylindrical system with azimuthal symmetric with a curved melt-solid interface shape which can be expressed as a linear combination of a series of Bessell's functions. The results show that melt-solid interface shape has a dominate effect on lateral composition distribution of these systems. For small values of b, the solute concentration at the melt-solid interface scales linearly with interface shape with a proportional constant of the product of b and (1 - k), where b = VR/D, with V as growth velocity, R as sample radius, D as diffusion constant and k as distribution constant. A detailed theory will be presented. A computer code has been developed and simulations have been performed and compared with experimental results. These will be published in another paper.

Wang, Jai-Ching↗
NASA NTRS1988

Microwave background anisotropy induced by gravitational waves

A cosmological background of gravitational waves induces redshift perturbations in light transversing it. Calculations of this Sachs-Wolfe effect on the microwave background are presented in an Omega = 1 Friedmann universe as a function of angular scale and gravitational wave spectrum. Blurriness of the last-scattering surface can cause nonnegligible dilution of the anisotropy for wavelengths less than about 100 Mpc. The limit implied for the energy density of the gravitational waves is given. A difficulty in associating a linear scale with an angular anisotropy, due to the clumpiness of the universe, is also pointed out.

Linder, Eric V.↗
NASA NTRS1999

Effect of a Nonplanar Melt-Solid Interface On Lateral Compositional Distribution During Unidirectional Solidification of a Binary Alloy with a Constant Growth Velocity V: Theory - Pt. 1

Infrared detected materials, such as Hg(1-x)Cd(x)Te, Hg(1-x)Zn(x)Te have energy gaps almost linearly proportional to their composition. Due to the wide separation of liquidus and solidus curves of their phase diagram, there are compositional segregation in both of the axial and radial directions of these crystals grown in the Bridgman system unidirectionally with constant growth rate. It is important to understand the mechanisms, which affect lateral segregation such that large radially uniform composition crystal can be produced. Following Coriel, etc's treatment, we have developed a theory to study the effect of a curved melt-solid interface shape on lateral composition distribution. The model is considered to be a cylindrical system with azimuthal symmetry and a curved melt-solid interface shape which can be expressed as a linear combination of a series of Bessell's functions. The results show that melt-solid interface shape has a dominant effect on the lateral composition distribution of these systems. For small values of beta, the solute concentration at the melt-solid interface scales linearly with interface shape with a proportional constant of the produce of beta and (1 -k), where beta = VR/D, with V as growth velocity, R as the sample radius, D as the diffusion constant and k as the distribution constant. A detailed theory will be presented. A computer code has been developed and simulations have been performed and compared with experimental results. These will be published in another paper.

Wang, Jai-Ching↗
NASA NTRS1994

Towards a Modern Theory of Multiphase Filtration Flow

An alternative theoretical model of joint filtration flow of immiscible incompressible fluids is presented. The model takes into account relaxation processes due to the interchange of the fluids between pores of difference sizes which is driven by capillary forces. The fluids occupy connected regions in a four-dimensional space formed by three coordinates and the pore length scale. When the fluid exchange between pores of given sizes is effected by way of successive flow through pores of all the intermediate sizes, the pressure within each region is governed by a hyperbolic equation, the role of time being played by the pore linear scale. Pressure jumps across hypersurfaces separating the regions equal corresponding values of the capillary pressure. A supplementary condition at any such hypersurface requires the speed of its displacement in the four-dimensional space to coincide with the normal velocity components of both the adjoining fluids. As a result, a principally new statement of multiphase filtration flow problems is gained with allowance for capillary relaxation in the porous space.

Buyevich, Yu A.↗
DOE OSTI2025

Co-designing Spectral Transformation Oracles with Hybrid Oscillator-Qubit Quantum Processors: From Algorithms to Compilation

We co-design a family of quantum eigenvalue transformation oracles that can be efficiently implemented on hybrid discrete- or continuous-variable (qubit or qumode) hardware. To illustrate the oracle’s representation-theoretic power and near-term experimental accessibility, we encode a Gaussian imaginary time-evolution spectral filter. As a result, we define a continuous linear combination of unitaries block encoding. Due to the ancillary qumode’s infinite-dimensional nature, continuous-variable qumodes constitute a powerful compilation tool for encoding continuous spectral functions without discretization errors while minimizing resource requirements. We then focus on the ubiquitous task of preparing eigenstates in quantum spin models. For completeness, we provide an end-to-end compilation which expresses high-level oracles in terms of an experimentally realizable instruction set architecture in both 1D and 2D. Finally, we examine the leading-order effects of physical errors and highlight open research directions. Our algorithms scale linearly with the spatial extent of the target system and are applicable to both near-term and large-scale quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Tracer-field cross-correlations with k -nearest neighbour distributions

ABSTRACT In astronomy and cosmology significant effort is devoted to characterizing and understanding spatial cross-correlations between points – e.g galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous fields – e.g. weak lensing meiand Cosmic Microwave Background maps. Recently, we introduced the k-nearest neighbour (kNN) formalism to better characterize the clustering of discrete (point) data sets. Here, we extend it to the point – field cross-correlations analysis. It combines kNN measurements of the point data set with measurements of the field smoothed at many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two data sets even when there are no linear (Gaussian) correlations between them. We further demonstrate that this framework is far more effective than the two point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlation between haloes and the underlying matter field in a cosmological simulation, between 10 h−1 Mpc and 30 h−1 Mpc, is detected at >5σ significance using the technique presented here, when the two-point cross-correlation significance is ∼1σ. Finally, we show that kNN cross-correlations of haloes and the matter field can be well modelled on quasi-linear scales using the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for 2-pt cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations using this method makes it a promising tool for various cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Improved CdTe PLR Estimates: Self-Shading and Spectral Mismatch

The RdTools year-on-year method of estimating performance loss rate (PLR) employs a simple normalization to remove the confounding effect of irradiance and temperature variation. However, the normalization's assumption that PV production scales linearly with in-plane broadband irradiance is a worse approximation for CdTe and other technologies with larger spectral sensitivities than it is for the more common c-Si technology. Additionally, CdTe systems using single-axis trackers (-20% of installed US utility-scale capacity) are subject to self-shading in the morning and afternoon, introducing another nonlinearity between PV output and broadband irradiance. Ignoring these effects may subject the estimated PLR to increased uncertainty, and perhaps bias, depending on the character of their short- and long-term variability. In this work we show that including self-shading and spectral mismatch models in the normalization for tracking CdTe systems can result not only in tighter PLR confidence intervals but different median PLRs as well. The shading and spectral models are kept simple to maintain consistency with the RdTools ethos of not requiring detailed system metadata or unusual measurements.

CdTe↗
DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗
NASA NTRS2012

Electron Beam Induced Luminescence of SiO2 Optical Coatings

Optical coatings of disordered thin film SiO,/SiO, samples undergoing electron beam bombardment exhibited a cathodoluminescence effect. Measurements were made of the absolute illuminance and emission spectra as functions of incident electron energy (0.2 keV to 30 keV). flux (approx. 0.1 nA/sq cm to > 100 nA/sq cm) and power (<50 nW/sq cm to > 3mW/sq cm) for coatings (approx.60 nm to approx.200 nm thick) on reflective metal substrates over a range of sample temperatures (approx.40 K to approx.400 K) and emission wavelengths (approx.260 nm to approx.5000 nm). Illuminance reached a plateau saturation intensity I(sub sat) approx. 10(exp -4)W/sq cm. scaling as [I(sub mc)/(I + I(sub mc)/I(sub sat)] as incident power I(sub mc) approached I(sub sat). Well below I(sub sat), illuminance scaled linearly with incident power for lower-energy non-penetrating radiation that deposited all of its energy in the dielectric film. At higher incident energies (> 1-5 keV) electrons penetrate the dielectric, depositing only a fraction of their energy in the dielectric film, where deposited energy scales as the ratio of film thickness to penetration range. Since electron range exhibits a peak below approx. 1 keV, the illuminance decreased with increasing energy for penetrating radiation. The temperature dependence of the illuminance. in general, increased with decreasing temperature. Changes in apparent color with temperature were also observed. Four bands in the UV/VIS/NIR range from 300 nm to 1000 nm were observed in spectral measurements; no significant luminescence was detected in a range from approx. 1000 nm to approx. 5000 nm. These different bands either increased or decreased in intensity and shifted slightly in energy with decreasing temperature. Our measurements are explained in terms of simple disordered band theory models based on atomic-scale stochastic transport of electrons in highly disordered insulating materials. The theoretical models provide a fundamental basis for understanding the dependence of cathodoluminescence on irradiation time, incident flux and energy. and sample thickness and temperature. The observed luminescence occurs when an incident high energy electron excites a valence band electron into the conduction band. Following the excitation. a complex series of transitions take place between the extended conduction states and the long-lived localized trapped states (or chromophores) associated with structural or compositional defects in the dielectric materials. Thc transport equations and models of thc spatial and energy distribution of the trapped states below the mobility edge, their occupancy, and transitions provide microscopic scale insights into this macroscopic phenomenon. Our observations and theoretical models are also compared with similar previous room temperature studies of similar SiO2 films.

Dennison, J. R.↗
NASA NTRS2011

A Flight Dynamics Model for a Multi-Actuated Flexible Rocket Vehicle

A comprehensive set of motion equations for a multi-actuated flight vehicle is presented. The dynamics are derived from a vector approach that generalizes the classical linear perturbation equations for flexible launch vehicles into a coupled three-dimensional model. The effects of nozzle and aerosurface inertial coupling, sloshing propellant, and elasticity are incorporated without restrictions on the position, orientation, or number of model elements. The present formulation is well suited to matrix implementation for large-scale linear stability and sensitivity analysis and is also shown to be extensible to nonlinear time-domain simulation through the application of a special form of Lagrange s equations in quasi-coordinates. The model is validated through frequency-domain response comparison with a high-fidelity planar implementation.

Orr, Jeb S.↗
DOE OSTI2023

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗