Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Graphical Model for Fusing Diverse Microbiome Data

This paper develops a Bayesian graphical model for fusing disparate types of count data. The motivating application is the study of bacterial communities from diverse high-dimensional features, in this case, transcripts, collected from different treatments. In such datasets, there are no explicit correspondences between the communities and each corresponds to different factors, making data fusion challenging. We introduce a flexible multinomial-Gaussian generative model for jointly modeling such count data. This latent variable model jointly characterizes the observed data through a common multivariate Gaussian latent space that parameterizes the set of multinomial probabilities of the transcriptome counts. The covariance matrix of the latent variables induces a covariance matrix of co-dependencies between all the transcripts, effectively fusing multiple data sources. We present a computationally scalable variational Expectation-Maximization (EM) algorithm for inferring the latent variables and the parameters of the model. Here, the inferred latent variables provide a common dimensionality reduction for visualizing the data and the inferred parameters provide a predictive posterior distribution. In addition to simulation studies that demonstrate the variational EM procedure, we apply our model to a bacterial microbiome dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Development of a Model Correlating Axial HEPA Filter Performance Between Elevated and Standard Conditions

The ASME AG-1 Code specifies minimum rated airflow based on standard cubic feet per minute (SCFM) for axial high efficiency particulate air (HEPA) filters and actual cubic feet per minute (ACFM) for radial HEPA filters with corresponding maximum pressure drop. A quantitative analysis can be performed to demonstrate the variation in filter efficiency and develop a correlation of performance between axial filters subjected to ACFM and SCFM. This study illustrates the variation in filter efficiency of axial HEPA filters subjected to ACFM and SCFM flow rate conditions at varying atmospheric pressures. It also presents the methods used to calculate correction factors to convert filter efficiency and pressure drop from elevated to standard conditions. The Axial Large Scale Test Stand (ALSTS) in the Institute for Clean Energy Technology (ICET) at Mississippi State University has been utilized to control aerosol concentration, flowrate, atmospheric pressure, relative humidity, and temperature. The filter test matrix consists of axial HEPA filters rated at 1,000 and 1,500 CFM with varying filter pack types. Each filter is subjected to varying pressures of 14.70 psi, 13.55 psi, 12.40 psi, and 11.25 psi, temperatures of 70 deg. F, 110 deg. F, and 150 deg. F, and relative humidities of 50% and 80%, resulting in twenty-four unique testing points. The corresponding SCFM flow rates are computed for each point in order to form a comparison. Evaluation of the effects of designating ACFM or SCFM will be based upon the filter efficiencies of generated Dioctyl Phthalate (DOP) particles at a diameter of 0.3 micrometers. The results show filter efficiency decreasing as atmospheric pressure decreases due to increased filter face velocity. A correlation model is developed to convert filter efficiency and pressure drop at elevated conditions to those of standard pressure, relative humidity, and temperature.

42 ENGINEERING↗
DOE OSTI2024

Low and moderate x gluon contribution to exclusive Compton scattering processes

We revisit the high energy semi-classical description of the exclusive processes DVCS, TCS, and Double DVCS by explicitly keeping track of the Feynman x dependence in both the hard and the hadronic matrix elements. This is achieved by a modification of the standard shock wave approximation to derive the effective Feynman rules, which leads to a generic expression on which we then perform a partial twist expansion to get rid of quantities suppressed by the proper physical scales. We obtain a compact factorized master formula that can be used to investigate the Bjorken limit at leading twist. In particular, we recover the full one-loop result in the collinear limit for pure gluon exchange with the target. Finally, we discuss the subtleties in taking the simultaneous collinear and small x limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Complex Electrical Conductivity of a Single‐Fractured Rock: Fracture‐and‐Matrix Coupling Mechanism and Aperture Size Predictions

Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2026

Chemical and Radiological Compatibility Testing of 3D Printed Materials

