Search NASA⌕ Search

SEARCH · Search NASA

Results for “Series Bosch”

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.

72 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stabilization Of The CN 3 5− Anion In Recoverable High‐pressure Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) Oxoguanidinates

Abstract A series of isostructural Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) oxoguanidinates was synthesized under high‐pressure (25–54 GPa) high‐temperature (2000–3000 K) conditions in laser‐heated diamond anvil cells. The crystal structure of this novel class of compounds was determined via synchrotron single‐crystal X‐ray diffraction (SCXRD) as well as corroborated by X‐ray absorption near edge structure (XANES) measurements and density functional theory (DFT) calculations. The Ln 3 O 2 (CN 3 ) solids are composed of the hitherto unknown CN 3 5− guanidinate anion—deprotonated guanidine. Changes in unit cell volumes and compressibility of Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) compounds are found to be dictated by the lanthanide contraction phenomenon. Decompression experiments show that Ln 3 O 2 (CN 3 ) compounds are recoverable to ambient conditions. The stabilization of the CN 3 5− guanidinate anion at ambient conditions provides new opportunities in inorganic and organic synthetic chemistry.

Chemistry↗
DOE OSTI2026

Are Deep Energy Retrofits in Commercial Buildings Including Window Upgrades?

U.S. Commercial buildings account for about 20% of total U.S. energy consumption. Because the thermal performance of windows significantly affects building energy efficiency and HVAC system performance, best practice guidance often includes window and envelope improvements in conjunction with HVAC upgrades to optimize energy use and improve occupant comfort. It is an open question, however, regarding how often these best practices are implemented in the field. This paper aims to address that gap by conducting a literature review and a series of interviews with commercial building auditing and management professionals to explore the factors that drive window retrofits in commercial buildings. The paper explores a range of case studies from deep energy retrofits across the globe, comparing projects with and without window retrofits. The primary goals of this review are to: (1) provide data from real-world case studies illustrating the role of windows in deep energy renovations and HVAC upgrades, (2) conduct retrofit cost analyses for windows and high-performance HVAC systems and (3) offer insights into how window upgrade decisions are made and when they are implemented as part of deep energy retrofits. Most of the retrofit studies focused exclusively on high performance HVAC upgrades without considering how window upgrades might further enhance the overall energy efficiency of commercial buildings. Interviews with building industry experts shed light on the key factors influencing deep energy retrofit decisions and what factors tip the scales in favor of including window measures with more comprehensive retrofit projects.

Cort, Katherine↗
DOE OSTI2026

Phase Stability and Electrochemical Performance of La-Site-Doped Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 High-Entropy Garnets

We investigate La-site substitution in the high-entropy garnet Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 (LLZNTH) using Ba2+, Sr2+, and Sm3+ to elucidate how dopant governs phase stability, Li-site distribution, and electrochemical behavior. X-ray diffraction shows that Sr2+ is incorporated homogeneously into the garnet lattice, whereas the larger Ba2+ and smaller Sm3+ ions partially exceed the structural tolerance, generating secondary phases. Nevertheless, the Sm-doped composition (x = 0.05) exhibits the highest room-temperature ionic conductivity (2.7 × 10–4 S cm–1). Neutron powder diffraction reveals that Sm substitution drives a redistribution of Li+ from the tetrahedral 24 d sites into the higher-mobility 96 h positions, enhancing the connectivity of the three-dimensional Li-ion migration network. A Sm-doping series (x = 0.01–0.05) further shows that only sufficiently high Sm levels induce this redistribution, whereas lower concentrations retain Li arrangements similar to the undoped garnet. Critical current density measurements demonstrate that La-site dopants also influence interfacial stability against Li metal, underscoring a trade-off between bulk transport enhancement and mechanical robustness. Collectively, these findings reveal that in high-entropy garnets improved ionic conductivity can originate not only from phase-pure structures but also from targeted modification of the Li sublattice, even when accompanied by secondary phases, offering a compositional design principle for garnet electrolytes.

