Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “working fluid characterization”

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.

616 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Permanently Porous Chalcogen-Bonded Organic Framework

The nature of connectivity between constituent atomic or molecular building blocks is fundamental in shaping the properties and functionality of materials. The extrapolation of emergent interatomic interactions to enable functional materials has driven transformative technological advancements. However, the bonding interactions used in material design have been largely static since the emergence of dynamic covalent chemistry ~30 years ago. Here we demonstrate that non-covalent chalcogen bonding (Ch-bonding) is a distinct mode of interatomic connectivity for constructing functional materials by design. This is established by leveraging self-complementary assembly of 1,2,5-telluradiazole moieties to construct a honeycomb-type permanently porous Ch-bonded organic framework, assembled and stabilized solely through non-covalent Te...N contacts. Empirical and computational studies of electronic structure, structural healing and lattice dynamics highlight the pi-type electronic communication, controlled assembly and modulated lattice dynamics in Trip3Tez-I arising directly from the unique nature of the Te...N Ch-bonding that holds substantial implications for next generation crystalline semiconductors. In addition to introducing a distinct class of permanently porous frameworks, this work establishes Ch-bonding as a programmable molecular tool for constructing functional materials with distinct properties.

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

Kinetic control of phase evolution and defect-mediated recombination in AACVD-grown copper antimony sulphide thin films

Copper antimony sulphide (CAS) is a multinary chalcogenide semiconductor in which small deviations from stoichiometry can drive phase competition and strong defect-mediated modulation of optoelectronic properties. However, the systematic roles of copper precursor fraction and deposition time in governing phase evolution, off-stoichiometry, and recombination dynamics in aerosol-assisted chemical vapour deposition (AACVD)-grown undoped CAS thin films remain insufficiently understood. In this work, undoped CAS thin films were deposited by AACVD using Cu(dedtc)2 and Sb(dedtc)3 single-source precursors at 550 °C and a carrier gas flow rate of 150 sccm, while the Cu(dedtc)2 mole fraction (x = 0.15 – 0.55) and deposition time (1 – 2 hours) were systematically varied to probe how growth kinetics influence phase composition, microstructure, and defect-mediated optical properties without post-deposition annealing or extrinsic doping. Increasing copper precursor content drives phase evolution toward tetrahedrite-dominant CAS films at intermediate Cu(dedtc)2 mole fractions, with the film deposited at x = 0.35 exhibiting the strongest tetrahedrite character within the parameter space examined. The films are also copper-rich, antimony-poor, and sulphur-deficient, consistent with off-stoichiometric growth and intrinsic defect formation, plausibly including copper interstitials, Cu-on-Sb antisites, and sulphur vacancies. These growth-dependent compositional deviations are accompanied by tunable indirect optical bandgaps of approximately 1.60 – 2.18 eV and weak visible photoluminescence governed by defect-mediated recombination. Time-resolved photoluminescence reveals bi-exponential decay behaviour with lifetimes of approximately 0.1 – 3.6 ns, with emission dominated by slower donor–acceptor pair recombination and a smaller contribution from faster trap-assisted pathways. Collectively, these results establish an explicit kinetic process–structure–defect–property relationship for AACVD-grown CAS thin films and provide a growth–structure–property framework relevant to future optimization of CAS-based optoelectronic and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

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

PACT Center: Perovskite PV Accelerator for Commercializing Technologies (Final Technical Report)

The Perovskite PV Accelerator for Commercializing Technologies (PACT) center was established in July 2021 as a national resource to accelerate the commercialization of perovskite photovoltaic (PV) technology in the United States. Since its inception, PACT has been led by Sandia National Laboratories (Sandia) in partnership with the National Laboratory of the Rockies (NLR), formerly known as NREL. From FY20-FY23, Los Alamos National Laboratory (LANL), CFV Labs, Black & Veatch (B&V), and the Electric Power Research Institute (EPRI) were part of the project team. LANL brought expertise in perovskite PV device designs and processing, CFV Labs (now GroundWork Renewables) provided initial indoor and outdoor measurement hardware technology, B&V led the initial effort on perovskite PV bankability, and EPRI worked on reviewing testing standards, identifying commercialization gaps, and helping to run PACT’s Industry Advisory Board, a group including representatives from commercial testing labs, independent engineering firms, insurance companies, state regulators, and electric utilities. To source perovskite PV module samples, PACT contracted with the University of North Carolina (UNC), the University of Toledo, the University of Washington, and SLAC/Stanford University to provide a steady stream of research-grade perovskite mini modules, enabling protocol development in advance of commercial module availability. The project period ran from July 1, 2021, through December 31, 2025, including a No Cost Extension. Starting in FY25, the project was continued as a Core Capability in the Lab Call portfolio and continues at a reduced budget with only Sandia and NLR as funded recipients. Notably, starting in FY25 PACT expanded its scope beyond MHP modules to accept all emerging PV mini module technologies for testing, including organic PV (OPV) and all-thin-film tandems, with the aim of supporting commercialization across the broader emerging PV ecosystem. With this change in scope the program was renamed the PV Accelerator for Commercializing Technologies, dropping perovskite from the name.

