Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor networks”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermalization at low temperatures via weakly damped multisite baths

We study the thermalization properties of one-dimensional open quantum systems coupled to baths at their boundary. The baths are driven to their thermal states via Lindblad operators, while the system undergoes Hamiltonian dynamics. We specifically consider multisite baths and investigate the extent to which the late-time steady state resembles a Gibbs state at some controllable temperature set by the baths. We study three models: a noninteracting fermion model accessible via free-fermion technology, and two interacting models, the XZ model and the chiral clock model, which are accessible via tensor network methods. We show that, by tuning towards the weak coupling and slow relaxation limits, one can engineer low temperatures in the bulk of the system provided the bath size is big enough. Here, we use this capability to study energy transport in the XZ model at lower temperatures than previously reported. Our work paves the way for future studies of interacting open quantum systems at low temperatures.

1-dimensional spin chains↗
DOE OSTI2023

Quantum criticality and topology in non-equilibrium systems

A key goal of condensed matter physics research is to identify new phases of matter, and to understand the universal features of the phase transitions between them. In the past decade, physicists have uncovered a wealth of new phases with interesting surface properties, exemplified by the theoretical prediction and subsequent experimental discovery of topological insulators and superconductors. Traditional condensed matter systems are usually in a thermal equilibrium state and typically at very low temperature. Very recently, experimental advances have sparked interest in the non-equilibrium setting. Non-equilibrium systems can host new phases and phenomena with no equilibrium counterpart, and could also enable robust ways to build quantum memory devices to store and manipulate quantum information in a coherent manner. These phases and phenomena are inherently “dynamical”: they are described not by changes in the arrangement or structure of the constituent particles, but instead marked by sharp distinctions in how the particles move and exchange energy or quantum information. The discovery of robust non-equilibrium phases raises many fundamental questions: Can we develop a systematic theory of states of matter and of dynamical transitions between such states? How can such states be realized and probed experimentally? The main goal of this project was to explore the emergence of topological phases and quantum criticality (two cornerstones of modern condensed matter physics in equilibrium) in such non-equilibrium quantum systems. Specific goals included (1) using tensor networks to efficiently represent non-equilibrium states of matter and their phase transitions; (2) studying and designing new probes for periodically driven systems; and (3) developing analytic and numerical tools to analyze non-equilibrium topological phase transitions. Advances in these directions were achieved using novel techniques appropriate to study the non-equilibrium dynamics in many-body quantum systems combining strong interactions and randomness. Taken together, these results provide a new conceptual framework for understanding the emergence of quantum critical and topological properties in quantum systems far from thermal equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Learning to classify quantum phases of matter with a few measurements

We study the identification of quantum phases of matter, at zero temperature, when only part of the phase diagram is known in advance. Following a supervised learning approach, we show how to use our previous knowledge to construct an observable capable of classifying the phase even in the unknown region. By using a combination of classical and quantum techniques, such as tensor networks, kernel methods, generalization bounds, quantum algorithms, and shadow estimators, we show that, in some cases, the certification of new ground states can be obtained with a polynomial number of measurements. An important application of our findings is the classification of the phases of matter obtained in quantum simulators, e.g. cold atom experiments, capable of efficiently preparing ground states of complex many-particle systems and applying simple measurements, e.g. single qubit measurements, but unable to perform a universal set of gates.

quantum machine learning↗
DOE OSTI2025

Almost strong zero modes at finite temperature

Interacting fermionic chains exhibit extended regions of topological degeneracy of their ground states as a result of the presence of Majorana or parafermionic zero modes localized at the edges. In the opposite limit of infinite temperature, the corresponding nonintegrable spin chains, obtained via generalized Jordan-Wigner mapping, are known to host so-called almost strong zero modes, which are long-lived with respect to any bulk excitations. Here we study the fairly unexplored territory that bridges these two extreme cases of zero and infinite temperature. We blend two established techniques for states, the Lanczos series expansion and a tensor network ansatz, uplifting them to the level of operator algebra. This allows us to efficiently simulate large system sizes for arbitrarily long timescales and to extract the temperature-dependent decay rates. We observe that for the Kitaev-Hubbard model, the decay rate of the edge mode depends exponentially on the inverse temperature 𝛽, and on an effective energy scale Δ eff that is greater than the thermodynamic gap of the system Δ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2025

