Source author record

Thomas Schuster

Thomas Schuster appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Tensor field tomography with attenuation and refraction: adjoint operators for the dynamic case and numerical experiments

This article is concerned with tensor field tomography in a fairly general setting, that takes refraction, attenuation and time-dependence of tensor fields into account. The mathematical model is given by attenuated ray transforms of the fields along geodesic curves corresponding to a Riemannian metric that is defined by the index of refraction. The data are given at the boundary tangent bundle of the domain and it is well-known that they can be characterized as boundary data of a transport equation turning tensor field tomography into an inverse source problem. This way the adjoint of the forward mapping can be computed using the integral representation or, equivalently, associated to a dual transport equation. The article offers and proves two different representations for the adjoint mappings both in the dynamic and static case. The numerical implementation is demonstrated and evaluated for static fields using the damped Landweber method with Nesterov acceleration applied to both, the integral and PDE-based formulations. The transport equations are solved using a viscosity approximation. The error analysis reveals that the integral representation significantly outperforms PDE-based methods in terms of computational efficiency while achieving comparable reconstruction accuracy. The impact of noise and deviations from straight-line trajectories are investigated confirming improved accuracy if refraction is taken into account. We conclude that the inclusion of refraction to the forward model pays in spite of increased numerical cost.

preprint2022arXiv

A note on $Γ$-convergence of Tikhonov functionals for nonlinear inverse problems

We consider variational regularization of nonlinear inverse problems in Banach spaces using Tikhonov functionals. This article addresses the problem of $Γ$-convergence of a family of Tikhonov functionals and assertions of the convergence of their respective infima. Such questions arise, if model uncertainties, inaccurate forward operators, finite dimensional approximations of the forward solutions and / or data, etc. make the evaluation of the original functional impossible and, thus, its minimizer not computable. But for applications it is of utmost importance that the minimizer of the replacement functional approximates the original minimizer. Under certain additional conditions this is satisfied if the approximated functionals converge to the original functional in the sense of $Γ$-convergence. We deduce simple criteria in different topologies which guarantee $Γ$-convergence as well as convergence of minimizing sequences.

preprint2022arXiv

Information-theoretic Hardness of Out-of-time-order Correlators

We establish that there are properties of quantum many-body dynamics which are efficiently learnable if we are given access to out-of-time-order correlators (OTOCs), but which require exponentially many operations in the system size if we can only measure time-ordered correlators. This implies that any experimental protocol which reconstructs OTOCs solely from time-ordered correlators must be, in certain cases, exponentially inefficient. Our proofs leverage and generalize recent techniques in quantum learning theory. Along the way, we elucidate a general definition of time-ordered versus out-of-time-order experimental measurement protocols, which can be considered as classes of adaptive quantum learning algorithms. Moreover, our results provide a theoretical foundation for novel applications of OTOCs in quantum simulations.

preprint2022arXiv

Learning quantum systems via out-of-time-order correlators

Learning the properties of dynamical quantum systems underlies applications ranging from nuclear magnetic resonance spectroscopy to quantum device characterization. A central challenge in this pursuit is the learning of strongly-interacting systems, where conventional observables decay quickly in time and space, limiting the information that can be learned from their measurement. In this work, we introduce a new class of observables into the context of quantum learning -- the out-of-time-order correlator -- which we show can substantially improve the learnability of strongly-interacting systems by virtue of displaying informative physics at large times and distances. We identify two general scenarios in which out-of-time-order correlators provide a significant advantage for learning tasks in locally-interacting systems: (i) when experimental access to the system is spatially-restricted, for example via a single "probe" degree of freedom, and (ii) when one desires to characterize weak interactions whose strength is much less than the typical interaction strength. We numerically characterize these advantages across a variety of learning problems, and find that they are robust to both read-out error and decoherence. Finally, we introduce a binary classification task that can be accomplished in constant time with out-of-time-order measurements. In a companion paper, we prove that this task is exponentially hard with any adaptive learning protocol that only involves time-ordered operations.

preprint2022arXiv

Many-body quantum teleportation via operator spreading in the traversable wormhole protocol