14 SOLAR ENERGY↗
DOE OSTI2026

Maintenance strategy, structural design, and site layout of the ST-E1 fusion power plant

An effective fusion reactor maintenance scheme enables safe operations and short downtimes. This in turn leads to high availability, which is critical to the commercial viability of a power-producing plant. In tokamak-based fusion power plants, the chosen maintenance approach has a significant impact on the spatial design of the tokamak, as well as the surrounding infrastructure, and therefore needs to be considered from the outset. Tokamak Energy has developed a pre-concept design of a fusion power plant, ST-E1. This work describes the major drivers and constraints that have been considered, presents the tokamak architecture and chosen maintenance regime, and discusses how this enables the plant’s two-phased approach to demonstrating commercial operations. It also shows the implications for the design of other systems areas, in particular the machine structural arrangement and bioshield and hot cell layout. The reactor core segmentation and removal scheme replaces entire toroidal segments radially through a large vacuum port, along a single axis only. The result is a change-tolerant machine and plant layout that can accommodate the evolving designs of the tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Landmark Publications in Analytical Atomic Spectrometry: Fundamentals and Instrumentation Development

The almost-two-centuries history of spectrochemical analysis has generated a body of literature so vast that it has become nearly intractable for experts, much less for those wishing to enter the field. Authoritative, focused reviews help to address this problem but become so granular that the overall directions of the field are lost. This broader perspective can be provided partially by general overviews but then the thinking, experimental details, theoretical underpinnings, and instrumental innovations of the original work must be sacrificed. In the present compilation, this dilemma is overcome by assembling the most impactful publications in the area of analytical atomic spectrometry. Each entry was proposed by at least one current expert in the field and supported by a narrative that justifies its inclusion. The entries were then assembled into a coherent sequence and returned to contributors for a round-robin review. A total of 48 scientists participated in this endeavor, contributing a combined list of 1055 individual articles spanning 17 sub-disciplines of spectrochemical analysis into what the current community views as “key” publications. Of these cited articles, 60 received nominations from four or more scientists, establishing them as the most indispensable reading materials. The outcome of this collaborative effort is intended to serve as a valuable resource not only for current practitioners in atomic spectroscopy but also for present and future students who represent coming generations of analytical atomic spectroscopists.

Instruments & Instrumentation↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2024

Lanthanum-Promoted Electrocatalyst for the Oxygen Evolution Reaction: Unique Catalyst or Oxide Deconstruction?

A conventional performance metric for electrocatalysts that promote the oxygen evolution reaction (OER) is the current density at a given overpotential. However, the assumption that increased current density at lower overpotentials indicates superior catalyst design is precarious for OER catalysts in the working environment, as the crystalline lattice is prone to deconstruction and amorphization, thus greatly increasing the concentration of catalytic active sites. We show this to be the case for La 3+ incorporation into Co 3 O 4 . Powder X-ray diffraction (PXRD), Raman spectroscopy and extended X-ray absorption fine structure (EXAFS) reveal smaller domain sizes with decreased long-range order and increased amorphization for La-modified Co 3 O 4 . This lattice deconstruction is exacerbated under the conditions of OER as indicated by operando spectroscopies. The overpotential for OER decreases with increasing La 3+ concentration, with maximum activity achieved at 17% La incorporation. HRTEM images and electron diffraction patterns clearly show the formation of an amorphous overlayer during OER catalysis that is accelerated with La 3+ addition. O 1s XPS spectra after OER show the loss of lattice-oxide and an increase in peak intensities associated with hydroxylated or defective O-atom environments, consistent with Co(O) x (OH) y species in an amorphous overlayer. Furthermore, our results suggest that improved catalytic activity of oxides incorporated with La 3+ ions (and likely other metal ions) is due to an increase in the number of terminal octahedral Co(O) x (OH) y edge sites upon Co 3 O 4 lattice deconstruction, rather than enhanced intrinsic catalysis.

Catalysts↗
DOE OSTI2026

Studies of laser stimulated photodetachment from nanoparticles for particle charge measurements