Thermal Radiation Transport with Tensor Trains

We present a novel tensor network algorithm to solve the time-dependent, gray thermal radiation transport equation. The method invokes a tensor train (TT) decomposition for the specific intensity. The efficiency of this approach is dictated by the rank of the decomposition. When the solution is “low rank,” the memory footprint of the specific intensity solution vector may be significantly compressed. The algorithm, following a step-then-truncate approach of a traditional discrete ordinates method, operates directly on the compressed state vector, thereby enabling large speedups for low-rank solutions. To achieve these speedups, we rely on a recently developed rounding approach based on the Gram-SVD. We detail how familiar S N algorithms for (gray) thermal transport can be mapped to this TT framework and present several numerical examples testing both the optically thick and thin regimes. The TT framework finds low-rank structure and supplies up to ≃60× speedups and ≃1000× compressions for problems demanding large angle counts, thereby enabling previously intractable SN calculations and supplying a promising avenue to mitigate ray effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Hollow-grams: generalized entanglement wedges from the gravitational path integral

Recently, Bousso and Penington (BP) made a proposal for the entanglement wedge associated to a gravitating bulk region. In this paper, we derive this proposal in time-reflection symmetric settings using the gravitational path integral. To do this, we exploit the connection between random tensor networks (RTNs) and fixed-geometry states in gravity. We define the entropy of a bulk region in an RTN by removing tensors in that region and computing the entropy of the open legs thus generated in the “hollowed” RTN. We thus derive the BP proposal for RTNs and hence, also for fixed-geometry states in gravity. By then expressing a general holographic state as a superposition over fixed-geometry states and using a diagonal approximation, we provide a general gravitational path integral derivation of the BP proposal. We demonstrate that the saddles computing the Rényi entropy Sn depend on how the bulk region is gauge-invariantly specified. Nevertheless, we show that the BP proposal is universally reproduced in the n → 1 limit.

2D Gravity↗
DOE OSTI2025

The design space of E(3)-equivariant atom-centred interatomic potentials

Abstract Molecular dynamics simulation is an important tool in computational materials science and chemistry, and in the past decade it has been revolutionized by machine learning. This rapid progress in machine learning interatomic potentials has produced a number of new architectures in just the past few years. Particularly notable among these are the atomic cluster expansion, which unified many of the earlier ideas around atom-density-based descriptors, and Neural Equivariant Interatomic Potentials (NequIP), a message-passing neural network with equivariant features that exhibited state-of-the-art accuracy at the time. Here we construct a mathematical framework that unifies these models: atomic cluster expansion is extended and recast as one layer of a multi-layer architecture, while the linearized version of NequIP is understood as a particular sparsification of a much larger polynomial model. Our framework also provides a practical tool for systematically probing different choices in this unified design space. An ablation study of NequIP, via a set of experiments looking at in- and out-of-domain accuracy and smooth extrapolation very far from the training data, sheds some light on which design choices are critical to achieving high accuracy. A much-simplified version of NequIP, which we call BOTnet (for body-ordered tensor network), has an interpretable architecture and maintains its accuracy on benchmark datasets.

Computer Science↗
DOE OSTI2024

Twice upon a time: timelike-separated quantum extremal surfaces

