On the Development of SCILAB Compatible Software for the Analysis and Control of Repetitive Processes

In this paper further results on the development of a SCILAB compatible software package for the analysis and control of repetitive processes is described. The core of the package consists of a simulation tool which enables the user to inspect the response of a given example to an input, design a control law for stability and/or performance, and also simulate the response of a controlled process to a specified reference signal

On the Control of Distributed Parameter Systems using a Multidimensional Systems Setting

The unique characteristic of a repetitive process is a series of sweeps, termed passes, through a set of dynamics defined over a finite duration with resetting before the start of the each new one. On each pass an output, termed the pass profile is produced which acts as a forcing function on, and hence contributes to, the dynamics of the next pass profile. This leads to the possibility that the output, i.e. the sequence of pass profiles, will contain oscillations which increase in amplitude in the pass-to-pass direction. Such behavior cannot be controlled by standard linear systems approach and instead they must be treated as a multidimensional system, i.e. information propagation in more than one independent direction. Physical examples of such processes include long-wall coal cutting and metal rolling. In this paper, stability analysis and control systems design algorithms are developed for a model where a plane, or rectangle, of information is propagated in the passto- pass direction. The possible use of these in the control of distributed parameter systems is then described using a fourthorder wavefront equation

Thermal evolution of the Schwinger model with Matrix Product Operators

We demonstrate the suitability of tensor network techniques for describing the thermal evolution of lattice gauge theories. As a benchmark case, we have studied the temperature dependence of the chiral condensate in the Schwinger model, using matrix product operators to approximate the thermal equilibrium states for finite system sizes with non-zero lattice spacings. We show how these techniques allow for reliable extrapolations in bond dimension, step width, system size and lattice spacing, and for a systematic estimation and control of all error sources involved in the calculation. The reached values of the lattice spacing are small enough to capture the most challenging region of high temperatures and the final results are consistent with the analytical prediction by Sachs and Wipf over a broad temperature range.Comment: 6 pages, 11 figure

A perturbative gadget for delaying the onset of barren plateaus in variational quantum algorithms

Variational quantum algorithms are being explored as a promising approach to finding useful applications for noisy intermediate-scale quantum computers. However, cost functions corresponding to many problems of interest are inherently global, defined by Hamiltonians with many-body interactions. Consequently, the optimization landscape can exhibit exponentially vanishing gradients, so-called barren plateaus, rendering optimal solutions difficult to find. Strategies for mitigating barren plateaus are therefore needed to make variational quantum algorithms trainable and capable of running on larger-scale quantum computers. In this work, we contribute the toolbox of perturbative gadgets to the portfolio of methods being explored in the quest for making noisy intermediate-scale quantum devices useful. We introduce a novel perturbative gadget, tailored to variational quantum algorithms, that can be used to delay the onset of barren plateaus. Our perturbative gadget encodes an arbitrary many-body Hamiltonian corresponding to a global cost function into the low-energy subspace of a three-body Hamiltonian. Our construction requires $rk$ additional qubits for a $k$-body Hamiltonian comprising $r$ terms. We provide guarantees on the closeness of global minima and prove that the local cost function defined by our three-body Hamiltonian exhibits non-vanishing gradients. We then provide numerical demonstrations to show the functioning of our approach and discuss heuristics that might aid its practical implementation.Comment: Added further discussion of the number of required measurement

Towards overcoming the Monte Carlo sign problem with tensor networks

The study of lattice gauge theories with Monte Carlo simulations is hindered by the infamous sign problem that appears under certain circumstances, in particular at non-zero chemical potential. So far, there is no universal method to overcome this problem. However, recent years brought a new class of non-perturbative Hamiltonian techniques named tensor networks, where the sign problem is absent. In previous work, we have demonstrated that this approach, in particular matrix product states in 1+1 dimensions, can be used to perform precise calculations in a lattice gauge theory, the massless and massive Schwinger model. We have computed the mass spectrum of this theory, its thermal properties and real-time dynamics. In this work, we review these results and we extend our calculations to the case of two flavours and non-zero chemical potential. We are able to reliably reproduce known analytical results for this model, thus demonstrating that tensor networks can tackle the sign problem of a lattice gauge theory at finite densit

Non-perturbative Test of the Witten-Veneziano Formula from Lattice QCD

We compute both sides of the Witten-Veneziano formula using lattice techniques. For the one side we perform dedicated quenched simulations and use the spectral projector method to determine the topological susceptibility in the pure Yang-Mills theory. The other side we determine in lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions including for the first time also the flavour singlet decay constant. The Witten-Veneziano formula represents a leading order expression in the framework of chiral perturbation theory and we also employ leading order chiral perturbation theory to relate the flavor singlet decay constant to the relevant decay constant parameters in the quark flavor basis and flavor non-singlet decay constants. After taking the continuum and the SU$(2)$ chiral limits we compare both sides and find good agreement within uncertainties.Comment: 30 pages, 7 figures, version accepted for publicatio

Phase structure of the (1+1)-dimensional massive Thirring model from matrix product states

Employing matrix product states as an ansatz, we study the non-thermal phase structure of the (1+1)-dimensional massive Thirring model in the sector of vanishing total fermion number with staggered regularization. In this paper, details of the implementation for this project are described. To depict the phase diagram of the model, we examine the entanglement entropy, the fermion bilinear condensate and two types of correlation functions. Our investigation shows the existence of two phases, with one of them being critical and the other gapped. An interesting feature of the phase structure is that the theory with non-zero fermion mass can be conformal. We also find clear numerical evidence that these phases are separated by a transition of the Berezinskii-Kosterlitz-Thouless type. Results presented in this paper establish the possibility of using the matrix product states for probing this type of phase transition in quantum field theories. They can provide information for further exploration of scaling behaviour, and serve as an important ingredient for controlling the continuum extrapolation of the model.Comment: 31 pages, 18 figures; minor changes to the text, typos corrected, references added; version published in Physical Review