To dramatically increase the adaptability, performance, and safety of processes in support of the Defense Waste Processing Facility (DWPF), Savannah River National Lab (SRNL) plans to perform chemical and radiological compatibility testing on a wide variety of 3D printed materials of interest. The 3D printing process provides numerous strategic operational benefits such as rapid prototyping of complex designs and geometry specific to the needs of the nuclear waste disposition process, as well as on-demand rapid prototyping and iteration with materials that aren’t as accessible through traditional manufacturing methods. Reaction chemistry in simulated waste batches can be matched closely to its radioactive counterpart, but glass reactor vessels have limitations. Vessel geometry can play a big factor in mixing transport limitations, process chemistry, and degradation reaction kinetics. In addition, additive manufacturing allows for much more detailed vessel design than traditional alternatives. Waste processing techniques in DWPF also encounter extreme chemical environments including high pH, strong acids, abrasive slurries, and significant irradiation. To meet these challenges, a matrix of various polymer, ceramic, and metal additive manufacturing materials have been exposed to a suite of chemical environments of interest as well as radioactive dose (such as gamma radiation from 60 Co) to properly test their durability under these conditions. Mass change has been monitored over a period of up to a week in these conditions, as well as added characterization for surface modification through Scanning Electron Microscopy/Electron Dispersive X-ray analysis (SEM/EDX). Further chemical characterization has been monitored through Fourier-Transform InfraRed Spectroscopy (FTIR), with planned investigation via thermal and tensile strength degradation. While the direct product of this research is identification of material(s) that can withstand specific hazardous environments encountered by the mercury water wash tank in DWPF process simulation experiments, the reference base of materials will be used for many other nuclear processes in the pursuit of rapidly developed, cost-efficient, and highly specific devices for environmental remediation and much more.