The Python’s Lunch conjecture for the complexity of bulk reconstruction involves two types of nonminimal quantum extremal surfaces (QESs): bulges and throats, which differ by their local properties. The conjecture relies on the connection between bulk spatial geometry and quantum codes: a constricting geometry from bulge to throat encodes the bulk state nonisometrically, and so requires an exponentially complex Grover search to decode. However, thus far, the Python’s Lunch conjecture is only defined for spacetimes where all QESs are spacelike-separated from one another. Here we explicitly construct (time-reflection symmetric) spacetimes featuring both timelike-separated bulges and timelike-separated throats. Interestingly, all our examples also feature a third type of QES, locally resembling a de Sitter bifurcation surface, which we name a bounce. By analyzing the Hessian of generalized entropy at a QES, we argue that this classification into throats, bulges and bounces is exhaustive. We then propose an updated Python’s Lunch conjecture that can accommodate general timelike-separated QESs and bounces. Notably, our proposal suggests that the gravitational analogue of a tensor network is not necessarily the time-reflection symmetric slice, even when one exists.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Entanglement negativity and replica symmetry breaking in general holographic states

The entanglement negativity $\mathcal{E}$(A : B) is a useful measure of quantum entanglement in bipartite mixed states. In random tensor networks (RTNs), which are related to fixed-area states, it was found in ref. [1] that the dominant saddles computing the even Rényi negativity $\mathcal{E}$ (2k) generically break the ℤ 2k replica symmetry. This calls into question previous calculations of holographic negativity using 2D CFT techniques that assumed ℤ 2k replica symmetry and proposed that the negativity was related to the entanglement wedge cross section. In this paper, we resolve this issue by showing that in general holographic states, the saddles computing $\mathcal{E}$ (2k) indeed break the ℤ 2k replica symmetry.

AdS-CFT Correspondence↗
DOE OSTI2024

Constrained HRT Surfaces and their Entropic Interpretation

Abstract Consider two boundary subregionsAandBthat lie in a common boundary Cauchy surface, and consider also the associated HRT surfaceγ B forB. In that context, the constrained HRT surfaceγ A:B can be defined as the codimension-2 bulk surface anchored toAthat is obtained by a maximin construction restricted to Cauchy slices containingγ B . As a result,γ A:B is the union of two pieces,$$ {\gamma}_{A:B}^B $$ γ A : B B and$$ {\gamma}_{A:B}^{\overline{B}} $$ γ A : B B ¯ lying respectively in the entanglement wedges ofBand its complement$$ \overline{B} $$ B ¯ . Unlike the area$$ \mathcal{A}\left({\gamma}_A\right) $$ A γ A of the HRT surfaceγ A , at least in the semiclassical limit, the area$$ \mathcal{A}\left({\gamma}_{A:B}\right) $$ A γ A : B ofγ A:B commutes with the area$$ \mathcal{A}\left({\gamma}_B\right) $$ A γ B ofγ B . To study the entropic interpretation of$$ \mathcal{A}\left({\gamma}_{A:B}\right) $$ A γ A : B , we analyze the Rényi entropies of subregionAin a fixed-area state of subregionB. We use the gravitational path integral to show that then ≈1 Rényi entropies are then computed by minimizing$$ \mathcal{A}\left({\gamma}_A\right) $$ A γ A over spacetimes defined by a boost angle conjugate to$$ \mathcal{A}\left({\gamma}_B\right) $$ A γ B . In the case where the pieces$$ {\gamma}_{A:B}^B $$ γ A : B B and$$ {\gamma}_{A:B}^{\overline{B}} $$ γ A : B B ¯ intersect at a constant boost angle, a geometric argument shows that then ≈1 Rényi entropy is then given by$$ \frac{\mathcal{A}\left({\gamma}_{A:B}\right)}{4G} $$ A γ A : B 4 G . We discuss how then ≈1 Rényi entropy differs from the von Neumann entropy due to a lack of commutativity of then→ 1 andG →0 limits. We also discuss how the behaviour changes as a function of the width of the fixed-area state. Our results are relevant to some of the issues associated with attempts to use standard random tensor networks to describe time dependent geometries.

Physics↗
DOE OSTI2025

Towards a real-time computation of timelike hadronic vacuum polarization and light-by-light scattering: Schwinger Model tests

