DOCSLIB.ORG

Symmetries and Self-Similarity of Many-Body Wavefunctions Arxiv:1412.6796V1 [Quant-Ph] 21 Dec 2014

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

A Thesis submitted for the degree of Doctor of Philosophy Symmetries and self-similarity of many-body wavefunctions Piotr Migdał advisor: Maciej Lewenstein co-advisor: Javier Rodríguez-Laguna submitted: Sept 2014, defended: Dec 2014 ICFO–Institut de Ciències Fotòniques, 08860 Castelldefels, Spain arXiv:1412.6796v1 [quant-ph] 21 Dec 2014 Contents 1 Introduction1 1.1 Background . .1 1.1.1 Historical background . .1 1.1.2 Properties of quantum mechanics . .3 1.1.3 Entanglement . .5 1.1.4 Quantum information with photons . .6 1.1.5 Entanglement invariants for symmetric states . .7 1.1.6 Quantum sequences and qubism . .8 1.1.7 Complexity and the Rényi entropy . .8 1.1.8 Quantum complex networks . 10 1.2 Structure . 12 1.3 Contribution . 13 2 Invariants for bosonic and symmetric states 15 2.1 Introduction . 15 2.2 Overview . 16 2.3 Preliminaries . 18 2.3.1 Representations and notation for symmetric states . 18 2.3.2 Local unitary equivalence . 19 2.3.3 Invariants . 22 2.3.4 Simple examples . 23 2.4 SLOCC and LOCC equivalence of permutationally symmetric pure states 26 2.4.1 Symmetric operations suffice . 28 2.4.2 Explicit formula for symmetrization . 31 2.4.3 Symmetry classes from single-particle stabilizers . 31 2.5 Invariants as functions of creation and annihilation operators . 36 2.5.1 Reduced density matrix . 37 2.5.2 Spectral method . 37 2.5.3 Symmetric component of tensor powers . 40 2.5.4 Examples . 43 2.5.5 Experimental recipe for tensor product in mode basis . 48 2.5.6 Schwinger representation of symmetric operators . 51 2.6 Singlet space for photons and information protection . 55 2.6.1 Example with logical qubits . 56 2.6.2 Experimental scheme . 58 2.6.3 General proof . 59 2.7 Further questions . 62 i Contents 3 Visualization of quantum sequences 65 3.1 Introduction . 65 3.1.1 Classical sequence analysis . 66 3.1.2 Data visualization . 67 3.1.3 Visualizing sequences . 68 3.2 Qubism . 69 3.2.1 Basic mapping . 69 3.2.2 Properties . 72 3.2.3 Examples . 73 3.2.4 General framework . 79 3.2.5 Mappings . 81 3.2.6 Entanglement visualization . 86 3.2.7 Fractal dimension of the state . 89 3.2.8 Qubism for mixed states and operators . 92 3.2.9 Discussion . 94 4 Quantum walks on complex networks 97 4.1 Introduction . 97 4.1.1 Networks and quantum walks . 97 4.1.2 Community detection . 98 4.1.3 Structure . 100 4.2 Classical and quantum walk . 101 4.2.1 Walks framework . 101 4.2.2 Degree distribution and quantumness . 105 4.2.3 Numerical results . 106 4.3 Community detection for quantum walks . 109 4.3.1 Community detection . 109 4.3.2 Quantum community detection . 111 4.3.3 Performance analysis . 115 4.4 Appendix . 120 4.4.1 Definitions . 120 4.4.2 Closeness measures . 123 5 Conclusion 127 ii Acknowledgments Science, even in its naked form, is a collective, social process [1]. For science any theory, observation or discovery is related to people and their relations. For a scientist it is all above plus much more. I am indebted to Maciek1 Lewenstein, especially for providing me a lot of freedom for pursuing my diverse scientific and educational interests, and his belief in me. I am convinced that there is no better gift for independence and creativity that never saying “no”. I am grateful to Javi Rodríguez-Laguna for his countless insights into anything, from scientific remarks on current projects (and unrelated ones), through pieces of advice on academic writing and workflow, to comments on education, society and, well, anything. It was encouraging, inspiring and fruitful. I would like to thank Jake Biamonte for inviting me to Turin for an intensive and fascinating research collaboration. Even tough it was a short stay, it was wonderful on so many axes. I am happy for hospitality and help of the administration: primarily my home institute, ICFO in Castelldefels, and also guest institutes, such as ISI Foundation in Turin and IFT in Madrid. It is thank to your work and attitude that I felt welcomed, and free from paperwork burden. I am grateful to my family and close friends for the support and encouragement. Out of many lessons learnt the most important one is that, in a long run, happiness is as important as intellectual prowess. The list only starts here. All coauthors, discussion partners, lecturers, fellow PhD students and friends — thank you! 1Q: Is it misspelled? A: No. Keys [j] and [k] are neighbors, but this time it is misleading: Maciek is a casual form of Maciej. iii Contents This PhD was supported by Spanish MINCIN/MINECO project TOQATA (FIS2008- 00784), EU Integrated Projects AQUTE and SIQS, CHISTERA project DIQUIP, ERC grants QUAGATUA and OSYRIS. iv Abstract The study of the structure of quantum states can provide insight into the possi- bilities of quantum mechanics applied to quantum communication, cryptography and computations, as well as the study of condensed matter systems. For example, it shows the physical restrictions on the ways how a quantum state can be used and allows us to tell which quantum states are equivalent up to local operations. Therefore, it is crucial for any analysis of the properties and applications of quantum states. This PhD thesis is dedicated to the study of the interplay between symmetries of quantum states and their self-similar properties. It consists of three connected threads of research: polynomial invariants for multiphoton states, visualization schemes for quantum many-body systems and a complex networks approach to quantum walks on a graph. First, we study the problem of which many-photon states are equivalent up to the action of passive linear optics. We prove that it can be converted into the problem of equivalence of two permutation-symmetric states, not necessarily restricted to the same operation on all parties. We show that the problem can be formulated in terms of symmetries of complex polynomials of many variables, and provide two families of invariants, which are straightforward to compute and provide analytical results. Fur- thermore, we prove that some highly symmetric states (singlet states implemented with photons) offer two degrees of robustness — both against collective decoherence and against a photon loss. Additionally, we provide two proposals for experiments, feasible with an optical setup and current technology: one related to the direct mea- surement of a family of invariants using photon-counting, and the other concerting the protection of transmitted quantum information employing the symmetries of the state. Second, we study a family of recursive visualization schemes for many-particle sys- tems, for which we have coined the name “qubism”. While all many-qudit states can be plotted with qubism, it is especially useful for spin chains and one-dimensional trans- lationally invariant states. This symmetry results in self-similarity of the plot, making it more comprehensible and allowing to discover certain structures. This visualization scheme allows to compare states of different particle numbers (which may be useful in numerical simulations when particle number is an open parameter) and puts emphasis on correlations between neighboring particles. The visualization scheme can be used to plot probability distribution of sequences, e.g. related to series of nucleotides in RNA and DNA or — aminoacids in proteins. However, unlike classical probabilistic ensem- bles of sequences, visualizing quantum states offers more — showing entanglement and allowing to observe quantum phase transitions. Third, we study quantum walks of a single particle on graphs, which are classical analogues of random walks. Our focus in on the long-time limit of the probability v Contents distribution. We define “quantumness” to be the difference between the probability distributions of the quantum and related classical random walks. Moreover, we study how (especially in the long-time limit) off-diagonal elements of the density matrix behave. That is, we measure coherence between different nodes, and we use them to perform quantum community detection — splitting of a graph into subgraphs in such a way that the coherence between them is small. We perform a bottom-up hierarchical aggregation, with a scheme similar to modularity maximization, which is a standard tool for the, so called, community detection for (classical) complex networks. However, our method captures properties that classical methods cannot — the impact of constructive and destructive interference, as well as the dependence of the results on the tunneling phase. vi 1 Introduction 1.1 Background This PhD thesis is divided into three chapters, each one describing a distinct thread of research: Invariants for bosonic and symmetric states, • Visualization of quantum sequences, • Quantum walks on complex networks. • Yet, these threads are connected through common concepts and methods related to study of entanglement, symmetry with respect to interchange of particles and self- similarity of quantum systems. An illustrative graph of these concepts and their rela- tions is depicted in Fig. 1.1. 1.1.1 Historical background While quantum mechanics dates back to the beginning of the 20th century [2], only in the last decades it was started being considered as a tool for information processing, that is, communication, cryptography and computation [3]. One of the first research in that line was a study of the information capacity for transmission of information with quantum [4] rather than classical states [5]. The result, now known as Holevo’s bound, is that regardless if we are operating with classical or quantum d-level states, we can transmit up to log2(d) bits of information per state. That is, for this particular task there is no advantage of using quantum states over classical states.
Recommended publications
  • Quantum Theory Cannot Consistently Describe the Use of Itself

    Quantum Theory Cannot Consistently Describe the Use of Itself

    ARTICLE DOI: 10.1038/s41467-018-05739-8 OPEN Quantum theory cannot consistently describe the use of itself Daniela Frauchiger1 & Renato Renner1 Quantum theory provides an extremely accurate description of fundamental processes in physics. It thus seems likely that the theory is applicable beyond the, mostly microscopic, domain in which it has been tested experimentally. Here, we propose a Gedankenexperiment 1234567890():,; to investigate the question whether quantum theory can, in principle, have universal validity. The idea is that, if the answer was yes, it must be possible to employ quantum theory to model complex systems that include agents who are themselves using quantum theory. Analysing the experiment under this presumption, we find that one agent, upon observing a particular measurement outcome, must conclude that another agent has predicted the opposite outcome with certainty. The agents’ conclusions, although all derived within quantum theory, are thus inconsistent. This indicates that quantum theory cannot be extrapolated to complex systems, at least not in a straightforward manner. 1 Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland. Correspondence and requests for materials should be addressed to R.R. (email: [email protected]) NATURE COMMUNICATIONS | (2018) 9:3711 | DOI: 10.1038/s41467-018-05739-8 | www.nature.com/naturecommunications 1 ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-05739-8 “ 1”〉 “ 1”〉 irect experimental tests of quantum theory are mostly Here, | z ¼À2 D and | z ¼þ2 D denote states of D depending restricted to microscopic domains. Nevertheless, quantum on the measurement outcome z shown by the devices within the D “ψ ”〉 “ψ ”〉 theory is commonly regarded as being (almost) uni- lab.
  • On the Linear Stability of Plane Couette Flow for an Oldroyd-B Fluid and Its

    On the Linear Stability of Plane Couette Flow for an Oldroyd-B Fluid and Its

    J. Non-Newtonian Fluid Mech. 127 (2005) 169–190 On the linear stability of plane Couette flow for an Oldroyd-B fluid and its numerical approximation Raz Kupferman Institute of Mathematics, The Hebrew University, Jerusalem, 91904, Israel Received 14 December 2004; received in revised form 10 February 2005; accepted 7 March 2005 Abstract It is well known that plane Couette flow for an Oldroyd-B fluid is linearly stable, yet, most numerical methods predict spurious instabilities at sufficiently high Weissenberg number. In this paper we examine the reasons which cause this qualitative discrepancy. We identify a family of distribution-valued eigenfunctions, which have been overlooked by previous analyses. These singular eigenfunctions span a family of non- modal stress perturbations which are divergence-free, and therefore do not couple back into the velocity field. Although these perturbations decay eventually, they exhibit transient amplification during which their “passive" transport by shearing streamlines generates large cross- stream gradients. This filamentation process produces numerical under-resolution, accompanied with a growth of truncation errors. We believe that the unphysical behavior has to be addressed by fine-scale modelling, such as artificial stress diffusivity, or other non-local couplings. © 2005 Elsevier B.V. All rights reserved. Keywords: Couette flow; Oldroyd-B model; Linear stability; Generalized functions; Non-normal operators; Stress diffusion 1. Introduction divided into two main groups: (i) numerical solutions of the boundary value problem defined by the linearized system The linear stability of Couette flow for viscoelastic fluids (e.g. [2,3]) and (ii) stability analysis of numerical methods is a classical problem whose study originated with the pi- for time-dependent flows, with the objective of understand- oneering work of Gorodtsov and Leonov (GL) back in the ing how spurious solutions emerge in computations.
  • A Highly Efficient Spectral-Galerkin Method Based on Tensor Product for Fourth-Order Steklov Equation with Boundary Eigenvalue

    A Highly Efficient Spectral-Galerkin Method Based on Tensor Product for Fourth-Order Steklov Equation with Boundary Eigenvalue

    An et al. Journal of Inequalities and Applications (2016)2016:211 DOI 10.1186/s13660-016-1158-1 R E S E A R C H Open Access A highly efficient spectral-Galerkin method based on tensor product for fourth-order Steklov equation with boundary eigenvalue Jing An1,HaiBi1 and Zhendong Luo2* *Correspondence: [email protected] Abstract 2School of Mathematics and Physics, North China Electric Power In this study, a highly efficient spectral-Galerkin method is posed for the fourth-order University, No. 2, Bei Nong Road, Steklov equation with boundary eigenvalue. By making use of the spectral theory of Changping District, Beijing, 102206, compact operators and the error formulas of projective operators, we first obtain the China Full list of author information is error estimates of approximative eigenvalues and eigenfunctions. Then we build a 1 ∩ 2 available at the end of the article suitable set of basis functions included in H0() H () and establish the matrix model for the discrete spectral-Galerkin scheme by adopting the tensor product. Finally, we use some numerical experiments to verify the correctness of the theoretical results. MSC: 65N35; 65N30 Keywords: fourth-order Steklov equation with boundary eigenvalue; spectral-Galerkin method; error estimates; tensor product 1 Introduction Increasing attention has recently been paid to numerical approximations for Steklov equa- tions with boundary eigenvalue, arising in fluid mechanics, electromagnetism, etc. (see, e.g.,[–]). However, most existing work usually treated the second-order Steklov equa- tions with boundary eigenvalue and there are relatively few articles treating fourth-order ones. The fourth-order Steklov equations with boundary eigenvalue also have been used in both mathematics and physics, for example, the main eigenvalues play a very key role in the positivity-preserving properties for the biharmonic-operator under the border conditions w = w – χwν =on∂ (see []).
  • Chebyshev and Fourier Spectral Methods 2000

    Chebyshev and Fourier Spectral Methods 2000

    Chebyshev and Fourier Spectral Methods Second Edition John P. Boyd University of Michigan Ann Arbor, Michigan 48109-2143 email: [email protected] http://www-personal.engin.umich.edu/jpboyd/ 2000 DOVER Publications, Inc. 31 East 2nd Street Mineola, New York 11501 1 Dedication To Marilyn, Ian, and Emma “A computation is a temptation that should be resisted as long as possible.” — J. P. Boyd, paraphrasing T. S. Eliot i Contents PREFACE x Acknowledgments xiv Errata and Extended-Bibliography xvi 1 Introduction 1 1.1 Series expansions .................................. 1 1.2 First Example .................................... 2 1.3 Comparison with finite element methods .................... 4 1.4 Comparisons with Finite Differences ....................... 6 1.5 Parallel Computers ................................. 9 1.6 Choice of basis functions .............................. 9 1.7 Boundary conditions ................................ 10 1.8 Non-Interpolating and Pseudospectral ...................... 12 1.9 Nonlinearity ..................................... 13 1.10 Time-dependent problems ............................. 15 1.11 FAQ: Frequently Asked Questions ........................ 16 1.12 The Chrysalis .................................... 17 2 Chebyshev & Fourier Series 19 2.1 Introduction ..................................... 19 2.2 Fourier series .................................... 20 2.3 Orders of Convergence ............................... 25 2.4 Convergence Order ................................. 27 2.5 Assumption of Equal Errors ...........................
  • A Critic Looks at Qbism Guido Bacciagaluppi

    A Critic Looks at Qbism Guido Bacciagaluppi

    A Critic Looks at QBism Guido Bacciagaluppi To cite this version: Guido Bacciagaluppi. A Critic Looks at QBism. 2013. halshs-00996289 HAL Id: halshs-00996289 https://halshs.archives-ouvertes.fr/halshs-00996289 Preprint submitted on 26 May 2014 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A Critic Looks at QBism Guido Bacciagaluppi∗ 30 April 2013 Abstract This chapter comments on that by Chris Fuchs on qBism. It presents some mild criticisms of this view, some based on the EPR and Wigner’s friend scenarios, and some based on the quantum theory of measurement. A few alternative suggestions for implementing a sub- jectivist interpretation of probability in quantum mechanics conclude the chapter. “M. Braque est un jeune homme fort audacieux. [...] Il m´eprise la forme, r´eduit tout, sites et figures et maisons, `ades sch´emas g´eom´etriques, `ades cubes. Ne le raillons point, puisqu’il est de bonne foi. Et attendons.”1 Thus commented the French art critic Louis Vauxcelles on Braque’s first one- man show in November 1908, thereby giving cubism its name. Substituting spheres and tetrahedra for cubes might be more appropriate if one wishes to apply the characterisation to qBism — the view of quantum mechanics and the quantum state developed by Chris Fuchs and co-workers (for a general reference see either the paper in this volume, or Fuchs (2010)).
  • Isospectral Families of High-Order Systems∗

    Isospectral Families of High-Order Systems∗

    Isospectral Families of High-order Systems∗ Peter Lancaster† Department of Mathematics and Statistics University of Calgary Calgary T2N 1N4, Alberta, Canada and Uwe Prells School of Mechanical, Materials, Manufacturing Engineering & Management University of Nottingham Nottingham NG7 2RD United Kingdom November 28, 2006 Keywords: Matrix Polynomials, Isospectral Systems, Structure Preserving Transformations, Palindromic Symmetries. Abstract Earlier work of the authors concerning the generation of isospectral families of second or- der (vibrating) systems is generalized to higher-order systems (with no spectrum at infinity). Results and techniques are developed first for systems without symmetries, then with Hermi- tian symmetry and, finally, with palindromic symmetry. The construction of linearizations which retain such symmetries is discussed. In both cases, the notion of strictly isospectral families of systems is introduced - implying that properties of both the spectrum and the sign-characteristic are preserved. Open questions remain in the case of strictly isospectral families of palindromic systems. Intimate connections between Hermitian and unitary sys- tems are discussed in an Appendix. ∗The authors are grateful to an anonymous reviewer for perceptive comments which led to significant improve- ments in exposition. †The research of Peter Lancaster was partly supported by a grant from the Natural Sciences and Engineering Research Council of Canada 1 1 Introduction ` By a high-order system we mean an n × n matrix polynomial L(λ) = A`λ + ··· + A1λ + A0 n×n with coefficients in C . Such a polynomial with A` 6= 0 is said to have degree `, (sometimes referred to as “order” `) and “high-order” implies that ` ≥ 2. It will be assumed throughout that the leading coefficient A` is nonsingular.
  • Contemporary Physics Uncollapsing the Wavefunction by Undoing

    Contemporary Physics Uncollapsing the Wavefunction by Undoing

    This article was downloaded by: [University of California, Riverside] On: 18 February 2010 Access details: Access Details: [subscription number 918975371] Publisher Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House, 37- 41 Mortimer Street, London W1T 3JH, UK Contemporary Physics Publication details, including instructions for authors and subscription information: http://www.informaworld.com/smpp/title~content=t713394025 Uncollapsing the wavefunction by undoing quantum measurements Andrew N. Jordan a; Alexander N. Korotkov b a Department of Physics and Astronomy, University of Rochester, Rochester, New York, USA b Department of Electrical Engineering, University of California, Riverside, CA, USA First published on: 12 February 2010 To cite this Article Jordan, Andrew N. and Korotkov, Alexander N.(2010) 'Uncollapsing the wavefunction by undoing quantum measurements', Contemporary Physics, 51: 2, 125 — 147, First published on: 12 February 2010 (iFirst) To link to this Article: DOI: 10.1080/00107510903385292 URL: http://dx.doi.org/10.1080/00107510903385292 PLEASE SCROLL DOWN FOR ARTICLE Full terms and conditions of use: http://www.informaworld.com/terms-and-conditions-of-access.pdf This article may be used for research, teaching and private study purposes. Any substantial or systematic reproduction, re-distribution, re-selling, loan or sub-licensing, systematic supply or distribution in any form to anyone is expressly forbidden. The publisher does not give any warranty express or implied or make any representation that the contents will be complete or accurate or up to date. The accuracy of any instructions, formulae and drug doses should be independently verified with primary sources. The publisher shall not be liable for any loss, actions, claims, proceedings, demand or costs or damages whatsoever or howsoever caused arising directly or indirectly in connection with or arising out of the use of this material.
  • Reformulation of Quantum Mechanics and Strong Complementarity from Bayesian Inference Requirements

    Reformulation of Quantum Mechanics and Strong Complementarity from Bayesian Inference Requirements

    Reformulation of quantum mechanics and strong complementarity from Bayesian inference requirements William Heartspring Abstract: This paper provides an epistemic reformulation of quantum mechanics (QM) in terms of inference consistency requirements of objective Bayesianism, which include the principle of maximum entropy under physical constraints. Physical constraints themselves are understood in terms of consistency requirements. The by-product of this approach is that QM must additionally be understood as providing the theory of theories. Strong complementarity - that dierent observers may live in separate Hilbert spaces - follows as a consequence, which resolves the rewall paradox. Other clues pointing to this reformu- lation are analyzed. The reformulation, with the addition of novel transition probability arithmetic, resolves the measurement problem completely, thereby eliminating subjectivity of measurements from quantum mechanics. An illusion of collapse comes from Bayesian updates by observer's continuous outcome data. Dark matter and dark energy pop up directly as entropic tug-of-war in the reformulation. Contents 1 Introduction1 2 Epistemic nature of quantum mechanics2 2.1 Area law, quantum information and spacetime2 2.2 State vector from objective Bayesianism5 3 Consequences of the QM reformulation8 3.1 Basis, decoherence and causal diamond complementarity 12 4 Spacetime from entanglement 14 4.1 Locality: area equals mutual information 14 4.2 Story of Big Bang cosmology 14 5 Evidences toward the reformulation 14 5.1 Quantum redundancy 15 6 Transition probability arithmetic: resolving the measurement problem completely 16 7 Conclusion 17 1 Introduction Hamiltonian formalism and Schrödinger picture of quantum mechanics are assumed through- out the writing, with time t 2 R. Whenever the word entropy is mentioned without addi- tional qualication, it refers to von Neumann entropy.
  • Scientific American

    Scientific American

    Take the Nature Publishing Group survey for the chance to win a MacBookFind out Air more Physics Scientific American (June 2013), 308, 46-51 Published online: 14 May 2013 | doi:10.1038/scientificamerican0613-46 Quantum Weirdness? It's All in Your Mind Hans Christian von Baeyer A new version of quantum theory sweeps away the bizarre paradoxes of the microscopic world. The cost? Quantum information exists only in your imagination CREDIT: Caleb Charland In Brief Quantum mechanics is an incredibly successful theory but one full of strange paradoxes. A recently developed model called Quantum Bayesianism (or QBism) combines quantum theory with probability theory in an effort to eliminate the paradoxes or put them in a less troubling form. QBism reimagines the entity at the heart of quantum paradoxes—the wave function. Scientists use wave functions to calculate the probability that a particle will have a certain property, such as being in one place and not another. But paradoxes arise when physicists assume that a wave function is real. QBism maintains that the wave function is solely a mathematical tool that an observer uses to assign his or her personal belief that a quantum system will have a specific property. In this conception, the wave function does not exist in the world—rather it merely reflects an individual's subjective mental state. ADDITIONAL IMAGES AND ILLUSTRATIONS [QUANTUM PHILOSOPHY] Four Interpretations of Quantum Mechanics [THOUGHT EXPERIMENT] The Fix for Quantum Absurdity Flawlessly accounting for the behavior of matter on scales from the subatomic to the astronomical, quantum mechanics is the most successful theory in all the physical sciences.
  • Introduction to Spectral Methods

    Introduction to Spectral Methods

    Introduction to spectral methods Eric Gourgoulhon Laboratoire de l’Univers et de ses Th´eories (LUTH) CNRS / Observatoire de Paris Meudon, France Based on a collaboration with Silvano Bonazzola, Philippe Grandcl´ement,Jean-Alain Marck & J´eromeNovak [email protected] http://www.luth.obspm.fr 4th EU Network Meeting, Palma de Mallorca, Sept. 2002 1 Plan 1. Basic principles 2. Legendre and Chebyshev expansions 3. An illustrative example 4. Spectral methods in numerical relativity 2 1 Basic principles 3 Solving a partial differential equation Consider the PDE with boundary condition Lu(x) = s(x); x 2 U ½ IRd (1) Bu(y) = 0; y 2 @U; (2) where L and B are linear differential operators. Question: What is a numerical solution of (1)-(2)? Answer: It is a function u¯ which satisfies (2) and makes the residual R := Lu¯ ¡ s small. 4 What do you mean by “small” ? Answer in the framework of Method of Weighted Residuals (MWR): Search for solutions u¯ in a finite-dimensional sub-space PN of some Hilbert space W (typically a L2 space). Expansion functions = trial functions : basis of PN : (Á0;:::;ÁN ) XN u¯ is expanded in terms of the trial functions: u¯(x) = u˜n Án(x) n=0 Test functions : family of functions (Â0;:::;ÂN ) to define the smallness of the residual R, by means of the Hilbert space scalar product: 8 n 2 f0;:::;Ng; (Ân;R) = 0 5 Various numerical methods Classification according to the trial functions Án: Finite difference: trial functions = overlapping local polynomials of low order Finite element: trial functions = local smooth functions (polynomial
  • Math 660-Lecture 15: Spectral Methods: Fourier Pseduo-Spectral Method

    Math 660-Lecture 15: Spectral Methods: Fourier Pseduo-Spectral Method

    Math 660-Lecture 15: Spectral methods: Fourier pseduo-spectral method The topics here will be in the book `Spectral methods in Matlab' (Chap- ter 3); Another reference book is `Spectral methods in chemistry and physics.' 1 Spectral methods v.s. pseudo-spectral methods 1.1 Spectral method A spectral method is to represent the solution of a differential equation in terms of a basis of some vector space and then reduce the differential equations to an ODE system for the coefficients. In this sense, the (contin- uous) Fourier Transform and the Finite Element methods are all spectral methods. In particular, suppose we have some equation T (u) = f(x) where T is some temporal-spatial operator for the differential equation. Then, under the basis we have the error: X X R(x) = T ( ci(t)φi(x)) − f~iφi i i Pick test functions fχjg, which is the dual basis of fφi(x)g (if the space is a Hilbert space like in FEM, these two are identical) and we require hχj;Ri = 0: This then gives a system of equations for cj. • Advantages: Reduce differential operators to multiplications, and solv- ing ODEs might be easier. • Disadvatange: The multiplication in the spatial variable becomes con- volution. Then, we may have a big matrix-vector multiplication. 1.2 Pseudo-spectral method Here, represent the solution in terms of the basis but impose the equa- tion only at discrete points. As we can see, the pseudo-spectral method in some sense is a spectral method in discrete space. There are some typical pseudo-spectral methods 1 • Fourier Pseudo-spectral method with trigonometric basis • Chebyshev pseudo-spectral method with Chebyshev polynomial basis • Collocation methods with polynomial basis (especially for ODEs solvers); (Some people also call the pseudo-spectral method the collocation method.) For discrete points, we often choose χj = δ(x − xj), so that we require R(xj) = 0: The pseudo-spectral method may have efficient algorithms, which makes solving PDEs fast.
  • Validated and Numerically Efficient Chebyshev Spectral Methods For

    Validated and Numerically Efficient Chebyshev Spectral Methods For

    Validated and numerically efficient Chebyshev spectral methods for linear ordinary differential equations Florent Bréhard, Nicolas Brisebarre, Mioara Joldes To cite this version: Florent Bréhard, Nicolas Brisebarre, Mioara Joldes. Validated and numerically efficient Chebyshev spectral methods for linear ordinary differential equations. ACM Transactions on Mathematical Soft- ware, Association for Computing Machinery, 2018, 44 (4), pp.44:1-44:42. 10.1145/3208103. hal- 01526272v3 HAL Id: hal-01526272 https://hal.archives-ouvertes.fr/hal-01526272v3 Submitted on 5 Jun 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Validated and numerically efficient Chebyshev spectral methods for linear ordinary differential equations Florent Bréhard∗ Nicolas Brisebarrey Mioara Joldeşz Abstract In this work we develop a validated numerics method for the solu- tion of linear ordinary differential equations (LODEs). A wide range of algorithms (i.e., Runge-Kutta, collocation, spectral methods) exist for nu- merically computing approximations of the solutions. Most of these come with proofs of asymptotic convergence, but usually, provided error bounds are non-constructive. However, in some domains like critical systems and computer-aided mathematical proofs, one needs validated effective error bounds. We focus on both the theoretical and practical complexity anal- ysis of a so-called a posteriori quasi-Newton validation method, which mainly relies on a fixed-point argument of a contracting map.