Li, Chang [Mechanical Engineering, School of Scien↗
DOE OSTI2026

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗
DOE OSTI2023

High‐Entropy Lithium Argyrodite Solid Electrolytes Enabling Stable All‐Solid‐State Batteries

Abstract Superionic solid electrolytes (SEs) are essential for bulk‐type solid‐state battery (SSB) applications. Multicomponent SEs are recently attracting attention for their favorable charge‐transport properties, however a thorough understanding of how configurational entropy (ΔS conf ) affects ionic conductivity is lacking. Here, we successfully synthesized a series of halogen‐rich lithium argyrodites with the general formula Li 5.5 PS 4.5 Cl x Br 1.5‐x (0≤x≤1.5). Using neutron powder diffraction and 31 P magic‐angle spinning nuclear magnetic resonance spectroscopy, the S 2− /Cl − /Br − occupancy on the anion sublattice was quantitatively analyzed. We show that disorder positively affects Li‐ion dynamics, leading to a room‐temperature ionic conductivity of 22.7 mS cm −1 (9.6 mS cm −1 in cold‐pressed state) for Li 5.5 PS 4.5 Cl 0.8 Br 0.7 (ΔS conf =1.98R). To the best of our knowledge, this is the first experimental evidence that configurational entropy of the anion sublattice correlates with ion mobility. Our results indicate the possibility of improving ionic conductivity in ceramic ion conductors by tailoring the degree of compositional complexity. Moreover, the Li 5.5 PS 4.5 Cl 0.8 Br 0.7 SE allowed for stable cycling of single‐crystal LiNi 0.9 Co 0.06 Mn 0.04 O 2 (s‐NCM90) composite cathodes in SSB cells, emphasizing that dual‐substituted lithium argyrodites hold great promise in enabling high‐performance electrochemical energy storage.

Chemistry↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

Open quantum system approach to inclusive jet production in heavy-ion collisions

We derive a factorization formula for inclusive jet production in heavy-ion collisions using the tools of Effective Field Theory (EFT). We show how physics at widely separated scales in this process can be systematically separated by matching to EFTs at successively lower virtualities. Owing to a strong scale separation, we recover a vacuum-like DGLAP evolution above the jet scale, while the additional low-energy scales induced by the medium effectively probe the internal structure of the jet. As a result, the cross section can be written as a series with an increasing number of subjets characterized by perturbative matching coefficients each of which is convolved with a distinct function. These functions encode broadening, medium-induced radiations as well as quantum interference such as the Landau-Pomeranchuk-Migdal effect and color coherence dynamics to all orders in perturbation theory. As a first application of this EFT framework, we investigate the case of an unresolved jet and show how the cross section can be factorized and fully separate the jet dynamics from the universal physics of the medium. To compare to the existing literature, we explicitly compute the medium jet function at next-to-leading order in the coupling and leading order in medium opacity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Enhanced Optical Contrast and Switching in Near‐Infrared Electrochromic Devices by Optimizing Conjugated Polymer Oligo(Ethylene Glycol) Sidechain Content and Gel Electrolyte Composition

Abstract A detailed investigation addressing the effects of functionalizing conjugated polymers with oligo(ethylene glycol) (EG n ) sidechains on the performance and polymer‐electrolyte compatibility of electrochromic devices (ECDs) is reported. The electrochemistry for a series of donor‐acceptor copolymers having near‐infrared (NIR)‐optical absorption, where the donor fragment is 3,4‐ethylenedioxythiophene (EDOT) or an EG n functionalized bithiophene (g2T) and the acceptor fragment is diketopyrrolopyrrole (DPP) functionalized with branched alkyl or EG n sidechains, is extensively probed. ECDs are next fabricated and it is found that EG n sidechain incorporation must be finely balanced to promote polymer‐electrolyte compatibility and provide efficient ion exchange. Proper electrolyte‐cation pairing and polymer structural tuning affords a 2x increase in optical contrast (from 12% to 24%) and >60x reduction in switching time (from 20 to 0.3 s). Atomic force microscopy (AFM)/grazing incidence wide‐angle X‐ray scattering (GIWAXS) characterization of the polymer film morphology/microstructure reveals that an over‐abundance of EG n sidechains generates large polymer crystallites, which can suppress ion exchange. Lastly, time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) indicates sidechain/electrolyte identity does not influence the electrolyte penetration depth into the films, and EG n sidechain inclusion increases electrolyte cation uptake. The material structural design insight and guidelines regarding the polymer‐electrolyte ion insertion/expulsion dynamics reported here should be of significant utility for developing next‐generation mixed ionic‐electronic conducting materials.

Chemistry↗
DOE OSTI2026

Adapter Signaling Evaluations on EVs [Slides]

In case J3400 and J3400/1 define a basic analog signaling approach for DC charging adapters to communicate over-temperature events to both the EV and EVSE. They also require the EV to implement mitigation and corrective actions in cases where the EVSE does not respond to adapter thermal signals. In this study, multiple production vehicles were evaluated for responsiveness to assess field readiness and ensure consistent performance in accordance with the SAE J1772, J3400 and J3400/1 standards. A series of test cases were developed and executed across different vehicles, using an adapter thermal breakout fixture designed by NLR for evaluations. The results revealed notable variations in behavior: some vehicles expected to comply with J1772 were found to be non-responsive under certain conditions. Additionally, while some EV OEMs appear to align with the J3400 standard, discrepancies exist due to its evolving nature, resulting in inconsistent implementation of the latest requirements. These findings highlight the need for alignment within standards organizations to ensure consistent interpretation, encourage compliance, and reduce potential confusion across implementations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Fingerprinting Uranium Oxides with Electron Energy Loss Spectroscopy Supported by Theoretical Computations