By leveraging shared entanglement between a pair of qubits, one can teleport a quantum state from one particle to another. Recent advances have uncovered an intrinsically many-body generalization of quantum teleportation, with an elegant and surprising connection to gravity. In particular, the teleportation of quantum information relies on many-body dynamics, which originate from strongly-interacting systems that are holographically dual to gravity; from the gravitational perspective, such quantum teleportation can be understood as the transmission of information through a traversable wormhole. Here, we propose and analyze a new mechanism for many-body quantum teleportation -- dubbed peaked-size teleportation. Intriguingly, peaked-size teleportation utilizes precisely the same type of quantum circuit as traversable wormhole teleportation, yet has a completely distinct microscopic origin: it relies upon the spreading of local operators under generic thermalizing dynamics and not gravitational physics. We demonstrate the ubiquity of peaked-size teleportation, both analytically and numerically, across a diverse landscape of physical systems, including random unitary circuits, the Sachdev-Ye-Kitaev model (at high temperatures), one-dimensional spin chains and a bulk theory of gravity with stringy corrections. Our results pave the way towards using many-body quantum teleportation as a powerful experimental tool for: (i) characterizing the size distributions of operators in strongly-correlated systems and (ii) distinguishing between generic and intrinsically gravitational scrambling dynamics. To this end, we provide a detailed experimental blueprint for realizing many-body quantum teleportation in both trapped ions and Rydberg atom arrays; effects of decoherence and experimental imperfections are analyzed.

preprint2022arXiv

Operator Growth in Open Quantum Systems

The spreading of quantum information in closed systems, often termed scrambling, is a hallmark of many-body quantum dynamics. In open systems, scrambling competes with noise, errors and decoherence. Here, we provide a universal framework that describes the scrambling of quantum information in open systems: we predict that the effect of open-system dynamics is fundamentally controlled by operator size distributions and independent of the microscopic error mechanism. This framework allows us to demonstrate that open quantum systems exhibit universal classes of information dynamics that fundamentally differ from their unitary counterparts. Implications for Loschmidt echo experiments and the classical simulability of open quantum dynamics will be discussed.

preprint2021arXiv

A method for determining the parameters in a rheological model for viscoelastic materials by minimizing Tikhonov functionals

Mathematical models describing the behavior of viscoelastic materials are often based on evolution equations that measure the change in stress depending on its material parameters such as stiffness, viscosity or relaxation time. In this article, we introduce a Maxwell-based rheological model, define the associated forward operator and the inverse problem in order to determine the number of Maxwell elements and the material parameters of the underlying viscoelastic material. We perform a relaxation experiment by applying a strain to the material and measure the generated stress. Since the measured data varies with the number of Maxwell elements, the forward operator of the underlying inverse problem depends on parts of the solution. By introducing assumptions on the relaxation times, we propose a clustering algorithm to resolve this problem. We provide the calculations that are necessary for the minimization process and conclude with numerical results by investigating unperturbed as well as noisy data. We present different reconstruction approaches based on minimizing a least squares functional. Furthermore, we look at individual stress components to analyze different displacement rates. Finally, we study reconstructions with shortened data sets to obtain assertions on how long experiments have to be performed to identify conclusive material parameters.

preprint2020arXiv

Speculative Dereferencing of Registers:Reviving Foreshadow

Since 2016, multiple microarchitectural attacks have exploited an effect that is attributed to prefetching. These works observe that certain user-space operations can fetch kernel addresses into the cache. Fetching user-inaccessible data into the cache enables KASLR breaks and assists various Meltdown-type attacks, especially Foreshadow. In this paper, we provide a systematic analysis of the root cause of this prefetching effect. While we confirm the empirical results of previous papers, we show that the attribution to a prefetching mechanism is fundamentally incorrect in all previous papers describing or exploiting this effect. In particular, neither the prefetch instruction nor other user-space instructions actually prefetch kernel addresses into the cache, leading to incorrect conclusions and ineffectiveness of proposed defenses. The effect exploited in all of these papers is, in fact, caused by speculative dereferencing of user-space registers in the kernel. Hence, mitigation techniques such as KAISER do not eliminate this leakage as previously believed. Beyond our thorough analysis of these previous works, we also demonstrate new attacks enabled by understanding the root cause, namely an address-translation attack in more restricted contexts, direct leakage of register values in certain scenarios, and the first end-to-end Foreshadow (L1TF) exploit targeting non-L1 data. The latter is effective even with the recommended Foreshadow mitigations enabled and thus revives the Foreshadow attack. We demonstrate that these dereferencing effects exist even on the most recent Intel CPUs with the latest hardware mitigations, and on CPUs previously believed to be unaffected, i.e., ARM, IBM, and AMD CPUs.

preprint2019arXiv

Floquet Hopf Insulators