Determining nanoparticle charge is more challenging than that for microparticles due to change in the particle size during the synthesis substantial plasma property variations, and difficulties in visualizing individual particles, rendering conventional microparticle charge diagnostics ineffective in dusty plasma. In this work, we utilized laser-stimulated photodetachment (LSPD) to deduce the mean charge of nanoparticles. Nanoparticles were grown in an Ar/C 2 H 2 mixture using a capacitively coupled RF discharge and the LSPD induced changes in the electron current monitored by a cylindrical Langmuir probe. LSPD signals were obtained and analyzed across different dust growth phases. The prolonged decay of electron current pulses was attributed to the presence of residual negative ions, caused by the effective electrostatic trapping of these ions and the potential post—LSPD re-formation of new ones. The particle charge was estimated by combining the laser-stimulated photodetachment signal from the probe with the dust density obtained from laser-light extinction using the measured nanoparticle size distribution. For a nanoparticles size range of approximately 100–250 nm and mean diameter of $d$ p ∼154.37 nm, the effective mean charge was estimated to be $\langle$$Q$$\rangle$ d $≈$ 37 elementary charge units. The measured charge values are lower than those predicted by orbital motion limited theory, which may be attributed to significant electron depletion in the nanodusty plasma. LSPD results in Ar/C 2 H 2 nano-dusty plasma confirm the applicability of this method for estimating individual nanoparticle charges. However, it has also been demonstrated that electron detachment from residual background negative ions can influence the detachment current decay and must be carefully considered.

dust particle charge↗
DOE OSTI2026

Boron Coordination in Multicomponent Glasses: Analytical Models and Machine Learning With Uncertainty

Borosilicate glasses are extensively used in a variety of applications from kitchenware to nuclear waste immobilization due to the strong network formed by the Si-O-B bond that makes it resistant to chemical corrosion and gives it a low thermal expansion. Boron, however, exists in both trigonal BO3 and tetrahedral BO4 bonds in glass systems, which impacts the chemical durability and thermal resistance of the glass, amongst other properties. Boron coordination (N4), or the ratio of the amount of BO4 to BO3 within a glass, may aid in predicting these properties but is difficult to derive without experimental data due to the complexity of impacts from varied glass compositions and processing factors. For this reason, compositional models have been developed to predict boron coordination, but the models typically include a limited number of glass components. To help fill this gap in the models, in this work, a diverse multicomponent glass dataset of 809 glasses is compiled from a literature search, and then a number of analytical and machine learning (ML) models are trained on the dataset. Previously developed modified Bernstein and modified Du Stebbins analytical models were fitted to update parameters with the new dataset. Then, partially Bayesian neural networks, Gaussian process regressor, and heteroskedastic deterministic neural networks were evaluated. The ML models examined all have different strategies to overcome the potential for overfitting as a result of a limited training dataset, and return results that account for model uncertainty, which can be valuable for understanding model reliability. For the first time, cooling rate is introduced as an input parameter for ML models, showing consistent improvements in performance and solidifying the importance of including parameters outside of composition alone for N4 prediction. The machine learning models examined here show promise in accurate predictions of boron coordination in borosilicate glasses, all achieving R2 values of 0.91.

boron coordination↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2025

Tunable high Néel temperature and large anomalous Hall response in antiferromagnetic Weyl semimetal Mn 3 Sn 1− x Ga x thin films

Antiferromagnetic Weyl semimetals based on Mn 3 X(X = Ge, Sn, Ga) kagome compounds exhibit the same ferromagnetic-like responses, including anomalous Hall, Nernst, and magneto-optical effects, as recently discussed for altermagnets. Driven by the Berry curvature due to Weyl fermions, these materials show a disproportionately large magnitude of electromagnetic effects even in the absence of large magnetization. For applications it is crucial to realize these responses in a wide range of temperatures both below and above 300 K. While stoichiometric Mn 3 X materials do not offer optimal performance, we show that Mn 3 Sn 1−x Ga x sputtered films with a variable composition offers a tunable Néel temperature, T N ≈ 425 ± 6–500 ± 15 K, which is crucial for device applications, together with a large tunable anomalous Hall effect. Our thin film growth method enables continuous and precise control over the film composition between x = 0 and x = 1. Through a detailed magnetization and Hall transport, we establish the magnetic phase diagram for the hexagonal Mn 3 Sn 1−x Ga x . Our results reveal an enhanced T N and antichiral magnetic phase in Ga-doped Mn 3 Sn and an enhanced anomalous Hall magnitude in Sn-doped Mn 3 Ga compared to their stoichiometric undoped forms. Our work demonstrates a route to optimize the technologically relevant antiferromagnets for various applications.

Magnetic properties and materials↗
DOE OSTI2023

Design and Synthesis of Cubic K 3−2 x Ba x SbSe 4 Solid Electrolytes for K–O 2 Batteries

