Inverse problems for hyperbolic PDEs
|
|
- Joel Adams
- 5 years ago
- Views:
Transcription
1 Inverse problems for hyperbolic PDEs Lauri Oksanen University College London
2 Example: inverse problem for the wave equation Let c be a smooth function on Ω R n and consider the wave equation t 2 u c 2 u = 0, in (0, T ) Ω, ν u (0,T ) Ω = f, u t=0 = t u t=0 = 0, and define the Neumann-to-Dirichlet map by Λf = u (0,T ) Ω. Theorem [Belishev 87]. Λ determines c uniquely for large enough T. Present challenges: Uniqueness for more general Lorentzian geometries Computational methods
3 More general Lorentzian geometries For a Lorentzian metric g on M = (0, T ) Ω, we define the wave operator g u = g 1/2 n ( x j g jk g 1/2 ) x k u. j,k=0 Writing x = (t, y) and g(x) = dt 2 + c 2 (y)dy 2, we have g u = 2 t u + c 2 (y) y u + lower order terms. Inverse problem for g. Determine g up to an isometry, given C(g) = {(u M, ν u M ); g u = 0, u C (M)}. The problem is open even when g is close to the Minkowski metric. The analogous elliptic inverse problem, with a Riemannian g close to the Euclidean metric, is also open.
4 The problem in a fixed conformal class Suppose now that g = ah where h is a known Lorentzian metric on M, and a is a strictly positive function on M. Inverse problem for a. Given h and the Cauchy data set C(g) = {(u M, ν u M ); g u = 0, u C (M)}. determine the conformal factor a. The gauge transformation, a (n+2)/4 g (a (n 2)/4 u) = h u + qu, q = a (n+2)/4 ( g a (n 2)/4 ), can be used to reduce determination of a conformal factor a to determination of a potential q.
5 Determination of a potential Suppose that h is of the product form h = dt 2 + g 0 (y), (t, y) M = (0, T ) Ω, where g 0 is a smooth Riemannian metric on Ω. Theorem [Kian-L.O. 17]. Suppose that (Ω, g 0 ) is simple. Then g 0 and C(q) = {(u M, ν u M ); h u + qu = 0, u C (M)}, determine the time-dependent potential q. Work in progress with Feizmohammadi, Ilmavirta and Kian: Non-smooth q Any g 0 such that the geodesic ray transform is injective on (Ω, g 0 ) Determination of a first order perturbation (i.e. magnetic potential)
6 Time independent problems Let (M 0, g 0 ) be a compact Riemannian manifold with boundary, and let E be a Hermitian vector bundle over M 0. Let be a connection on E and let Q C (M 0 ; End(E)). Suppose that and Q are symmetric: d u, v E = u, v E + u, v E, Qu, v E = u, Qv E for u, v C (M 0 ; E). Consider the local Cauchy data C(M 0, g 0, E,, Q) = {(u Γ, ν u Γ ); 2 t u u + Qu = 0, u M\Γ = 0}, where M = R M 0, Γ = R Γ 0, with Γ 0 M 0 open and non-empty, and is the adjoint of with respect to g 0. Theorem [Kurylev-L.O.-Paternain 18]. The set C(M 0, g 0, E,, Q) determines the structure (M 0, g 0, E,, Q) up to the gauge invariances.
7 Two proof strategies Theorem [Kian-L.O. 17]. Suppose that (M 0, g 0 ) is simple. Then g 0 and C(q) determine the time-dependent potential q. Use geometric optics to recover the light ray transform of q Reduce the light ray transform to the geodesic ray transform on (M 0, g 0 ) (uses the product structure) Invert the geodesic ray transform (uses simplicity of (M 0, g 0 )) Theorem [Kurylev-L.O.-Paternain 18]. The set C(M 0, g 0, E,, Q) determines the time-independent structure (M 0, g 0, E,, Q). The Boundary Control method [Belishev 87] based on hyperbolic unique continuation [Tataru 95] 1 1 Fails for smooth time-dependent coefficients [Alinhac 83]
8 The Minkowski space Consider R 1+3 with the Minkowski metric g = dt 2 + dy 2. As a matrix g = y θ β y,θ A vector v = (v 0, v 1, v 2, v 3 ) is lightlike if the scalar product (v, v) g = v j g jk v k vanishes, and a line (i.e. a geodesic) is called a light ray if its tangent vector is lightlike. Writing S 2 = {θ R 3 ; θ = 1}, the light rays can be parametrized by β y,θ (s) = (s, y + sθ), (y, θ) R 3 S 2. Example: taking y = 0 and θ = (1, 0, 0) gives β(s) = (s, s, 0, 0).
9 Example: geometric optics in the Minkowski space Let g be the Minkowski metric on R 1+3 and consider the ansatz u(t, x) e iσ(t x1) (a(t, x) + iσ 1 b(t, x)), σ 1. The first order terms w.r.t. σ in g u + qu = 0 vanish when a satisfies the transport equation t a + x 1a = 0. Choosing a(t, x) = h(t x 1 )χ(x 2, x 3 ), with h δ and χ δ, focuses a on the light ray β(s) = (s, s, 0, 0). The zeroth order terms vanish when t b + x 1b = g a + qa. After eliminating the known term g a in t b + x 1b = g a + qa, we can read from b the integral (qa)(s, s t + x 1, x 2, x 3 )ds q(s, s, 0, 0)ds = q. R R β
10 Example: light ray transform in the Minkowski space Using the parametrization β y,θ (s) = (s, y + sθ), (y, θ) R 3 S 2, of light rays, the light ray transform is given by Lq(y, θ) = q = q(s, y + sθ)ds, (y, θ) R 3 S 2. β y,θ R Writing (s, y + sθ) = x = (x 0, x ), we have the Fourier slicing e iη y Lq(y, θ) dy = e iη (x x0θ) q(x) dx = q( η θ, η). R 3 R 4 Here η θ η, and choosing θ as a suitable linear combination of η and a vector orthogonal to it, we recover q(r η, η) for any r [ 1, 1]. The Fourier transform q(ξ) can be recovered for spacelike and lightlike ξ, that is, for ξ = (ξ 0, ξ ) satisfying ξ 0 ξ. If q is compactly supported, it can be recovered using analytic continuation [Stefanov 89].
11 Two proof strategies and convexity For time-dependent potential q. Recover the light ray transform Lq(β) = β q, with β a light ray. Then invert L. For time-independent structure (M 0, g 0, E,, Q). Use hyperbolic unique continuation that is sharp in time.
12 Inverse problems for non-linear wave equations Consider a dimensional, globally hyperbolic Lorentzian manifold without boundary (M, g), and let V M be a bounded neighbourhood of a timelike curve µ : [0, 1] M. Moreover, let W be the intersection of the chronological future of µ(0) and past of µ(1). Theorem (IP with u 2 non-linearity) [Kurylev-Lassas-Uhlmann 17]. For any k N and r > 0 the data set D k,r (M, g) = {(u V, f ); g u + u 2 = f, u t<0 = 0, f C k 0 (V ), f r} determines (W, g) up to an isometry and a conformal rescaling.
13 Inverse problem for the Einstein-scalar field equations Theorem [Kurylev-Lassas-L.O.-Uhlmann]. Let M = R M 0 and ĝ = β(t, y)dt 2 + κ(t, y), where β is strictly positive and κ(t, ) is a Riemannian metric on M 0. Let ˆφ be scalar fields such that the Einstein-scalar field equations Ein(g) = T(g, φ) + f 1, g φ = f 2, (1) hold with (g, φ, f ) = (ĝ, ˆφ, 0), and suppose that d ˆφ l, l = 1, 2, 3, 4, are linearly independent near V. Then, for any k N and r > 0, and for a suitable pull-back to Fermi coordinates Φ g, the data set D k,r (M, ĝ, ˆφ) = Φ g {(g V, φ V, f ); (1) holds, g t<0 = ĝ, φ t<0 = ˆφ, f C k 0 (V ), f r} determines (W, ĝ) up to an isometry and a conformal rescaling.
14 Non-linearity in the Einstein-scalar field equations In the wave gauge, the Einstein-scalar field equations can be written as Ricĝ (g) dφ dφ = f 1, in M, g φ = f 2, in M, g = ĝ, φ = φ, in (, 0) M 0, where the reduced Ricci tensor Ricĝ (g) is given locally by Ricĝ (g) jk = 1 2 g pq p q g jk + g pq g rs (Γ prj Γ qsk + Γ prj Γ qks + Γ prk Γ qjs ) + j Ĥ k, Γ pkq = 1 2 ( pg kq + q g pk k g pq ), Ĥ k = g kr g pqˆγ r pq. Here j is the covariant derivative with respect to g in the coordinate direction j, and Γ r pq are the Christoffel symbols of ĝ.
15 Strategy of the proof Use again geometric optics but replace the WKB ansatz e iσφ a with conormal distributions. Linearize u(ɛ) = (g(ɛ), φ(ɛ)) with respect to a multi-dimensional family of sources f = J j=1 ɛ jf j. Here ɛ = (ɛ 1,..., ɛ J ). Choose each f j so that u j = ɛj u ɛ=0 is a geometric optics solution. Analyse the interaction of u j, j = 1,..., J, via the non-linearity, and show that singularities propagating back to V can be generated. Schematic. Geometric optics solutions u j focus near null geodesics (black curves). When they overlap, the non-linear terms cause new propagating singularities to appear (in red).
16 Example: quadratic non-linearity Consider the solution u(ɛ) = u of Then u j = ɛj u ɛ=0 satisfies and u 12 = ɛ1 ɛ2 u ɛ=0 satisfies g u u2 2 = ɛ 1f 1 + ɛ 2 f 2, u t<0 = 0. g u j = f j, u j t<0 = 0, g u 12 = u 1 u 2, u j t<0 = 0. How nice singularities in u 1 and u 2 behave under the product u 1 u 2?
17 Products of conormal distributions [Greenleaf-Uhlmann 93] Let K M be a submanifold. Recall [Hörmander 71]: a distribution 1/2-density u D (M) is in I (K) if it is smooth away from K and in local coordinates x = (x, x ) R N+(n N), with K = {x = 0}, it is of the form u, φ = R N iξ x e R n a(x, ξ )φ(x)dxdξ. We write σ[u] = a for the symbol, a 1/2-density on the conormal bundle, N K = {(x, ξ) T M; x K, ξv = 0 for v T x K} = {x = 0, ξ = 0}. Let K 1, K 2 M be transversal submanifolds, and let u j I (K j ). Then away from N K 1 N K 2 it holds that u 1 u 2 I (K 1 K 2 ) and σ[u 1 u 2 ](x, ξ) = σ[u 1 ](x, ξ (1) )σ[u 2 ](x, ξ (2) ), where ξ = ξ (1) + ξ (2) N x (K 1 K 2 ) = N x K 1 N x K 2.
18 Example: 4th order non-linearity Recall that we have assumed dim M = Let K j M, j = 1, 2, 3, 4, be of codim 1. If x 4 j=1 K j and u j I (K j ), then generically T x M = N x K 1 N x K 4, σ[u 1 u 2 u 3 u 4 ](x, ξ) = The easiest non-linearity to analyse is 4 g u u4 4! = ɛ j f j, u t<0 = 0. j=1 Writing u j = ɛj u ɛ=0 and w = ɛ1... ɛ4 u ɛ=0 we have g u j = f j, g w = u 1 u 2 u 3 u 4. 4 σ[u j ](x, ξ (j) ). Near x, outside a small exceptional set, WF (u 1 u 2 u 3 u 4 ) = WF (δ x ) = T x M. j=1
19 Summary of uniqueness results The inverse problem to recover q given the Cauchy data C(q) = {(u M, ν u M ); g u + qu = 0, u C (M)} is open on (M, g) with M = (0, T ) Ω R 1+3 and g near the Minkowski metric. Time-independent problems are well-understood in theory, however, they are computationally challenging. The method based on interaction of singularities has proven to be flexible. It works for several non-linear wave equations, and it is under active further study. Non-linearity is used in an essential way, and there are (almost) no results for linear wave equations outside the case of product geometry.
20 On computational methods Let us consider again the basic wave equation 2 t u c 2 u = 0, in (0, T ) Ω, ν u (0,T ) Ω = f, u t=0 = t u t=0 = 0, and define the Neumann-to-Dirichlet map by Λf = u (0,T ) Ω. Λ determines c uniquely for large enough T [Belishev 87] Current computational methods in geophysical imaging are based on local techniques such as various linearizations together with least squares data fitting, not on the mathematical theory Convergence to local minima is a practical problem
21 Example: an acoustic lens The speed of propagation typically increases with the depth in geophysical imaging Low velocity regions cause problems for local data fitting methods due to focusing effects [Huang-Nammour-Symes 17] Speed of sound c(x). Surface is on top.
22 Example: an acoustic lens The speed of propagation typically increases with the depth in geophysical imaging Low velocity regions cause problems for local data fitting methods due to focusing effects [Huang-Nammour-Symes 17] Geodesics emanating from a point at the surface.
23 Reconstruction of the front face of the lens Reconstruction of c given simulated measurements corresponding to 241 point like sources on the top edge [de Hoop-Kepley-L.O.]. Our method is guaranteed to converge to the true solution We use ideas from the Boundary Control method [Belishev 87]
24 Reconstruction of the front face of the lens Reconstruction of c given simulated measurements corresponding to 241 point like sources on the top edge [de Hoop-Kepley-L.O.]. Our method is guaranteed to converge to the true solution We use ideas from the Boundary Control method [Belishev 87]
25 What about the back face of the lens? Due to the focusing, our method can not immediately recover the back face of the lens In theory, the data can be extended across the region where the speed of sound is already known by using unique continuation Then the speed of sound can be reconstructed deeper down by repeating the local reconstruction step We have developed stabilized finite element methods for unique continuation problems Heat equation [Burman-L.O. 18] Helmholtz equation [Burman-Nechita-L.O.] Time domain wave equation is a work in progress with Burman and Feizmohammadi
26 Unique continuation for the Helmholtz equation Unique continuation problem. Consider sets ω B Ω. Let u be a solution to u + k 2 u = 0 in Ω. Given u ω determine u B. Due to the lack of boundary conditions, the problem is ill-posed, and the optimal stability is of conditional Hölder type. B B ω ω Two computational solutions to the unique continuation problem. Left. Exact u. Middle. Solution in a convex geometry, that is, B belongs to the convex hull of ω. Right. Solution in a non-convex geometry.
27 Stabilized finite element method Using G(u, z) = ( u, z) L 2 (Ω) k 2 (u, z) L 2 (Ω), the weak formulation of the Helmholtz equation can be written: G(u, z) = 0 for all z H 1 0 (Ω). Our stabilized FEM is obtained as the equation for critical points of L(u h, z h ) = 1 2 u h u ω 2 L 2 (ω) + G(u h, z h ) s(u h, u h ) 1 2 s (z h, z h ). We consider L on V h W h, where V h is the H 1 conformal, P 1 finite element space over a quasi-uniform triangulation of polygonal Ω, and W h = V h H0 1 (Ω). Here h is the mesh size, and s(u h, u h ) = h n u h 2 F ds + hk 2 u h 2 L 2 (Ω), F F h s (z h, z h ) = z h 2 L 2 (Ω), F where n u h F is the jump of the normal derivative across F and F h is the set of internal faces of the triangulation.
28 Convergence estimate Theorem [Burman-Nechita-L.O.]. Let ω B Ω be as in the convex case. Let (u h, z h ) V h be the critical point of the Lagrangian L, and u the solution to the unique continuation problem. There are C > 0 and α (0, 1) such that for all k 0 and h > 0 satisfying kh 1, it holds that ( ) u h u L 2 (B) Chα k 2α 2 u H 2 (Ω) + k2 u L 2 (Ω). Contrary to the previous methods for unique continuation problems, e.g. the quasi-reversibility method [Lattès-Lions 67], we regularize on the discrete spaces, and h is the only asymptotic parameter In the non-convex case, we do not derive an explicit dependence on k but we show that if there is a line that intersects B but not ω, then the estimate must blow up faster than any polynomial in k
29 Convergence: the convex case Circles: H 1 -error, rate 0.64; squares: L 2 -error, rate 0.66; triangles : h 1 J(u h, u h ), rate 1; triangles : s (z h, z h ) 1/2, rate 1.3.
30 Convergence: the effect of noise in the convex case Perturbation O(h) Perturbation O(h 2 )
Unique continuation for the Helmholtz equation using a stabilized finite element method
Unique continuation for the Helmholtz equation using a stabilized finite element method Lauri Oksanen University College London Based on a joint work with Erik Burman and Mihai Nechita Motivation: recovering
More informationHyperbolic inverse problems and exact controllability
Hyperbolic inverse problems and exact controllability Lauri Oksanen University College London An inverse initial source problem Let M R n be a compact domain with smooth strictly convex boundary, and let
More informationRecovery of anisotropic metrics from travel times
Purdue University The Lens Rigidity and the Boundary Rigidity Problems Let M be a bounded domain with boundary. Let g be a Riemannian metric on M. Define the scattering relation σ and the length (travel
More informationA SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY
A SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY PLAMEN STEFANOV 1. Introduction Let (M, g) be a compact Riemannian manifold with boundary. The geodesic ray transform I of symmetric 2-tensor fields f is
More informationThe X-ray transform for a non-abelian connection in two dimensions
The X-ray transform for a non-abelian connection in two dimensions David Finch Department of Mathematics Oregon State University Corvallis, OR, 97331, USA Gunther Uhlmann Department of Mathematics University
More informationMicrolocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries
Microlocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries David Dos Santos Ferreira LAGA Université de Paris 13 Wednesday May 18 Instituto de Ciencias Matemáticas,
More informationDiffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)
Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D
More informationON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES
ON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES CLAUS GERHARDT Abstract. We prove that the mean curvature τ of the slices given by a constant mean curvature foliation can be used
More informationLECTURE NOTES ON GEOMETRIC OPTICS
LECTURE NOTES ON GEOMETRIC OPTICS PLAMEN STEFANOV We want to solve the wave equation 1. The wave equation (1.1) ( 2 t c 2 g )u = 0, u t=0 = f 1, u t t=0 = f 2, in R t R n x, where g is a Riemannian metric
More informationAdiabatic limits and eigenvalues
Adiabatic limits and eigenvalues Gunther Uhlmann s 60th birthday meeting Richard Melrose Department of Mathematics Massachusetts Institute of Technology 22 June, 2012 Outline Introduction 1 Adiabatic metrics
More informationStability and Instability of Black Holes
Stability and Instability of Black Holes Stefanos Aretakis September 24, 2013 General relativity is a successful theory of gravitation. Objects of study: (4-dimensional) Lorentzian manifolds (M, g) which
More informationInverse Scattering with Partial data on Asymptotically Hyperbolic Manifolds
Inverse Scattering with Partial data on Asymptotically Hyperbolic Manifolds Raphael Hora UFSC rhora@mtm.ufsc.br 29/04/2014 Raphael Hora (UFSC) Inverse Scattering with Partial data on AH Manifolds 29/04/2014
More informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
More informationNull Cones to Infinity, Curvature Flux, and Bondi Mass
Null Cones to Infinity, Curvature Flux, and Bondi Mass Arick Shao (joint work with Spyros Alexakis) University of Toronto May 22, 2013 Arick Shao (University of Toronto) Null Cones to Infinity May 22,
More informationA brief introduction to Semi-Riemannian geometry and general relativity. Hans Ringström
A brief introduction to Semi-Riemannian geometry and general relativity Hans Ringström May 5, 2015 2 Contents 1 Scalar product spaces 1 1.1 Scalar products...................................... 1 1.2 Orthonormal
More informationA connection between Lorentzian distance and mechanical least action
A connection between Lorentzian distance and mechanical least action Ettore Minguzzi Università Degli Studi Di Firenze Non-commutative structures and non-relativistic (super)symmetries, LMPT Tours, June
More informationGlobal stability problems in General Relativity
Global stability problems in General Relativity Peter Hintz with András Vasy Murramarang March 21, 2018 Einstein vacuum equations Ric(g) + Λg = 0. g: Lorentzian metric (+ ) on 4-manifold M Λ R: cosmological
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
More informationTHERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING
THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING PLAMEN STEFANOV AND GUNTHER UHLMANN Abstract. We study the mathematical model of thermoacoustic and photoacoustic tomography when the sound speed has
More informationHyperbolic Geometric Flow
Hyperbolic Geometric Flow Kefeng Liu Zhejiang University UCLA Page 1 of 41 Outline Introduction Hyperbolic geometric flow Local existence and nonlinear stability Wave character of metrics and curvatures
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationPoincaré Duality Angles on Riemannian Manifolds with Boundary
Poincaré Duality Angles on Riemannian Manifolds with Boundary Clayton Shonkwiler Department of Mathematics University of Pennsylvania June 5, 2009 Realizing cohomology groups as spaces of differential
More informationQuasi-conformal minimal Lagrangian diffeomorphisms of the
Quasi-conformal minimal Lagrangian diffeomorphisms of the hyperbolic plane (joint work with J.M. Schlenker) January 21, 2010 Quasi-symmetric homeomorphism of a circle A homeomorphism φ : S 1 S 1 is quasi-symmetric
More informationON THE INVERSE PROBLEM OF FINDING COSMIC STRINGS AND OTHER TOPOLOGICAL DEFECTS
ON THE INVERSE PROBLEM OF FINDING COSMIC STRINGS AND OTHER TOPOLOGICAL DEFECTS MATTI LASSAS, LAURI OKSANEN, PLAMEN STEFANOV, AND GUNTHER UHLMANN Abstract. We consider how microlocal methods developed for
More informationNon-existence of time-periodic dynamics in general relativity
Non-existence of time-periodic dynamics in general relativity Volker Schlue University of Toronto University of Miami, February 2, 2015 Outline 1 General relativity Newtonian mechanics Self-gravitating
More informationIntroduction to Differential Geometry
More about Introduction to Differential Geometry Lecture 7 of 10: Dominic Joyce, Oxford University October 2018 EPSRC CDT in Partial Differential Equations foundation module. These slides available at
More informationÜbungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.
Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση
More informationTHE GEODESIC RAY TRANSFORM ON RIEMANNIAN SURFACES WITH CONJUGATE POINTS
THE GEODESIC RAY TRANSFORM ON RIEMANNIAN SURFACES WITH CONJUGATE POINTS FRANÇOIS MONARD, PLAMEN STEFANOV, AND GUNTHER UHLMANN Abstract. We study the geodesic X-ray transform X on compact Riemannian surfaces
More informationPseudodifferential calculus and Hadamard states
Pseudodifferential calculus and Hadamard states Local Quantum Physics and beyond - in memoriam Rudolf Haag Hamburg, Sept. 26-27 2016 Christian Gérard Département of Mathématiques Université Paris-Sud 1
More informationPractice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009.
Practice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009. Solutions (1) Let Γ be a discrete group acting on a manifold M. (a) Define what it means for Γ to act freely. Solution: Γ acts
More informationGeneral Relativity and Cosmology Mock exam
Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers
More informationLOCAL LENS RIGIDITY WITH INCOMPLETE DATA FOR A CLASS OF NON-SIMPLE RIEMANNIAN MANIFOLDS
LOCAL LENS RIGIDITY WITH INCOMPLETE DATA FOR A CLASS OF NON-SIMPLE RIEMANNIAN MANIFOLDS PLAMEN STEFANOV AND GUNTHER UHLMANN Abstract. Let σ be the scattering relation on a compact Riemannian manifold M
More informationMicrolocal Methods in X-ray Tomography
Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods
More informationAn Overview of Mathematical General Relativity
An Overview of Mathematical General Relativity José Natário (Instituto Superior Técnico) Geometria em Lisboa, 8 March 2005 Outline Lorentzian manifolds Einstein s equation The Schwarzschild solution Initial
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
More informationAn inverse source problem in optical molecular imaging
An inverse source problem in optical molecular imaging Plamen Stefanov 1 Gunther Uhlmann 2 1 2 University of Washington Formulation Direct Problem Singular Operators Inverse Problem Proof Conclusion Figure:
More informationLp Bounds for Spectral Clusters. Compact Manifolds with Boundary
on Compact Manifolds with Boundary Department of Mathematics University of Washington, Seattle Hangzhou Conference on Harmonic Analysis and PDE s (M, g) = compact 2-d Riemannian manifold g = Laplacian
More informationMathematical Relativity, Spring 2017/18 Instituto Superior Técnico
Mathematical Relativity, Spring 2017/18 Instituto Superior Técnico 1. Starting from R αβµν Z ν = 2 [α β] Z µ, deduce the components of the Riemann curvature tensor in terms of the Christoffel symbols.
More informationTHE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY
THE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY SAM KAUFMAN Abstract. The (Cauchy) initial value formulation of General Relativity is developed, and the maximal vacuum Cauchy development theorem is
More informationRECENT PROGRESS ON THE BOUNDARY RIGIDITY PROBLEM
ELECTRONIC RESEARCH ANNOUNCEMENTS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Pages 000 000 (Xxxx XX, XXXX) S 1079-6762(XX)0000-0 RECENT PROGRESS ON THE BOUNDARY RIGIDITY PROBLEM PLAMEN STEFANOV AND
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationCurved Spacetime I. Dr. Naylor
Curved Spacetime I Dr. Naylor Last Week Einstein's principle of equivalence We discussed how in the frame of reference of a freely falling object we can construct a locally inertial frame (LIF) Space tells
More information1. Geometry of the unit tangent bundle
1 1. Geometry of the unit tangent bundle The main reference for this section is [8]. In the following, we consider (M, g) an n-dimensional smooth manifold endowed with a Riemannian metric g. 1.1. Notations
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationHyperbolic Geometric Flow
Hyperbolic Geometric Flow Kefeng Liu CMS and UCLA August 20, 2007, Dong Conference Page 1 of 51 Outline: Joint works with D. Kong and W. Dai Motivation Hyperbolic geometric flow Local existence and nonlinear
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationTransport Continuity Property
On Riemannian manifolds satisfying the Transport Continuity Property Université de Nice - Sophia Antipolis (Joint work with A. Figalli and C. Villani) I. Statement of the problem Optimal transport on Riemannian
More informationThe inverse conductivity problem with power densities in dimension n 2
The inverse conductivity problem with power densities in dimension n 2 François Monard Guillaume Bal Dept. of Applied Physics and Applied Mathematics, Columbia University. June 19th, 2012 UC Irvine Conference
More informationMICROLOCAL ANALYSIS METHODS
MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis
More informationRiemannian Curvature Functionals: Lecture I
Riemannian Curvature Functionals: Lecture I Jeff Viaclovsky Park City athematics Institute July 16, 2013 Overview of lectures The goal of these lectures is to gain an understanding of critical points of
More informationRigidity of Black Holes
Rigidity of Black Holes Sergiu Klainerman Princeton University February 24, 2011 Rigidity of Black Holes PREAMBLES I, II PREAMBLE I General setting Assume S B two different connected, open, domains and
More informationCALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =
CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationTHE CALDERÓN PROBLEM AND NORMAL FORMS
THE CALDERÓN PROBLEM AND NORMAL FORMS MIKKO SALO Abstract. We outline an approach to the inverse problem of Calderón that highlights the role of microlocal normal forms and propagation of singularities
More informationThe Klein-Gordon Equation Meets the Cauchy Horizon
Enrico Fermi Institute and Department of Physics University of Chicago University of Mississippi May 10, 2005 Relativistic Wave Equations At the present time, our best theory for describing nature is Quantum
More informationSingularity formation in black hole interiors
Singularity formation in black hole interiors Grigorios Fournodavlos DPMMS, University of Cambridge Heraklion, Crete, 16 May 2018 Outline The Einstein equations Examples Initial value problem Large time
More informationMicrolocal Analysis : a short introduction
Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction
More informationA local characterization for constant curvature metrics in 2-dimensional Lorentz manifolds
A local characterization for constant curvature metrics in -dimensional Lorentz manifolds Ivo Terek Couto Alexandre Lymberopoulos August 9, 8 arxiv:65.7573v [math.dg] 4 May 6 Abstract In this paper we
More informationStrichartz Estimates in Domains
Department of Mathematics Johns Hopkins University April 15, 2010 Wave equation on Riemannian manifold (M, g) Cauchy problem: 2 t u(t, x) gu(t, x) =0 u(0, x) =f (x), t u(0, x) =g(x) Strichartz estimates:
More informationarxiv: v1 [math.dg] 24 Feb 2017
GEODESIC X-RAY TOMOGRAPHY FOR PIECEWISE CONSTANT FUNCTIONS ON NONTRAPPING MANIFOLDS JOONAS ILMAVIRTA, JERE LEHTONEN, AND MIKKO SALO arxiv:1702.07622v1 [math.dg] 24 Feb 2017 Abstract. We show that on a
More informationConstructing compact 8-manifolds with holonomy Spin(7)
Constructing compact 8-manifolds with holonomy Spin(7) Dominic Joyce, Oxford University Simons Collaboration meeting, Imperial College, June 2017. Based on Invent. math. 123 (1996), 507 552; J. Diff. Geom.
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationEigenvalues and eigenfunctions of the Laplacian. Andrew Hassell
Eigenvalues and eigenfunctions of the Laplacian Andrew Hassell 1 2 The setting In this talk I will consider the Laplace operator,, on various geometric spaces M. Here, M will be either a bounded Euclidean
More informationNon-existence of time-periodic vacuum spacetimes
Non-existence of time-periodic vacuum spacetimes Volker Schlue (joint work with Spyros Alexakis and Arick Shao) Université Pierre et Marie Curie (Paris 6) Dynamics of self-gravitating matter workshop,
More informationLECTURE 8: THE SECTIONAL AND RICCI CURVATURES
LECTURE 8: THE SECTIONAL AND RICCI CURVATURES 1. The Sectional Curvature We start with some simple linear algebra. As usual we denote by ( V ) the set of 4-tensors that is anti-symmetric with respect to
More informationOn isometries of Finsler manifolds
On isometries of Finsler manifolds International Conference on Finsler Geometry February 15-16, 2008, Indianapolis, IN, USA László Kozma (University of Debrecen, Hungary) Finsler metrics, examples isometries
More informationBlack Holes and Thermodynamics I: Classical Black Holes
Black Holes and Thermodynamics I: Classical Black Holes Robert M. Wald General references: R.M. Wald General Relativity University of Chicago Press (Chicago, 1984); R.M. Wald Living Rev. Rel. 4, 6 (2001).
More informationA REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS
A REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS DAN-ANDREI GEBA Abstract. We obtain a sharp local well-posedness result for an equation of wave maps type with variable coefficients.
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationQuasi-local mass and isometric embedding
Quasi-local mass and isometric embedding Mu-Tao Wang, Columbia University September 23, 2015, IHP Recent Advances in Mathematical General Relativity Joint work with Po-Ning Chen and Shing-Tung Yau. The
More informationThe stability of Kerr-de Sitter black holes
The stability of Kerr-de Sitter black holes András Vasy (joint work with Peter Hintz) July 2018, Montréal This talk is about the stability of Kerr-de Sitter (KdS) black holes, which are certain Lorentzian
More informationCausality, hyperbolicity, and shock formation in Lovelock theories
Causality, hyperbolicity, and shock formation in Lovelock theories Harvey Reall DAMTP, Cambridge University HSR, N. Tanahashi and B. Way, arxiv:1406.3379, 1409.3874 G. Papallo, HSR arxiv:1508.05303 Lovelock
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationSpectral theory of first order elliptic systems
Spectral theory of first order elliptic systems Dmitri Vassiliev (University College London) 24 May 2013 Conference Complex Analysis & Dynamical Systems VI Nahariya, Israel 1 Typical problem in my subject
More informationThe Geometrization Theorem
The Geometrization Theorem Matthew D. Brown Wednesday, December 19, 2012 In this paper, we discuss the Geometrization Theorem, formerly Thurston s Geometrization Conjecture, which is essentially the statement
More informationAn introduction to General Relativity and the positive mass theorem
An introduction to General Relativity and the positive mass theorem National Center for Theoretical Sciences, Mathematics Division March 2 nd, 2007 Wen-ling Huang Department of Mathematics University of
More informationInitial-Value Problems in General Relativity
Initial-Value Problems in General Relativity Michael Horbatsch March 30, 2006 1 Introduction In this paper the initial-value formulation of general relativity is reviewed. In section (2) domains of dependence,
More informationDETERMINATION OF THE SPACETIME FROM LOCAL TIME MEASUREMENTS
DETERMINATION OF THE SPACETIME FROM LOCAL TIME MEASUREMENTS MATTI LASSAS, LAURI OKSANEN, AND YANG YANG Abstract. We consider an inverse problem for a Lorentzian spacetime (M, g), and show that time measurements,
More informationA loop of SU(2) gauge fields on S 4 stable under the Yang-Mills flow
A loop of SU(2) gauge fields on S 4 stable under the Yang-Mills flow Daniel Friedan Rutgers the State University of New Jersey Natural Science Institute, University of Iceland MIT November 3, 2009 1 /
More informationStrichartz Estimates for the Schrödinger Equation in Exterior Domains
Strichartz Estimates for the Schrödinger Equation in University of New Mexico May 14, 2010 Joint work with: Hart Smith (University of Washington) Christopher Sogge (Johns Hopkins University) The Schrödinger
More informationA CONSTRUCTION OF TRANSVERSE SUBMANIFOLDS
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XLI 2003 A CONSTRUCTION OF TRANSVERSE SUBMANIFOLDS by J. Szenthe Abstract. In case of Riemannian manifolds isometric actions admitting submanifolds
More informationOn the mass of asymptotically hyperbolic manifolds
On the mass of asymptotically hyperbolic manifolds Mattias Dahl Institutionen för Matematik, Kungl Tekniska Högskolan, Stockholm April 6, 2013, EMS/DMF Joint Mathematical Weekend Joint work with Romain
More informationOn the exponential map on Riemannian polyhedra by Monica Alice Aprodu. Abstract
Bull. Math. Soc. Sci. Math. Roumanie Tome 60 (108) No. 3, 2017, 233 238 On the exponential map on Riemannian polyhedra by Monica Alice Aprodu Abstract We prove that Riemannian polyhedra admit explicit
More informationA Brief Introduction to Mathematical Relativity
A Brief Introduction to Mathematical Relativity Arick Shao Imperial College London Arick Shao (Imperial College London) Mathematical Relativity 1 / 31 Special Relativity Postulates and Definitions Einstein
More informationPAPER 309 GENERAL RELATIVITY
MATHEMATICAL TRIPOS Part III Monday, 30 May, 2016 9:00 am to 12:00 pm PAPER 309 GENERAL RELATIVITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
More informationarxiv: v2 [math.dg] 26 Feb 2017
LOCAL AND GLOBAL BOUNDARY RIGIDITY AND THE GEODESIC X-RAY TRANSFORM IN THE NORMAL GAUGE PLAMEN STEFANOV, GUNTHER UHLMANN AND ANDRÁS VASY arxiv:170203638v2 [mathdg] 26 Feb 2017 Abstract In this paper we
More informationTravel Time Tomography and Tensor Tomography, I
Travel Time Tomography and Tensor Tomography, I Plamen Stefanov Purdue University Mini Course, MSRI 2009 Plamen Stefanov (Purdue University ) Travel Time Tomography and Tensor Tomography, I 1 / 17 Alternative
More informationTwo dimensional manifolds
Two dimensional manifolds We are given a real two-dimensional manifold, M. A point of M is denoted X and local coordinates are X (x, y) R 2. If we use different local coordinates, (x, y ), we have x f(x,
More informationRigidity and Non-rigidity Results on the Sphere
Rigidity and Non-rigidity Results on the Sphere Fengbo Hang Xiaodong Wang Department of Mathematics Michigan State University Oct., 00 1 Introduction It is a simple consequence of the maximum principle
More informationSyllabus. May 3, Special relativity 1. 2 Differential geometry 3
Syllabus May 3, 2017 Contents 1 Special relativity 1 2 Differential geometry 3 3 General Relativity 13 3.1 Physical Principles.......................................... 13 3.2 Einstein s Equation..........................................
More informationAsymptotic Behavior of Marginally Trapped Tubes
Asymptotic Behavior of Marginally Trapped Tubes Catherine Williams January 29, 2009 Preliminaries general relativity General relativity says that spacetime is described by a Lorentzian 4-manifold (M, g)
More informationQuasi-local Mass in General Relativity
Quasi-local Mass in General Relativity Shing-Tung Yau Harvard University For the 60th birthday of Gary Horowtiz U. C. Santa Barbara, May. 1, 2015 This talk is based on joint work with Po-Ning Chen and
More informationLOCAL EXISTENCE THEORY OF THE VACUUM EINSTEIN EQUATIONS
LOCAL EXISTENCE THEORY OF THE VACUUM EINSTEIN EQUATIONS TSZ KIU AARON CHOW Abstract. In this paper, we treat the vacuum Einstein equations as a local Cauchy problem. We show that it reduces to a system
More informationNumerical Methods for geodesic X-ray transforms and applications to open theoretical questions
Numerical Methods for geodesic X-ray transforms and applications to open theoretical questions François Monard Department of Mathematics, University of Washington. Nov. 13, 2014 UW Numerical Analysis Research
More informationThe generally covariant locality principle and an extension to gauge symmetries
The generally covariant locality principle and an extension to gauge symmetries Jochen Zahn Universität Wien Seminar Mathematische Physik, November 2012 Outline We review the framework of locally covariant
More informationPhoto-Acoustic imaging in layered media
Photo-Acoustic imaging in layered media Faouzi TRIKI Université Grenoble-Alpes, France (joint works with Kui Ren, Texas at Austin, USA.) 4 juillet 2016 Photo-Acoustic effect Photo-Acoustic effect (Bell,
More informationFefferman Graham Expansions for Asymptotically Schroedinger Space-Times
Copenhagen, February 20, 2012 p. 1/16 Fefferman Graham Expansions for Asymptotically Schroedinger Space-Times Jelle Hartong Niels Bohr Institute Nordic String Theory Meeting Copenhagen, February 20, 2012
More information