Applications of the Hodge Decomposition to Biological Structure and Function Modeling

Total Page:16

File Type:pdf, Size:1020Kb

Applications of the Hodge Decomposition to Biological Structure and Function Modeling Applications of the Hodge Decomposition to Biological Structure and Function Modeling Andrew Gillette joint work with Chandrajit Bajaj and John Luecke Department of Mathematics, Institute of Computational Engineering and Sciences University of Texas at Austin, Austin, Texas 78712, USA http://www.math.utexas.edu/users/agillette university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 1/32 University Introduction Molecular dynamics are governed by electrostatic forces of attraction and repulsion. These forces are described as the solutions of a PDE over the molecular surfaces. Molecular surfaces may have complicated topological features affecting the solution. The Hodge Decomposition relates topological properties of the surface to solution spaces of PDEs over the surface. ∼ (space of forms) = (solutions to ∆u = f 6≡ 0) ⊕ (non-trivial deRham classes)university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 2/32 University Outline 1 The Hodge Decomposition for smooth differential forms 2 The Hodge Decomposition for discrete differential forms 3 Applications of the Hodge Decomposition to biological modeling university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 3/32 University Outline 1 The Hodge Decomposition for smooth differential forms 2 The Hodge Decomposition for discrete differential forms 3 Applications of the Hodge Decomposition to biological modeling university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 4/32 University Differential Forms Let Ω denote a smooth n-manifold and Tx (Ω) the tangent space of Ω at x. A k-form ω is a mapping from Ω to the space of alternating k-tensors on the tangent space of Ω at the input point: k k ∗ ω : Ω → Λ [Tx (Ω) ], ω(x) : Tx (Ω) → R, Mi=1 where ω(x) is an alternating k-tensor. The space of k-forms is denoted Λk (Ω). Alternate Characterization of k-forms A k-form represents an intrinsically k-dimensional phenomena and can be integrated over a k-dimensional region. q u E B electric electric fields are magnetic fields are charge density is potential is valued based on a are dual to electric fields valued over a volumeuniversity-logo point−valued linear current flow and valued on planes Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 5/32 University The Exterior Derivative Operator The exterior derivative operator denoted by d is d : Λk (Ω) → Λk+1(Ω), defined as follows. Let I := {i1,..., ik } denote an increasing sequence of k indices (ij < ij+1) and let dxI = dx1 ∧···∧ dxk . Given ω = I aIdxI define P ∂aI dω := daI ∧ dxI where daI := dxi ∂xi XI Xi∈I The exterior derivative generalizes the familiar operators grad, curl, and div for Ω= R3: d d d 0 −−−−−→ Λ0(Ω) −−−−−→ Λ1(Ω) −−−−−→ Λ2(Ω) −−−−−→ Λ3(Ω) −−−−−→ 0 grad curl div Example: Given f ∈ Λ0(R3), we have f : R3 → R. Accordingly: ∂f ∂f ∂f ∂f ∂f ∂f grad f := , , = dx + dy + dz =: df „ ∂x ∂y ∂z « ∂x ∂y ∂z university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 6/32 University Stokes’ Theorem Stokes’ Theorem Given a compact, oriented n-dimensional manifold Ω with boundary ∂Ω and a smooth (n − 1) form ω on Ω, the following equality holds: ω = dω Z∂Ω ZΩ Alternate Characterization of d Let ω be a k-form on a compact oriented n-manifold Ω (0 ≤ k < n). Then dω is the unique (k + 1)-form such that on any (k + 1)-dimensional submanifold Π ⊂ Ω the following holds dω = ω ZΠ Z∂Π university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 7/32 University Inner Products We would like to construct an operator going the other way in the sequence: δ δ δ δ 0 ←−−−−− Λ0(Ω) ←−−−−− Λ1(Ω) ←−−−−− Λ2(Ω) ←−−−−− Λ3(Ω) ←−−−−− 0 To do so, we must first fix an inner product on our space. An inner product (·, ·)V on a vector space V is linear in each variable, symmetric, and non-degenerate, i.e. given α ∈ V : (α, β)= 0 ∀β ⇐⇒ α = 0. An inner product (·, ·)V on V extends to an inner product (·, ·) on Λ(V ) Let {e1,..., en} be an orthonormal basis of V . We say that elements of the form ei1 ∧ . ∧ eik have grading k. Define the inner product of elements with different gradings to be zero. Define (ei1 , ej1 )V · · · (ei1 , ejk )V e e e e det 0 . .. 1 1 0 1 ( i1 ∧ . ∧ ik , j1 ∧ . ∧ jk ) := . ∈ {− , , } B (e , e ) · · · (e , e ) C @ ik j1 V ik jk V A Extend the above definitions bilinearly for arbitrary elements of Λ(V ). university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 8/32 University The Hodge Star Operator The Hodge Star operator denoted by ⋆ is ⋆ : Λk (Ω) → Λn−k (Ω). Let {e1,..., en} be an orthonormal basis for V with positive orientation. For 0 < k < n, let σ ∈ Sn satisfy σ(1) < · · · <σ(k) and σ(k + 1) < · · · <σ(n). Define ⋆ eσ(1) ∧ . ∧ eσ(k) := sign(σ) eσ(k+1) ∧ . ∧ eσ(n) . ` ´ ` ´ For k = 0 and k = n, ⋆ is defined by ⋆(1)= ±e1 ∧···∧ en and ⋆(e1 ∧···∧ en)= ±1. Example: ⋆dxdz = ⋆dx1 ∧ dx3 = sign (1 3 2)dx2 = −dx2 = −dy Alternate Characterization of ⋆ The Hodge Star is the unique operator satisfying the relationship α ∧ ⋆β =(α, β)µ ∀α, β ∈ Λk (V ), where µ := e1 ∧···∧ en is the volume element of Λ(V ). university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 9/32 University The Codifferential Operator The (exterior) codifferential operator denoted by δ is δ : Λk (Ω) → Λk−1(Ω) defined as the formal adjoint to d with respect to the inner product (·, ·) on Λ(Ω). That is, given η ∈ Λk (Ω), δη is the unique element of Λk−1(Ω) satisfying (dω, η)=(ω,δη) for all ω ∈ Λk−1(Ω). Alternate Characterization of δ Given an inner product on Λ(V ), where V has dimension n, the action of δ on a k-form is given by δ =(−1)nk+1 ⋆ d⋆ Example: δ(xyzdxdz)=(−1)3(2)+1 ⋆ d(−xyzdy)= ⋆d(xyzdy) = ⋆(yzdxdy − xydydz)= yzdz − xydx university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 10/32 University The Laplacian Operator The Laplacian operator denoted by ∆ is ∆ : Λk (Ω) → Λk (Ω) given by ∆ := dδ + δd Note that ∆ reduces to the typical Laplacian operator when k = 0. By the definitions, d0 = ∇ and δ1 = div. Hence n ∂2 ∆0 = d 1δ0 + δ1d0 = δ1d0 = div∇ = ∇·∇ = . − ∂x 2 Xi=1 i Therefore, understanding solutions to the differential form problem ∆ω = α for a given form α will shed light on the Laplace and Poisson PDE problems. university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 11/32 University The Hodge Decomposition The space of harmonic k-forms is Hk := {ω ∈ Λk (Ω) : ∆ω = 0} The Hodge Decomposition Theorem For each integer k with 0 ≤ k ≤ n, Hk is finite dimensional and we have the following orthogonal direct sum decompositions of Λk (Ω) of smooth k-forms on Ω: Λk (Ω) = ∆(Λk ) ⊕ Hk = dδ(Λk ) ⊕ δd(Λk ) ⊕ Hk = d(Λk−1) ⊕ δ(Λk+1) ⊕ Hk . Consequently, the equation ∆ω = α has a solution ω ∈ Λk (Ω) if and only if the k-form α is orthogonal to the space of harmonic k-forms. university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 12/32 University Relation to Cohomology Recall that the deRham complex on Ω is d d d − 0 −−−−−→ Λ0(Ω) −−−−−→0 Λ1(Ω) −−−−−→·1 · · −−−−−→n 1 Λn(Ω) −−−−−→ 0. k The kth deRham cohomology group is denoted HdR(Ω) and is given by k HdR := ker dk /imdk−1 Lemma Let Ω be a compact oriented Riemannian n-manifold Ω without boundary and 0 ≤ k ≤ n. Then k k HdR(Ω) =∼ H . Abbreviated Proof: k k−1 k k−1 k HdR = d(Λ ) ⊕ H /d(Λ ) =∼ H h i university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 13/32 University Summary of Smooth Hodge Decomposition Facts Summary of results in the smooth setting: k-forms are meant to be integrated over k dimensional regions d : Λk → Λk+1 is characterized by Stokes’ theorem ⋆ : Λk → Λn−k is characterized by α ∧ ⋆β =(α, β)µ δ : Λk → Λk−1 is characterized by δ =(−1)nk+1 ⋆ d⋆ The Hodge Decomposition is Λk = dΛk−1 ⊕ δΛk+1 ⊕ Hk Solutions to the Poisson equation ∆ω = α lie in dΛk−1 ⊕ δΛk+1 Solutions to the Laplace equation ∆ω = 0 are exactly the harmonic forms Hk Each k-cohomology class of Ω has a unique harmonic k-form representative. How can we recreate these results on a discrete level for computational modeling? d0 d1 d2 d3 smooth: 0 o Λ0 o / Λ1 o / Λ2 o / Λ3 / 0 δ0 O δ1 O δ2 O δ3 O discrete: 0 / 0 / 1 / 2 / 3 / 0 C C C C university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 14/32 University Outline 1 The Hodge Decomposition for smooth differential forms 2 The Hodge Decomposition for discrete differential forms 3 Applications of the Hodge Decomposition to biological modeling university-logo Computational Visualization Center , I C E S ( DepartmentThe University of Mathematics, of Texas atInstitute Austin of Computational EngineeringDec2008 and Sciences 15/32 University Discrete Differential Forms k-forms are meant to be integrated over k dimensional regions Let T be a triangulation of a smooth compact n-manifold Ω.
Recommended publications
  • General Relativity Fall 2018 Lecture 10: Integration, Einstein-Hilbert Action
    General Relativity Fall 2018 Lecture 10: Integration, Einstein-Hilbert action Yacine Ali-Ha¨ımoud (Dated: October 3, 2018) σ σ HW comment: T µν antisym in µν does NOT imply that Tµν is antisym in µν. Volume element { Consider a LICS with primed coordinates. The 4-volume element is dV = d4x0 = 0 0 0 0 dx0 dx1 dx2 dx3 . If we change coordinates, we have 0 ! @xµ d4x0 = det d4x: (1) µ @x Now, the metric components change as 0 0 0 0 @xµ @xν @xµ @xν g = g 0 0 = η 0 0 ; (2) µν @xµ @xν µ ν @xµ @xν µ ν since the metric components are ηµ0ν0 in the LICS. Seing this as a matrix operation and taking the determinant, we find 0 !2 @xµ det(g ) = − det : (3) µν @xµ Hence, we find that the 4-volume element is q 4 p 4 p 4 dV = − det(gµν )d x ≡ −g d x ≡ jgj d x: (4) The integral of a scalar function f is well defined: given any coordinate system (even if only defined locally): Z Z f dV = f pjgj d4x: (5) We can only define the integral of a scalar function. The integral of a vector or tensor field is mean- ingless in curved spacetime. Think of the integral as a sum. To sum vectors, you need them to belong to the same vector space. There is no common vector space in curved spacetime. Only in flat spacetime can we define such integrals. First parallel-transport the vector field to a single point of spacetime (it doesn't matter which one).
    [Show full text]
  • Arxiv:2009.07259V1 [Math.AP] 15 Sep 2020
    A GEOMETRIC TRAPPING APPROACH TO GLOBAL REGULARITY FOR 2D NAVIER-STOKES ON MANIFOLDS AYNUR BULUT AND KHANG MANH HUYNH Abstract. In this paper, we use frequency decomposition techniques to give a direct proof of global existence and regularity for the Navier-Stokes equations on two-dimensional Riemannian manifolds without boundary. Our techniques are inspired by an approach of Mattingly and Sinai [15] which was developed in the context of periodic boundary conditions on a flat background, and which is based on a maximum principle for Fourier coefficients. The extension to general manifolds requires several new ideas, connected to the less favorable spectral localization properties in our setting. Our argu- ments make use of frequency projection operators, multilinear estimates that originated in the study of the non-linear Schr¨odingerequation, and ideas from microlocal analysis. 1. Introduction Let (M; g) be a closed, oriented, connected, compact smooth two-dimensional Riemannian manifold, and let X(M) denote the space of smooth vector fields on M. We consider the incompressible Navier-Stokes equations on M, with viscosity coefficient ν > 0, @ U + div (U ⊗ U) − ν∆ U = − grad p in M t M ; (1) div U = 0 in M with initial data U0 2 X(M); where I ⊂ R is an open interval, and where U : I ! X(M) and p : I × M ! R represent the velocity and pressure of the fluid, respectively. Here, the operator ∆M is any choice of Laplacian defined on vector fields on M, discussed below. arXiv:2009.07259v1 [math.AP] 15 Sep 2020 The theory of two-dimensional fluid flows on flat spaces is well-developed, and a variety of global regularity results are well-known.
    [Show full text]
  • Lecture Notes in Mathematics
    Lecture Notes in Mathematics Edited by A. Dold and B. Eckmann Subseries: Mathematisches Institut der Universit~it und Max-Planck-lnstitut for Mathematik, Bonn - voL 5 Adviser: E Hirzebruch 1111 Arbeitstagung Bonn 1984 Proceedings of the meeting held by the Max-Planck-lnstitut fur Mathematik, Bonn June 15-22, 1984 Edited by E Hirzebruch, J. Schwermer and S. Suter I IIII Springer-Verlag Berlin Heidelberg New York Tokyo Herausgeber Friedrich Hirzebruch Joachim Schwermer Silke Suter Max-Planck-lnstitut fLir Mathematik Gottfried-Claren-Str. 26 5300 Bonn 3, Federal Republic of Germany AMS-Subject Classification (1980): 10D15, 10D21, 10F99, 12D30, 14H10, 14H40, 14K22, 17B65, 20G35, 22E47, 22E65, 32G15, 53C20, 57 N13, 58F19 ISBN 3-54045195-8 Springer-Verlag Berlin Heidelberg New York Tokyo ISBN 0-387-15195-8 Springer-Verlag New York Heidelberg Berlin Tokyo CIP-Kurztitelaufnahme der Deutschen Bibliothek. Mathematische Arbeitstagung <25. 1984. Bonn>: Arbeitstagung Bonn: 1984; proceedings of the meeting, held in Bonn, June 15-22, 1984 / [25. Math. Arbeitstagung]. Ed. by E Hirzebruch ... - Berlin; Heidelberg; NewYork; Tokyo: Springer, 1985. (Lecture notes in mathematics; Vol. 1t11: Subseries: Mathematisches I nstitut der U niversit~it und Max-Planck-lnstitut for Mathematik Bonn; VoL 5) ISBN 3-540-t5195-8 (Berlin...) ISBN 0-387q5195-8 (NewYork ...) NE: Hirzebruch, Friedrich [Hrsg.]; Lecture notes in mathematics / Subseries: Mathematischee Institut der UniversitAt und Max-Planck-lnstitut fur Mathematik Bonn; HST This work ts subject to copyright. All rights are reserved, whether the whole or part of the material is concerned, specifically those of translation, reprinting, re~use of illustrations, broadcasting, reproduction by photocopying machine or similar means, and storage in data banks.
    [Show full text]
  • Vector Calculus and Multiple Integrals Rob Fender, HT 2018
    Vector Calculus and Multiple Integrals Rob Fender, HT 2018 COURSE SYNOPSIS, RECOMMENDED BOOKS Course syllabus (on which exams are based): Double integrals and their evaluation by repeated integration in Cartesian, plane polar and other specified coordinate systems. Jacobians. Line, surface and volume integrals, evaluation by change of variables (Cartesian, plane polar, spherical polar coordinates and cylindrical coordinates only unless the transformation to be used is specified). Integrals around closed curves and exact differentials. Scalar and vector fields. The operations of grad, div and curl and understanding and use of identities involving these. The statements of the theorems of Gauss and Stokes with simple applications. Conservative fields. Recommended Books: Mathematical Methods for Physics and Engineering (Riley, Hobson and Bence) This book is lazily referred to as “Riley” throughout these notes (sorry, Drs H and B) You will all have this book, and it covers all of the maths of this course. However it is rather terse at times and you will benefit from looking at one or both of these: Introduction to Electrodynamics (Griffiths) You will buy this next year if you haven’t already, and the chapter on vector calculus is very clear Div grad curl and all that (Schey) A nice discussion of the subject, although topics are ordered differently to most courses NB: the latest version of this book uses the opposite convention to polar coordinates to this course (and indeed most of physics), but older versions can often be found in libraries 1 Week One A review of vectors, rotation of coordinate systems, vector vs scalar fields, integrals in more than one variable, first steps in vector differentiation, the Frenet-Serret coordinate system Lecture 1 Vectors A vector has direction and magnitude and is written in these notes in bold e.g.
    [Show full text]
  • Mathematical Theorems
    Appendix A Mathematical Theorems The mathematical theorems needed in order to derive the governing model equations are defined in this appendix. A.1 Transport Theorem for a Single Phase Region The transport theorem is employed deriving the conservation equations in continuum mechanics. The mathematical statement is sometimes attributed to, or named in honor of, the German Mathematician Gottfried Wilhelm Leibnitz (1646–1716) and the British fluid dynamics engineer Osborne Reynolds (1842–1912) due to their work and con- tributions related to the theorem. Hence it follows that the transport theorem, or alternate forms of the theorem, may be named the Leibnitz theorem in mathematics and Reynolds transport theorem in mechanics. In a customary interpretation the Reynolds transport theorem provides the link between the system and control volume representations, while the Leibnitz’s theorem is a three dimensional version of the integral rule for differentiation of an integral. There are several notations used for the transport theorem and there are numerous forms and corollaries. A.1.1 Leibnitz’s Rule The Leibnitz’s integral rule gives a formula for differentiation of an integral whose limits are functions of the differential variable [7, 8, 22, 23, 45, 55, 79, 94, 99]. The formula is also known as differentiation under the integral sign. H. A. Jakobsen, Chemical Reactor Modeling, DOI: 10.1007/978-3-319-05092-8, 1361 © Springer International Publishing Switzerland 2014 1362 Appendix A: Mathematical Theorems b(t) b(t) d ∂f (t, x) db da f (t, x) dx = dx + f (t, b) − f (t, a) (A.1) dt ∂t dt dt a(t) a(t) The first term on the RHS gives the change in the integral because the function itself is changing with time, the second term accounts for the gain in area as the upper limit is moved in the positive axis direction, and the third term accounts for the loss in area as the lower limit is moved.
    [Show full text]
  • Differential Forms Diff Geom II, WS 2015/16
    J.M. Sullivan, TU Berlin B: Differential Forms Diff Geom II, WS 2015/16 B. DIFFERENTIAL FORMS instance, if S has k elements this gives a k-dimensional vector space with S as basis. We have already seen one-forms (covector fields) on a Given vector spaces V and W, let F be the free vector space over the set V × W. (This consists of formal sums manifold. In general, a k-form is a field of alternating k- P linear forms on the tangent spaces of a manifold. Forms ai(vi, wi) but ignores all the structure we have on the set are the natural objects for integration: a k-form can be in- V × W.) Now let R ⊂ F be the linear subspace spanned by tegrated over an oriented k-submanifold. We start with ten- all elements of the form: sor products and the exterior algebra of multivectors. (v + v0, w) − (v, w) − (v0, w), (v, w + w0) − (v, w) − (v, w0), (av, w) − a(v, w), (v, aw) − a(v, w). B1. Tensor products These correspond of course to the bilinearity conditions Recall that, if V, W and X are vector spaces, then a map we started with. The quotient vector space F/R will be the b: V × W → X is called bilinear if tensor product V ⊗ W. We have started with all possible v ⊗ w as generators and thrown in just enough relations to b(v + v0, w) = b(v, w) + b(v0, w), make the map (v, w) 7→ v ⊗ w be bilinear. b(v, w + w0) = b(v, w) + b(v, w0), The tensor product is commutative: there is a natural linear isomorphism V⊗W → W⊗V such that v⊗w 7→ w⊗v.
    [Show full text]
  • [Math.GT] 7 Mar 2013 Ae Stems Ftetinuaintnst Zero
    GEOMETRIC STRUCTURES ON THE COCHAINS OF A MANIFOLD SCOTT O. WILSON Abstract. In this paper we develop several algebraic structures on the sim- plicial cochains of a triangulated manifold that are analogues of objects in differential geometry. We study a cochain product and prove several state- ments about its convergence to the wedge product on differential forms. Also, for cochains with an inner product, we define a combinatorial Hodge star oper- ator, and describe some applications, including a combinatorial period matrix for surfaces. We show that for a particularly nice cochain inner product, these combinatorial structures converge to their continuum analogues as the mesh of a triangulation tends to zero. 1. Introduction In this paper we develop combinatorial analogues of several objects in differential and complex geometry, including the Hodge star operator and the period matrix of a Riemann surface. We define these structures on the appropriate combinatorial analogue of differential forms, namely simplicial cochains. As we recall in section 3, the two essential ingredients to the smooth Hodge star operator are Poincar´eDuality and a metric, or inner product. We’ll define the combinatorial star operator in much the same way, using both an inner product and Poincar´eDuality, expressed on cochains in the form of a (graded) commutative product. Using the inner product introduced in [6], we prove the following: Theorem 1.1. The combinatorial star operator, defined on the simplicial cochains of a triangulated Riemannian manifold, converges to the smooth Hodge star operator as the mesh of the triangulation tends to zero. We show in section 7 that, on a closed surface, this combinatorial star operator gives rise to a combinatorial period matrix and prove: arXiv:math/0505227v4 [math.GT] 7 Mar 2013 Theorem 1.2.
    [Show full text]
  • Operator-Valued Tensors on Manifolds: a Framework for Field
    Operator-Valued Tensors on Manifolds: A Framework for Field Quantization 1 2 H. Feizabadi , N. Boroojerdian ∗ 1,2Department of pure Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University of Technology, No. 424, Hafez Ave., Tehran, Iran. Abstract In this paper we try to prepare a framework for field quantization. To this end, we aim to replace the field of scalars R by self-adjoint elements of a commu- tative C⋆-algebra, and reach an appropriate generalization of geometrical concepts on manifolds. First, we put forward the concept of operator-valued tensors and extend semi-Riemannian metrics to operator valued metrics. Then, in this new geometry, some essential concepts of Riemannian geometry such as curvature ten- sor, Levi-Civita connection, Hodge star operator, exterior derivative, divergence,... will be considered. Keywords: Operator-valued tensors, Operator-Valued Semi-Riemannian Met- rics, Levi-Civita Connection, Curvature, Hodge star operator MSC(2010): Primary: 65F05; Secondary: 46L05, 11Y50. 0 Introduction arXiv:1501.05065v2 [math-ph] 27 Jan 2015 The aim of the present paper is to extend the theory of semi-Riemannian metrics to operator-valued semi-Riemannian metrics, which may be used as a framework for field quantization. Spaces of linear operators are directly related to C∗-algebras, so C∗- algebras are considered as a basic concept in this article. This paper has its roots in Hilbert C⋆-modules, which are frequently used in the theory of operator algebras, allowing one to obtain information about C⋆-algebras by studying Hilbert C⋆-modules over them. Hilbert C⋆-modules provide a natural generalization of Hilbert spaces arising when the field of scalars C is replaced by an arbitrary C⋆-algebra.
    [Show full text]
  • Hodge Theory
    HODGE THEORY PETER S. PARK Abstract. This exposition of Hodge theory is a slightly retooled version of the author's Harvard minor thesis, advised by Professor Joe Harris. Contents 1. Introduction 1 2. Hodge Theory of Compact Oriented Riemannian Manifolds 2 2.1. Hodge star operator 2 2.2. The main theorem 3 2.3. Sobolev spaces 5 2.4. Elliptic theory 11 2.5. Proof of the main theorem 14 3. Hodge Theory of Compact K¨ahlerManifolds 17 3.1. Differential operators on complex manifolds 17 3.2. Differential operators on K¨ahlermanifolds 20 3.3. Bott{Chern cohomology and the @@-Lemma 25 3.4. Lefschetz decomposition and the Hodge index theorem 26 Acknowledgments 30 References 30 1. Introduction Our objective in this exposition is to state and prove the main theorems of Hodge theory. In Section 2, we first describe a key motivation behind the Hodge theory for compact, closed, oriented Riemannian manifolds: the observation that the differential forms that satisfy certain par- tial differential equations depending on the choice of Riemannian metric (forms in the kernel of the associated Laplacian operator, or harmonic forms) turn out to be precisely the norm-minimizing representatives of the de Rham cohomology classes. This naturally leads to the statement our first main theorem, the Hodge decomposition|for a given compact, closed, oriented Riemannian manifold|of the space of smooth k-forms into the image of the Laplacian and its kernel, the sub- space of harmonic forms. We then develop the analytic machinery|specifically, Sobolev spaces and the theory of elliptic differential operators|that we use to prove the aforementioned decom- position, which immediately yields as a corollary the phenomenon of Poincar´eduality.
    [Show full text]
  • Vector Calculus Operators
    Vector Calculus Operators July 6, 2015 1 Div, Grad and Curl The divergence of a vector quantity in Cartesian coordinates is defined by: −! @ @ @ r · ≡ x + y + z (1) @x @y @z −! where = ( x; y; z). The divergence of a vector quantity is a scalar quantity. If the divergence is greater than zero, it implies a net flux out of a volume element. If the divergence is less than zero, it implies a net flux into a volume element. If the divergence is exactly equal to zero, then there is no net flux into or out of the volume element, i.e., any field lines that enter the volume element also exit the volume element. The divergence of the electric field and the divergence of the magnetic field are two of Maxwell's equations. The gradient of a scalar quantity in Cartesian coordinates is defined by: @ @ @ r (x; y; z) ≡ ^{ + |^ + k^ (2) @x @y @z where ^{ is the unit vector in the x-direction,| ^ is the unit vector in the y-direction, and k^ is the unit vector in the z-direction. Evaluating the gradient indicates the direction in which the function is changing the most rapidly (i.e., slope). If the gradient of a function is zero, the function is not changing. The curl of a vector function in Cartesian coordinates is given by: 1 −! @ @ @ @ @ @ r × ≡ ( z − y )^{ + ( x − z )^| + ( y − x )k^ (3) @y @z @z @x @x @y where the resulting function is another vector function that indicates the circulation of the original function, with unit vectors pointing orthogonal to the plane of the circulation components.
    [Show full text]
  • A New Definition of the Representative Volument Element in Numerical
    UNIVERSITE´ DE MONTREAL´ A NEW DEFINITION OF THE REPRESENTATIVE VOLUMENT ELEMENT IN NUMERICAL HOMOGENIZATION PROBLEMS AND ITS APPLICATION TO THE PERFORMANCE EVALUATION OF ANALYTICAL HOMOGENIZATION MODELS HADI MOUSSADDY DEPARTEMENT´ DE GENIE´ MECANIQUE´ ECOLE´ POLYTECHNIQUE DE MONTREAL´ THESE` PRESENT´ EE´ EN VUE DE L'OBTENTION DU DIPLOME^ DE PHILOSOPHIÆ DOCTOR (GENIE´ MECANIQUE)´ AVRIL 2013 © Hadi Moussaddy, 2013. UNIVERSITE´ DE MONTREAL´ ECOLE´ POLYTECHNIQUE DE MONTREAL´ Cette th`eseintitul´ee : A NEW DEFINITION OF THE REPRESENTATIVE VOLUMENT ELEMENT IN NUMERICAL HOMOGENIZATION PROBLEMS AND ITS APPLICATION TO THE PERFORMANCE EVALUATION OF ANALYTICAL HOMOGENIZATION MODELS pr´esent´ee par : MOUSSADDY Hadi en vue de l'obtention du dipl^ome de : Philosophiæ Doctor a ´et´ed^ument accept´eepar le jury d'examen constitu´ede : M. GOSSELIN Fr´ed´eric, Ph.D., pr´esident M. THERRIAULT Daniel, Ph.D., membre et directeur de recherche M. LE´VESQUE Martin, Ph.D., membre et codirecteur de recherche M. LESSARD Larry, Ph.D., membre M. BOHM¨ Helmut J., Techn. Doct., membre iii To my parents Safwan and Aisha To my wife Nour To my whole family... iv ACKNOWLEDGEMENTS I would like to thank Pr. Daniel Therriault and Pr. Martin L´evesque for their guidance throughout this project. Your insightful criticism guided me to complete a rigorous work, while gaining scientific maturity. Also, I want to thank Pr. B¨ohm, Pr. Lessard and Pr. Henry for being my committee members and Pr. Gosselin for being the Jury president. I owe special gratitude to the members of Laboratory for Multiscale Mechanics (LM2) research group. I would like to paricularly acknowledge the contributions of Maryam Pah- lavanPour by conducting parts of the analytical modeling included in this thesis, and for numerous fruitful dicussions all along this project.
    [Show full text]
  • Electric Flux Density, Gauss's Law, and Divergence
    www.getmyuni.com CHAPTER3 ELECTRIC FLUX DENSITY, GAUSS'S LAW, AND DIVERGENCE After drawinga few of the fields described in the previous chapter and becoming familiar with the concept of the streamlines which show the direction of the force on a test charge at every point, it is difficult to avoid giving these lines a physical significance and thinking of them as flux lines. No physical particle is projected radially outward from the point charge, and there are no steel tentacles reaching out to attract or repel an unwary test charge, but as soon as the streamlines are drawn on paper there seems to be a picture showing``something''is present. It is very helpful to invent an electric flux which streams away symmetri- cally from a point charge and is coincident with the streamlines and to visualize this flux wherever an electric field is present. This chapter introduces and uses the concept of electric flux and electric flux density to solve again several of the problems presented in the last chapter. The work here turns out to be much easier, and this is due to the extremely symmetrical problems which we are solving. www.getmyuni.com 3.1 ELECTRIC FLUX DENSITY About 1837 the Director of the Royal Society in London, Michael Faraday, became very interested in static electric fields and the effect of various insulating materials on these fields. This problem had been botheringhim duringthe past ten years when he was experimentingin his now famous work on induced elec- tromotive force, which we shall discuss in Chap. 10.
    [Show full text]