Developing K-ion conducting solid-state electrolytes (SSEs) plays a critical role in the safe implementation of potassium batteries. In this work, a chalcogenide-based potassium ion SSE is reported, K 3 SbSe 4 , which adopts a trigonal structure at room temperature. Single-crystal structural analysis reveals a trigonal-to-cubic phase transition at the low temperature of 50 °C, which is the lowest among similar compounds and thus provides easy access to the cubic phase. The substitution of barium for potassium in K 3 SbSe 4 leads to the creation of potassium vacancies, expansion of lattice parameters, and a transformation from a trigonal phase to a cubic phase. As a result, the maximum conductivity of K 3−2x Ba x SbSe 4 reaches around 0.1 mS cm −1 at 40 °C for K 2.2 Ba 0.4 SbSe 4 , which is over two orders of magnitude higher than that of undoped K 3 SbSe 4 . This novel SSE is successfully employed in a K–O 2 battery operating at room temperature where a polymer-laminated K 2.2 Ba 0.4 SbSe 4 pellet serves as a separator between the oxygen cathode and the potassium metal anode. Effective protection of the K metal anode against corrosion caused by O 2 is demonstrated.

25 ENERGY STORAGE↗
DOE OSTI2024

The Active Asteroids Citizen Science Program: Overview and First Results

Abstract We present the Citizen Science program Active Asteroids and describe discoveries stemming from our ongoing project. Our NASA Partner program is hosted on the Zooniverse online platform and launched on 2021 August 31, with the goal of engaging the community in the search for active asteroids—asteroids with comet-like tails or comae. We also set out to identify other unusual active solar system objects, such as active Centaurs, active quasi-Hilda asteroids (QHAs), and Jupiter-family comets (JFCs). Active objects are rare in large part because they are difficult to identify, so we ask volunteers to assist us in searching for active bodies in our collection of millions of images of known minor planets. We produced these cutout images with our project pipeline that makes use of publicly available Dark Energy Camera data. Since the project launch, roughly 8300 volunteers have scrutinized some 430,000 images to great effect, which we describe in this work. In total, we have identified previously unknown activity on 15 asteroids, plus one Centaur, that were thought to be asteroidal (i.e., inactive). Of the asteroids, we classify four as active QHAs, seven as JFCs, and four as active asteroids, consisting of one main-belt comet (MBC) and three MBC candidates. We also include our findings concerning known active objects that our program facilitated, an unanticipated avenue of scientific discovery. These include discovering activity occurring during an orbital epoch for which objects were not known to be active, and the reclassification of objects based on our dynamical analyses.

Astronomy & Astrophysics↗
DOE OSTI2025

Direct observation of distinct bulk and edge nonequilibrium spin accumulation in ultrathin MoTe 2

Low-symmetry two-dimensional (2D) topological materials such as MoTe 2 host efficient charge-to-spin conversion (CSC) mechanisms that can be harnessed for novel electronic and spintronic devices. However, the nature of the various CSC mechanisms and their correlation with underlying crystal symmetries remain unsettled. In this work, we use local spin-sensitive electrochemical potential measurements to directly probe the spatially dependent nonequilibrium spin accumulation in MoTe 2 flakes down to four atomic layers. We are able to clearly disentangle contributions originating from the spin Hall and Rashba-Edelstein effects and uncover an abundance of unconventional spin polarizations that develop uniquely in the sample bulk and edges with decreasing thickness. Using ab-initio calculations, we construct a unified understanding of all the observed CSC components in relation to the material dimensionality and stacking arrangement. Our findings not only illuminate previous CSC results on MoTe 2 but also have important ramifications for future devices that can exploit the local and layer-dependent spin properties of this 2D topological material.

Science & Technology - Other Topics↗
DOE OSTI2025

Considerations for Maritime Nuclear Technologies, Economic Viability and Public Acceptance

The Maritime Nuclear Application Group (MNAG) is a working group convened by the National Reactor Innovation Center at Idaho National Laboratory (INL), the American Bureau of Shipping, and Morgan, Lewis, and Bockius LLP. This report documents an MNAG examination of considerations relevant to implementing nuclear technology in commercial maritime applications. In general, two types of use case are examined: maritime nuclear power plants and nuclear reactors used on board shipping vessels for propulsion and other ship needs. The report finds that there may be economic benefits related to maritime nuclear technologies, including the flexible deployment of maritime nuclear reactors, which would allow them to complement land-based nuclear projects, and operational differences for nuclear cargo ships that may lead to an overall increase in revenue. High-level analyses in this report show that, based on general small modular reactor and microreactor cost estimates developed by INL, maritime nuclear reactors may be economically competitive for electricity production in remote regions and for use in the propulsion of large cargo ships. Besides economic viability, public acceptance will be key to implementing maritime nuclear technologies. The report discusses the public’s current perception of nuclear technologies. Engaging with the public will be important to improving this perception. The report discusses some key benefits and risks associated with maritime nuclear technologies. Benefits include the creation of jobs, the production of reliable energy, and the potential to improve air quality. Risks that concern the public are the potential for radioactive releases during operation and decommissioning, as well as those related to waste management. Communicating the benefits and the risks of maritime nuclear technologies will be essential to improving public perception.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
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↗