We predict the existence of a novel Floquet topological insulator in three-dimensional two-band systems, the Floquet Hopf insulator, which possesses two distinct topological invariants. One is the Hopf $\mathbb{Z}$ invariant, a linking number characterizing the (non-driven) Hopf topological insulator. The second invariant is an intrinsically Floquet $\mathbb{Z}_2$ invariant, and represents a condensed matter realization of the topology underlying the Witten anomaly in particle physics. Both invariants arise from topological defects in the system's time-evolution, subject to a process in which defects at different quasienergy exchange even amounts of topological charge. Their contrasting classifications lead to a measurable physical consequence, namely, an unusual bulk-boundary correspondence where gapless edge modes are topologically protected, but may exist at either 0- or $π$-quasienergy. Our results represent a phase of matter beyond the conventional classification of Floquet topological insulators.

preprint2016arXiv

A CG-type method in Banach spaces with an application to computerized tomography

Conjugate Gradient (CG) methods are one of the most effective iterative methods to solve linear equations in Hilbert spaces. So far, they have been inherently bound to these spaces since they make use of the inner product structure. In more general Banach spaces one of the most prominent iterative solvers are Landweber-type methods that essentially resemble the Steepest Descent method applied to the normal equation. More advanced are subspace methods that take up the idea of a Krylov-type search space, wherein an optimal solution is sought. However, they do not share the conjugacy property with CG methods. In this article we propose that the Sequential Subspace Optimization (SESOP) method can be considered as an extension of CG methods to Banach spaces. We employ metric projections to orthogonalize the current search direction with respect to the search space from the last iteration. For the l2-space our method then exactly coincides with the Polak-Ribière type of the CG method when applied to the normal equation. We show that such an orthogonalized search space still leads to weak convergence of the subspace method. Moreover, numerical experiments on a random matrix toy problem and 2D computerized tomography on lp-spaces show superior convergence properties over all p compared to non-orthogonalized search spaces. This especially holds for lp-spaces with small p. We see that the closer we are to an l2-space, the more we recover of the conjugacy property that holds in these spaces, i. e., as expected, the more the convergence behaves independently of the size of the truncated search space.

preprint2016arXiv

A modified algebraic reconstruction technique taking refraction into account with an application in terahertz tomography