Hadronic vacuum polarization (HVP) and light-by-light scattering (HLBL) are crucial for evaluating the Standard Model predictions concerning the muon’s anomalous magnetic moment. However, direct first-principle lattice gauge theory-based calculations of these observables in the timelike region remain challenging. Discrepancies persist between lattice quantum chromodynamics (QCD) calculations in the spacelike region and dispersive approaches relying on experimental data parametrization from the timelike region. Here, we introduce a methodology employing 1+1-dimensional quantum electrodynamics (QED), i.e. the Schwinger Model, to investigate the HVP and HLBL. To that end, we use both tensor network techniques, specifically matrix product states, and classical emulators of digital quantum computers. Demonstrating feasibility in a simplified model, our approach sets the stage for future endeavors leveraging digital quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗
DOE OSTI2026

Quantum dynamics of cosmological particle production: interacting quantum field theories with matrix product states

Understanding real-time dynamics of interacting quantum fields in curved space-time remains a major theoretical challenge. We employ tensor network methods to study such dynamics using interacting scalar and gauge theories in 1+1 spacetime dimensions, subject to a quench modeling a homogeneously expanding gravitational background. The models considered are the scalar λϕ 4 theory and the Schwinger model, i.e. a Dirac fermion coupled to a U(1) gauge field which is equivalent via bosonization to a scalar field with a cosine self-interaction. In the free scalar limit, both theories reproduce known analytical results, providing a nontrivial numerical validation of bosonization in curved spacetime for the Schwinger model. Our central finding is that self-interactions lead to a suppression of gravitational particle production compared to the free-field case, as evidenced by two-point functions and the spectra of produced particles. We further examine the behavior of entanglement generation and find that interactions suppress entanglement growth in the λϕ 4 theory, while in the Schwinger model, the interplay between suppressed particle production and enhanced inter-particle correlations leads to more complex entanglement behavior. Our results pave the way for further explorations of nonperturbative quantum real-time dynamics of interacting scalar and gauge theories in arbitrary gravitational backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Entanglement Renormalization for Quantum Field Theories with Discrete Wavelet Transforms

We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Simple holography in general spacetimes

The simple or “outermost” wedge in AdS is the portion of the entanglement wedge that can be reconstructed with sub-exponential effort from CFT data. Here we furnish a definition in arbitrary spacetimes: given an input wedge a analogous to a CFT boundary region, the simple wedge z(a) is the largest wedge accessible by a “zigzag,” a certain sequence of antinormal lightsheets. We show that z(a) is a throat, and that it is contained in every other throat. This implies that z(a) is unique; that it is contained in the generalized entanglement wedge; and that it reduces to the AdS prescription as a special case. The zigzag explicitly constructs a preferred Cauchy slice that renders the simple wedge accessible from a; thus it adds a novel structure even in AdS. So far, no spacelike construction is known to reproduce these results, even in time-symmetric settings. This may have implications for the modeling of holographic encoding by tensor networks.

AdS-CFT Correspondence↗
DOE OSTI2020

Islands in cosmology

A quantum extremal island suggests that a region of spacetime is encoded in the quantum state of another system, like the encoding of the black hole interior in Hawking radiation. We study conditions for islands to appear in general spacetimes, with or without black holes. They must violate Bekenstein’s area bound in a precise sense, and the boundary of an island must satisfy several other information-theoretic inequalities. These conditions combine to impose very strong restrictions, which we apply to cosmological models. We find several examples of islands in crunching universes. In particular, in the four-dimensional FRW cosmology with radiation and a negative cosmological constant, there is an island near the turning point when the geometry begins to recollapse. In a two-dimensional model of JT gravity in de Sitter spacetime, there are islands inside crunches that are encoded at future infinity or inside bubbles of Minkowski spacetime. Finally, we discuss simple tensor network toy models for islands in cosmology and black holes.

79 ASTRONOMY AND ASTROPHYSICS↗