Uranium oxides occur in a variety of phases that differ in crystal structure and uranium oxidation states. Electron energy loss spectroscopy (EELS) is one of the few techniques that has sufficient spatial resolution and sensitivity to electronic structure to distinguish amongst phases at the nanoscale. However, beam-sensitive materials such as uranium oxides are subject to spectral modification due to interactions with the electron beam. Therefore, theory support is essential to reliably exclude the impact of beam damage and generate true reference datasets. Here we use a comparison of theoretical and experimental spectra to probe the impact of beam damage on O K-edge and U N-edge (N6,7 and N4,5) EELS spectra of various single-valent and mixed-valence uranium oxide bulk phases. Using a low-dose experimental set-up, we show that the O K-edge theoretical spectra are in excellent agreement with experiment for both peak positions and relative intensities of respective peaks. In contrast, U N-edge features are less distinguishing due to the partially localized nature of the U 5f orbitals and overlapping multiplet and spin–orbit coupling effects. This work demonstrates that O K-edge EELS is sufficiently diagnostic to distinguish a wide range of uranium oxides and that the experimental approach used here minimizes beam damage and allows valence state discrimination across the U(IV), U(V) and U(VI) series. When combined with imaging modes available in electron mi-croscopy, the work enables detailed investigation and characterization of uranium redox transformations at the nanoscale.

Carbone, Jacopo↗
DOE OSTI2025

Tunable Metasurface External-Cavity Quantum Cascade Lasers up to 5.74 THz

Vertical-external-cavity surface-emitting lasers based on amplifying quantum-cascade metasurfaces are demonstrated in the 5–6 THz range for the first time. Enhanced parasitic coupling to lossy surface modes requires updated design guidelines, while elevated losses and lower available material gain limit the scaling of the metasurface period from previous designs. A series of metasurface devices with varying periods employing both uniform and focusing metasurfaces is fabricated and characterized. Reducing the metasurface period below 70% of the free-space wavelength enables devices with superior performance, achieving pulsed operation up to 5.74 THz with peak output powers of 1.3 mW, maximum operating temperatures up to 83 K, and continuous-wave operation up to 54 K with 23 μW of output power. In the best case, single-mode tuning from 5.23–5.73 THz, which corresponds to a 9.1% fractional tuning, is realized with near-Gaussian far-field profiles maintained across the entire range. Finally, these results establish quantum cascade vertical-external-cavity surface-emitting lasers as promising candidates for high-frequency local oscillator applications in heterodyne spectroscopy.

gallium arsenide↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2026

Microstructure-dependent particulate populations in ODS steels revealed by single-particle ICP-TOF-MS

In this work, we evaluated the ability of single-particle inductively coupled plasma time-of-flight mass spectrometry (SP-ICP-TOF-MS) to characterize particulates generated from 300 series stainless steels and oxide dispersion–strengthened (ODS) steels. Particulates were produced using two different rotary tools (cutting disk and grinding bit) to assess the influence of particle generation on effluent particulates. Across all materials and preparation techniques, SP-ICP-TOF-MS detected a substantial population of particles containing the expected constituent elements, with average compositions closely matching the nominal alloy values. A secondary population of Y-bearing particles was evident in ODS steels; their formation is strongly favored in particulates generated using the grinding bit. Specifically, in PM2000, Y-bearing particles were predominantly colocated with Al, whereas in 14YWT, Y-bearing particles were less frequently Fe-free but showed stronger colocation with Ti. These trends are consistent with the known microstructural differences between Y–Al–O oxide particles in PM2000 and finer Y–Ti–O nanoclusters in 14YWT. Together, these results demonstrate that SP-ICP-TOF-MS can resolve distinct particulate populations in complex, multicomponent steels at the femtogram mass scale, enabling particle-resolved insight into microstructure-driven phase separation and preparation-dependent particle generation.

Brubaker, Zach E. [Oak Ridge National Laboratory (↗
DOE OSTI2026

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗
DOE OSTI2026

Spectrograph stabilization using a single-delay interferometer on the Hale Telescope

We describe a technique for spectrograph stabilization useful when conventional mitigation techniques of vacuum tanks, thermal insulation, and laser frequency comb may be impractical, expensive, heavy, or bulky. This includes spectrographs on airborne platforms or mounted on telescopes where they suffer a changing gravity vector or other drifts. Placing a fixed-delay interferometer in series with a spectrograph forms an externally dispersed interferometer (EDI). This produces a uniform sinusoidal comb multiplying input spectrum, creating (through heterodyning) beats (moiré patterns). In Fourier space for low frequencies up to the comb frequency, the moiré generated signal counter-rotates to ordinary spectra under an unknown disperser wavenumber drift Δx. This generates a large negative feedback signal useful in a conceptual control loop, to converge rapidly to a stable spectrum and yield Δx. A modified EDI data analysis algorithm (“crossfading”) combines frequency-weighted moiré with conventional spectrum to cancel net output spectrum reaction to Δx. Needing only a single-delay, this is a practical improvement over prior crossfading analyses requiring multiple delays. We test crossfading on ThAr data near 4850 cm−1 taken on Hale telescope in an earlier project. In a single pass, we reduce drift 20 times. Using seven iterations, we reduce 0.5 cm−1 (31 km/s Doppler equivalent) drift to 4×10−7 cm−1 (2.5 cm/s). The interferometer delay can wander, because linearity of phase versus wavenumber interpolates science features between bracketing calibrating spectral references. Second, mathematically reversing the heterodyning effect doubles effective spectral resolution without changing disperser slit.

Erskine, David J [Lawrence Livermore National Labo↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