Terahertz (THz) tomography is a rather novel technique for nondestructive testing that is particularly suited for the testing of plastics and ceramics. Previous publications showed a large variety of conventional algorithms adapted from computed tomography or ultrasound tomography which were directly applied to THz tomography. Conventional algorithms neglect the specific nature of THz radiation, i.e. refraction at interfaces, reflection losses and the beam profile (Gaussian beam), which results in poor reconstructions. The aim is the efficient reconstruction of the complex refractive index, since it indicates inhomogeneities in the material. A hybrid algorithm has been developed based on the algebraic reconstruction technique (ART). ART is adapted by including refraction (Snell's law) and reflection losses (Fresnel equations). Our method uses a priori information about the interface and layer geometry of the sample. This results in the 'Modified ART for THz tomography', which reconstructs simultaneously the complex refractive index from transmission coefficient and travel time measurements.

preprint2016arXiv

An improved exact inversion formula for cone beam vector tomography

In this article we present an improved exact inversion formula for the 3D cone beam transform of vector fields. It is well known that only the solenoidal part of a vector field can be determined by the longitudinal ray transform of a vector field in cone beam geometry. The exact inversion formula, as it was developed in A. Katsevich and T. Schuster, 'An exact inversion formula for cone beam vector tomography', Inverse Problems 29 (2013), consists of two parts. the first part is of filtered backprojection type, whereas the second part is a costly 4D integration and very inefficient. In this article we tackle this second term and achieve an improvement which is easily to implement and saves one order of integration. The theory says that the first part contains all information about the curl of the field, whereas the second part presumably has information about the boundary values. This suggestion is supported by the fact that the second part vanishes if the exact field is divergence free and tangential at the boundary. A number of numerical tests, that are also subject of this article, confirm the theoretical results and the exactness of the formula.

preprint2016arXiv

Dissipationless conductance in a topological coaxial cable

We present a dynamical mechanism leading to dissipationless conductance, whose quantized value is controllable, in a (3+1)-dimensional electronic system. The mechanism is exemplified by a theory of Weyl fermions coupled to a Higgs field--also known as an axion insulator. We show that the insertion of an axial gauge flux can induce vortex lines in the Higgs field, similarly to the development of vortices in a superconductor upon the insertion of magnetic flux. We further show that the necessary axial gauge flux can be generated using Rashba spin-orbit coupling or a magnetic field. Vortex lines in the Higgs field are known to bind chiral fermionic modes, each of which serves as a one-way channel for electric charge with conductance $e^2/h$. Combining these elements, we present a physical picture, the "topological coaxial cable," illustrating how the value of the quantized conductance could be controlled in such an axion insulator.

preprint2016arXiv

Non-Abelian braiding of light

Many topological phenomena first proposed and observed in the context of electrons in solids have recently found counterparts in photonic and acoustic systems. In this work, we demonstrate that non-Abelian Berry phases can arise when coherent states of light are injected into "topological guided modes" in specially-fabricated photonic waveguide arrays. These modes are photonic analogues of topological zero modes in electronic systems. Light traveling inside spatially well-separated topological guided modes can be braided, leading to the accumulation of non-Abelian phases, which depend on the order in which the guided beams are wound around one another. Notably, these effects survive the limit of large photon occupation, and can thus also be understood as wave phenomena arising directly from Maxwell's equations, without resorting to the quantization of light. We propose an optical interference experiment as a direct probe of this non-Abelian braiding of light.

preprint2016arXiv

Sequential subspace optimization for nonlinear inverse problems

In this work we discuss a method to adapt sequential subspace optimization (SESOP), which has so far been developed for linear inverse problems in Hilbert and Banach spaces, to the case of nonlinear inverse problems. We start by revising the well-known technique for Hilbert spaces. In a next step, we introduce a method using multiple search directions that are especially designed to fit the nonlinearity of the forward operator. To this end, we iteratively project the initial value onto stripes whose shape is determined by the search direction, the nonlinearity of the operator and the noise level. We additionally propose a fast algorithm that uses two search directions. Finally we will show convergence and regularization properties for the presented method.

preprint2015arXiv

An iterative method to reconstruct the refractive index of a medium from time-of-flight measurements

The article deals with a classical inverse problem: the computation of the refractive index of a medium from ultrasound time-of-flight (TOF) measurements. This problem is very popular in seismics but also for tomographic problems in inhomogeneous media. For example ultrasound vector field tomography needs a priori knowledge of the sound speed. According to Fermat's principle ultrasound signals travel along geodesic curves of a Riemannian metric which is associated with the refractive index. The inverse problem thus consists of determining the index of refraction from integrals along geodesics curves associated with the integrand leading to a nonlinear problem. In this article we describe a numerical solver for this problem scheme based on an iterative minimization method for an appropriate Tikhonov functional. The outcome of the method is a stable approximation of the sought index of refraction as well as a corresponding set of geodesic curves. We prove some analytical convergence results for this method and demonstrate its performance by means of several numerical experiments. Another novelty in this article is the explicit representation of the backprojection operator for the ray transform in Riemannian geometry and its numerical realization relying on a corresponding phase function that is determined by the metric. This gives a natural extension of the conventional backprojection from 2D computerized tomography to inhomogeneous geometries.

preprint2015arXiv

On the linearization of identifying the stored energy function of a hyperelastic material from full knowledge of the displacement field

We consider the nonlinear,inverse problem of computing the stored energy function of a hyperelastic material from the full knowledge of the displacement field. The displacement field is described as solution of the nonlinear, dynamic, elastic wave equation, where the first Piola-Kirchhoff stress tensor is given as the gradient of the stored energy function. We assume that we have a dictionary at hand such that the energy function is given as a conic combination of the dictionary's elements. In that sense the mathematical model of the direct problem is the nonlinear operator that maps the vector of expansion coefficients to the solution of the hyperelastic wave equation. In this article we summarize some continuity results for this operator and deduce its Fréchet derivative as well as the adjoint of this derivative. Since the stored energy function encodes mechanical properties of the underlying, hyperelastic material, the considered inverse problem is of highest interest for structural health monitoring systems where defects are detected from boundary measurements of the displacement field. For solving the inverse problem iteratively by the Landweber method or Newton type methods, the knowledge of the Fréchet derivative and its adjoint is of utmost importance.

preprint2014arXiv

Uniqueness and stability result for Cauchy's equation of motion for a certain class of hyperelastic materials

We consider Cauchy's equation of motion for hyperelastic materials. The solution of this nonlinear initial-boundary value problem is the vector field which discribes the displacement which a particle of this material perceives when exposed to stress and external forces. This equation is of greatest relevance when investigating the behaviour of elastic, anisotropic composites and for the detection of defects in such materials from boundary measurements. This is why results on unique solvability and continuous dependence from the initial values are of large interest in materials research and structural health monitoring. In this article we present such a result, provided that reasonable smoothness assumptions for the displacement field and the boundary of the domain are satisfied, for a certain class of hyperelastic materials where the first Piola-Kirchhoff tensor is written as a conic combination of finitely many given tensors.