Wilson, Nathan W. [Savannah River National Laborat↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A review of fluids under nanoconfinement: Reactivity, geomechanics, phase transitions, and flow

Due to its many industrial applications, the physics associated with the nanoconfinement of fluids has garnered great attention from the scientific community in roughly the last fifteen years. One pioneer of that research has been the Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials (MUSE) Energy Frontier Research Center, which produced more than 50 journal articles from 2018 to 2023. We structure the findings of MUSE into four categories and contextualize the work with the broader literature. Here, the four categories are (1) the dependence of chemical reactions on the degree of confinement, which given the lower freedom of motion of molecules yet improved availability of reaction sites generally lead to enhancement from bulk reactions; (2) the mechanical properties of the porous matrix, which are dependent on the heterogeneity of the system in terms of both composition and structure; (3) the phase transitions under nanoconfinement, which are influenced by the pore size and morphology, surface chemistry, and substrate composition; and (4) the physical controls of slip flow, which demonstrate the importance of several factors, such as wettability. These aspects of nanoconfinement are becoming clearer, but there are other elements such as mathematical theory that are still being developed.

58 GEOSCIENCES↗
DOE OSTI2022

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Bioproduction, bioprotection, and biocontainment in multi-kingdom microbial systems with 3D spatial control

Engineered living materials (ELMs) are a class of hybrid materials that include engineered microbes encapsulated by a polymer matrix. The biotic and abiotic components define the ELMs design space and can be altered to improve performance and function. While current synthetic materials in the field display robust biocompatibility with both native and engineered living systems, we have a limited understanding of how to leverage three-dimensional (3D) form factors to spatially organize and control microbial dynamics within the material. Motivated by this knowledge gap, we employed extrusion-based 3D printing to fabricate multi-kingdom hydrogel constructs for the encapsulation of both single and multi-kingdom microbial systems. Core–shell cubic constructs enabled the spatial organization of a constitutive multi-kingdom system of levodopa (L-DOPA)-producing E. coli and betaxanthins (BXN)-producing S. cerevisiae. This spatial organization in 3D materials can introduce precise control over bioproduction, bioprotection, and biocontainment features that are critical to the efficacy of current ELMs. The relative spatial organization of the organisms, as well as the surface area-to-volume ratio were investigated to determine how these design elements impact microbial behavior (metabolite production, growth, expression, and cell distribution) over time. We demonstrated that F127-bis-urethane methacrylate (F127-BUM) core–shell geometries enable the hierarchical 3D printing of multi-kingdom constructs, offering customizable control over bioproduction, bioprotection, and biocontainment. With the optimization of these core–shell structures for continuous bioproduction, these ELMs could be deployed as compact and sustainable bioreactors in remote environments.

additive manufacturing↗
DOE OSTI2023

New insights into the natural history of bronchopulmonary dysplasia from proteomics and multiplexed immunohistochemistry

Bronchopulmonary dysplasia (BPD) is a disease of prematurity related to the arrest of normal lung development. The objective of this study was to better understand how proteome modulation and cell-type shifts are noted in BPD pathology. Pediatric human donors aged 1–3 yr were classified based on history of prematurity and histopathology consistent with “healed” BPD (hBPD, n = 3) and “established” BPD (eBPD, n = 3) compared with respective full-term born (n = 6) age-matched term controls. Proteins were quantified by tandem mass spectroscopy with selected Western blot validations. Multiplexed immunofluorescence (MxIF) microscopy was performed on lung sections to enumerate cell types. Protein abundances and MxIF cell frequencies were compared among groups using ANOVA. Cell type and ontology enrichment were performed using an in-house tool and/or EnrichR. Proteomics detected 5,746 unique proteins, 186 upregulated and 534 downregulated, in eBPD versus control with fewer proteins differentially abundant in hBPD as compared with age-matched term controls. Cell-type enrichment suggested a loss of alveolar type I, alveolar type II, endothelial/capillary, and lymphatics, and an increase in smooth muscle and fibroblasts consistent with MxIF. Histochemistry and Western analysis also supported predictions of upregulated ferroptosis in eBPD versus control. Finally, several extracellular matrix components mapping to angiogenesis signaling pathways were altered in eBPD. Despite clear parsing by protein abundance, comparative MxIF analysis confirms phenotypic variability in BPD. This work provides the first demonstration of tandem mass spectrometry and multiplexed molecular analysis of human lung tissue for critical elucidation of BPD trajectory-defining factors into early childhood.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include a factor of 10^4 improvement in dependency resolution speed, sorting of user objects, ability to compute residuals and Jacobians together, transfer fixes, support for the mortar method in finite volume discretizations, addition of generalized advection schemes for fluid simulations, 10^2 speedup in some Griffin simulations due to a new matrix-only solve type, and much more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fully Homomorphic Encryption

This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.

Lois, Roberts [Idaho National Laboratory (INL), Id↗
DOE OSTI2022

Preexisting Type 1 Diabetes Mellitus Blunts the Development of Posttraumatic Osteoarthritis

Type 1 diabetes mellitus (T1DM) affects 9.5% of the population. T1DM is characterized by severe insulin deficiency that causes hyperglycemia and leads to several systemic effects. T1DM has been suggested as a risk factor for articular cartilage damage and loss, which could expedite the development of osteoarthritis (OA). OA represents a major public health challenge by affecting 300 million people globally, yet very little is known about the correlation between T1DM and OA. In addition, current studies that have looked at the interaction between diabetes mellitus and OA have reported conflicting results with some suggesting a positive correlation whereas others did not. In this study, we aimed to evaluate whether T1DM exacerbates the development of spontaneous OA or accelerates the progression of posttraumatic osteoarthritis (PTOA) after joint injury. Histological evaluation of T1DM and control joints determined that T1DM mice displayed cartilage degeneration measurements consistent with mild OA phenotypes. RNA sequencing analyses identified significantly upregulated genes in T1DM corresponding to matrix-degrading enzymes known to promote cartilage matrix degradation, suggesting a role of these enzymes in OA development. Next, we assessed whether preexisting T1DM influences PTOA development subsequent to trauma. At 6 weeks post-injury, T1DM injured joints displayed significantly less cartilage damage and joint degeneration than injured non-diabetic joints, suggesting a significant delay in PTOA disease progression. At the single-cell resolution, we identified increased number of cells expressing the chondrocyte markers Col2a1, Acan, and Cytl1 in the T1DM injured group. Our findings demonstrate that T1DM can be a risk factor for OA but not for PTOA. This study provides the first account of single-cell resolution related to T1DM and the risk for OA and PTOA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Interplay of surface energy and rheology in biopolymer soil enhancement

Biopolymers such as xanthan gum (XG) and locust bean gum (LBG) hold great potential as eco-friendly alternative soil binders. In this work, we investigated the impact of XG/LBG mixtures on the unconfined compressive strength (UCS) of sand. The high strength of dry biopolymer/sand arises from the cohesion between solid polymer films and sand particles which supported by work of adhesion calculation and soil mechanics measurement. LBG exhibits much lower sand reinforcement efficacy because polymers unevenly distributed within sand matrix. The formation of a core-shell structure in LBG/sand is an interplay of surface free energy and viscoelastic properties of polymer solutions. This structure is altered when LBG mixed with XG at varying ratios as those physical properties changed due to the complexity of polymer chains association. By probing these factors, we aim to elucidate the role of surface energies and polymer physics in governing the strength of the sand/polymer network, thereby contributing to a more comprehensive understanding polymer-sand interface. The low strength of gels (G’ ∼10Pa) cannot solely account for the increased UCS of wet sand over 10 kPa. Instead, the high strength of biopolymer/sand is more likely derived from the granular particles with biopolymers as solid glue.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fiber Sorbents – A Versatile Platform for Sorption-Based Gas Separations

Increasing demand for high-purity fine chemicals and a drive for process intensification of large-scale separations have driven significant work on the development of highly engineered porous materials with promise for sorption-based separations. While sorptive separations in porous materials offer energy-efficient alternatives to longstanding thermal-based methods, the particulate nature of many of these sorbents has sometimes limited their large-scale deployment in high-throughput applications such as gas separations, for which the necessary high feed flow rates and gas velocities accrue prohibitive operational costs. These processability limitations have been historically addressed through powder shaping methods aimed at the fabrication of structured sorbent contactors based on pellets, beads or monoliths, commonly obtained as extrudates. These structures overcome limitations such as elevated pressure drops commonly recorded across powder adsorption beds but often accrue thermal limitations arising from elevated particle density and aggregation, which ultimately cap their maximum separation performance. Furthermore, the harsh mechanical strain to which powder particles are subjected during contactor fabrication, in the form of extrusion/compression forces, can result in partial pore occlusion and framework degradation, further limiting their performance. Here, we present the development of porous fiber sorbents as an alternative sorbent contactor design capable of addressing sorbent processability limitations while enabling an array of performance-maximizing heat integration capabilities. This new sorbent form factor leverages pre-existing know-how from hollow fiber spinning to produce fiber-shaped sorbent contactors through the phase inversion of known polymers in a process known as dry-jet/wet quenching. The process of phase inversion allows microporous sorbent particles to be latched onto a macroporous polymer matrix under mild processing conditions, thus making it compatible with soft porous materials prone to amorphization under traditional pelletization conditions. Sorbent fibers can be created with different geometries through control of the spinning apparatus and process, offering the possibility to produce monolithic and hollow fibers alike, the latter of which can be integrated with thermalization fluid flows. In this Account, we summarize our progress in the field of fiber sorbents from both design and application standpoints. We further guide the reader through the evolution of this field from the early inceptive work on zeolite hollow fibers to recent developments on MOF fibers. We highlight the versatile nature of fiber sorbents, both from the composition, fabrication and structure points of view, and further demonstrate how fiber sorbents offer alternative paths in tackling new and challenging chemical separation challenges like direct air capture (DAC), with a final perspective on the future of the field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Grain boundary facilitated dissolution of nanocrystalline NpO 2 (s) from legacy waste processing

Dissolution of actinide dioxides, including neptunium dioxide (NpO 2 (s)), is paramount for the prediction of the environmental fate of nuclear materials. Quantifying dissolution rates, as well as understanding qualitative dissolution mechanisms, informs performance assessment for geologic disposal of spent nuclear fuel and management of legacy radioactive waste. The aim of this research was to measure the dissolution rate of nanocrystalline NpO 2 (s), produced through legacy nuclear waste processing, under oxidizing conditions, as well as to characterize surface alteration to the material. The solid phase was characterized using electron microscopy techniques (SEM/STEM) and X-ray photoelectron spectroscopy (XPS), indicating preferential dissolution of Np-hydroxide contained in the grain boundaries of NpO 2 (s) and fragmentation of grains from the matrix. The oxidative dissolution was monitored over 40 weeks, yielding a two-step kinetic dissolution model involving hydration of NpO 2 (s) and subsequent oxidation and dissolution of the hydroxide phase. Here, the proposed dissolution models for nanocrystalline NpO 2 (s) suggest that microstructural features such as grain boundaries are key factors affecting dissolution, including release of colloidal particles, and ultimately, environmental fate and transport of nuclear materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A multisite decomposition of the tensor network path integrals

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