Bubbles Bad; Ripples Good

01/11/2011

Continuity of the infimum

Just realised (two seeks ago, but only gotten around to finish this blog posting now) that an argument used to prove a proposition in a project I am working on is wrong. After reducing the problem to its core I found that it is something quite elementary. So today’s post would be of a different flavour from the ones of recent past.

Question Let X,Y be topological spaces. Let f:X\times Y\to\mathbb{R} be a bounded, continuous function. Is the function g(x) = \inf_{y\in Y}f(x,y) continuous?

Intuitively, one may be tempted to say “yes”. Indeed, there are plenty of examples where the answer is in the positive. The simplest one is when we can replace the infimum with the minimum:

Example Let the space Y be a finite set with the discrete topology. Then g(x) = \min_{y\in Y} f(x,y) is continuous.
Proof left as exercise.

But in fact, the answer to the question is “No”. Here’s a counterexample:

Example Let X = Y = \mathbb{R} with the standard topology. Define

\displaystyle f(x,y) = \begin{cases} 1 & x > 0 \\ 0 & x < -e^{y} \\ 1 + x e^{-y} & x\in [-e^{y},0]  \end{cases}

which is clearly continuous. But the infimum function g(x) is roughly the Heaviside function: g(x) = 1 if x \geq 0, and g(x) = 0 if x < 0.

So what is it about the first example that makes the argument work? What is the different between the minimum and the infimum? A naive guess maybe that in the finite case, we are taking a minimum, and therefore the infimum is attained. This guess is not unreasonable: there are a lot of arguments in analysis where when the infimum can be assumed to be attained, the problem becomes a lot easier (when we are then allowed to deal with a minimizer instead of a minimizing sequence). But sadly that is not (entirely) the case here: for every x_0, we can certainly find a y_0 such that f(x_0,y_0) = g(x_0). So attaining the infimum point-wise is not enough.

What we need, here, is compactness. In fact, we have the following

Theorem If X,Y are topological spaces and Y is compact. Then for any continuous f:X\times Y\to\mathbb{R}, the function g(x) := \inf_{y\in Y} f(x,y) is well-defined and continuous.

Proof usually proceeds in three parts. That g(x) > -\infty follows from the fact that for any fixed x\in X, f(x,\cdot):Y\to\mathbb{R} is a continuous function defined on a compact space, and hence is bounded (in fact the infimum is attained). Then using that the sets (-\infty,a) and (b,\infty) form a subbase for the topology of \mathbb{R}, it suffices to check that g^{-1}((-\infty,a)) and g^{-1}((b,\infty)) are open.

Let \pi_X be the canonical projection \pi_X:X\times Y\to X, which we recall is continuous and open. It is easy to see that g^{-1}((-\infty,a)) = \pi_X \circ f^{-1}((-\infty,a)). So continuity of f implies that this set is open. (Note that this part does not depend on compactness of Y. In fact, a minor modification of this proof shows that for any family of upper semicontinuous functions \{f_c\}_C, the pointwise infimum \inf_{c\in C} f_c is also upper semicontinuous, a fact that is very useful in convex analysis. And indeed, the counterexample function given above is upper semicontinuous.)

It is in this last part, showing that g^{-1}((b,\infty)) is open, that compactness is crucially used. Observe that g(x) > b \implies f(x,y) > b~ \forall y. In other words g(x) > b \implies \forall y, (x,y) \in f^{-1}((b,\infty)) an open set. This in particular implies that \forall x\in g^{-1}((b,\infty)) \forall y\in Y there exists a “box” neighborhood U_{(x,y)}\times V_{(x,y)} contained in f^{-1}((b,\infty)). Now using compactness of Y, a finite subset \{(x,y_i)\} of all these boxes cover \{x\}\times Y. And in particular we have

\displaystyle \{x\}\times Y \subset \left(\cap_{i = 1}^k U_{(x,y_i)}\right)\times Y \subset f^{-1}((b,\infty))

and hence g^{-1}((b,\infty)) = \cup_{x\in g^{-1}((b,\infty))} \cap_{i = 1}^{k(x)} U_{x,y_i} is open. Q.E.D.

One question we may ask is how sharp is the requirement that Y is compact. As with most things in topology, counterexamples abound.

Example Let Y be any uncountably infinite set equipped with the co-countable topology. That is, the collection of open subsets are precisely the empty set and all subsets whose complement is countable. The two interesting properties of this topology are (a) Y is not compact and (b) Y is hyperconnected. (a) is easy to see: let C be some countably infinite subset of Y. For each c\in C let U_c = \{c\}\cup (Y\setminus C). This forms an open cover with not finite sub-cover. Hyperconnected spaces are, roughly speaking, spaces in which all open nonempty sets are “large”, in the sense that they mutually overlap a lot. In particular, a continuous map from a hyperconnected space to a Hausdorff space must be constant. In our case we can see this directly: suppose h:Y\to \mathbb{R} is a continuous map. Fix y_1,y_2\in Y. Let N_{1,2}\subset \mathbb{R} be open neighborhoods of f(y_{1,2}). Since h is continuous, h^{-1}(N_1)\cap h^{-1}(N_2) is open and non-empty (by the co-countable assumption). Therefore N_1\cap N_2\neq \emptyset for any pairs of neighborhoods. Since \mathbb{R} is Hausdorff, this forces h to be the constant map. This implies that for any topological space X, a continuous function f:X\times Y\to\mathbb{R} is constant along Y, and hence for any y_0\in Y, we have \inf_{y\in Y} f(x,y) =: g(x) = f(x,y_0) is continuous.

One can try to introduce various regularity/separation assumptions on the spaces X,Y to see at what level compactness becomes a crucial requirement. As an analyst, however, I really only care about topological manifolds. In which case the second counterexample up top can be readily used. We can slightly weaken the assumptions and still prove the following partial converse in essentially the same way.

Theorem Let X be Tychonoff, connected, and first countable, such that X contains a non-trivial open subset whose closure is not the entire space; and let Y be paracompact, Lindelof. Then if Y is noncompact, there exists a continuous function f:X\times Y\to\mathbb{R} such that \inf_{y\in Y}f:X\to \mathbb{R} is not continuous.

Remark Connected (nontrivial) topological manifolds automatically satisfy the conditions on X and Y except for non-compactness. The conditions given are not necessary for the theorem to hold; but they more or less capture the topological properties used in the construction of the second counterexample above.

Remark If X is such that every open set’s closure is the entire space, we must have that it is hyperconnected (let C\subset X be a closed set. Suppose D\subset X is another closed set such that C\cup D = X. Then C\subset D^c and vice versa, but D^c is open, so C = X. Hence X cannot be written as the union of two proper closed subsets). And if it is Tychonoff, then X is either the empty-set or the one-point set.

Lemma For a paracompact Lindelof space that is noncompact, there exists a countably infinite open cover \{U_k\} and a sequence of points y_k \in U_k such that \{y_k\}\cap U_j = \emptyset if j\neq k.

Proof: By noncompactness, there exists an open cover that is infinite. By Lindelof, this open cover can be assumed to be countable, which we enumerate by \{V_k\} and assume WLOG that \forall k, V_k \setminus \cup_{j =1}^{k-1} V_j \neq \emptyset. Define \{U_k\} and \{y_k\} inductively by: U_k = V_k \setminus \cup_{j = 1}^{k-1} \{ y_j\} and choose y_k \in U_k \setminus \cup_{j=1}^{k-1}U_j.

Proof of theorem: We first construct a sequence of continuous functions on X. Let G\subset X be a non-empty open set such that its closure-complement H = (\bar{G})^c is a non-empty open set (G exists by assumption). By connectedness \bar{G}\cap \bar{H} \neq \emptyset, so we can pick x_0 in the intersection. Let \{x_j\}\subset H be a sequence of points converging to x_0, which exists by first countability. Using Tychonoff, we can get a sequence of continuous functions f_jon X such that f_j|_{\bar{G}} = 0 and f_j(x_j) = -1.

On Y, choose an open cover \{U_k\} and points \{y_k\} per the previous Lemma. By paracompactness we have a partition of unity \{\psi_k\} subordinate to U_k, and by the conclusion of the Lemma we have that \psi_k(y_k) = 1. Now we define the function

\displaystyle f(x,y) = \sum_{k} f_k(x)\psi_k(y)

which is continuous, and such that f|_{\bar{G}\times Y} = 0. But by construction \inf_{y\in Y}f(x,y) \leq f(x_k,y_k) = f_k(x_k) = -1, which combined with the fact that x_k \to x_0 \in \bar{G} shows the desired result. q.e.d.

30/09/2011

Gauge invariance, geometrically

A somewhat convoluted chain of events led me to think about the geometric description of partial differential equations. And a question I asked myself this morning was

Question
What is the meaning of gauge invariance in the jet-bundle treatment of partial differential equations?

The answer, actually, is quite simple.

Review of geometric formulation PDE
We consider here abstract PDEs formulated geometrically. All objects considered will be smooth. For more about the formal framework presented here, a good reference is H. Goldschmidt, “Integrability criteria for systems of nonlinear partial differential equations”, JDG (1967) 1:269–307.

A quick review: the background manifold X is assumed (here we take a slightly more restrictive point of view) to be a connected smooth manifold. The configuration space \mathcal{C} is defined to be a fibred manifold p:\mathcal{C}\to X. By J^r\mathcal{C} we refer to the fibred manifold of r-jets of \mathcal{C}, whose projection p^r = \pi^r_0 \circ p where for r > s we use \pi^r_s: J^r\mathcal{C}\to J^s\mathcal{C} for the canonical projection.

A field is a (smooth) section \phi \subset \Gamma \mathcal{C}. A simple example that capture most of the usual cases: if we are studying mappings between manifolds \phi: X\to N, then we take \mathcal{C} = N\times X the trivial fibre bundle. The s-jet operator naturally sends j^s: \Gamma\mathcal{C} \ni \phi \mapsto j^s\phi \in \Gamma J^r\mathcal{C}.

A partial differential equation of order r is defined to be a fibred submanifold J^r\mathcal{C} \supset R^r \to X. A field is said to solve the PDE if j^r\phi \subset R^r.

In the usual case of systems of PDEs on Euclidean space, X is taken to be \mathbb{R}^d and \mathcal{C} = \mathbb{R}^n\times X the trivial vector bundle. A system of m PDEs of order r is usually taken to be F(x,\phi, \partial\phi, \partial^2\phi, \ldots, \partial^r\phi) = 0 where

\displaystyle F: X\times \mathbb{R}^n \times \mathbb{R}^{dn} \times \mathbb{R}^{\frac{1}{2}d(d+1)n} \times \cdots \times \mathbb{R}^{{d+r-1 \choose r} n} \to \mathbb{R}^m

is some function. We note that the domain of F can be identified in this case with J^r\mathcal{C}, We can then extend F to \tilde{F}: J^r\mathcal{C} \ni c \mapsto (F(c),p^r(c)) \in \mathbb{R}^m\times X a fibre bundle morphism.

If we assume that \tilde{F} has constant rank, then \tilde{F}^{-1}(0) is a fibred submanifold of J^r\mathcal{C}, and this is our differential equation.

Gauge invariance
In this frame work, the gauge invariance of a partial differential equation relative to certain symmetry groups can be captured by requiring R^r be an invariant submanifold.

More precisely, we take

Definition
A symmetry/gauge group \mathcal{G} is a subgroup of \mathrm{Diff}(\mathcal{C}), with the property that for any g\in\mathcal{G}, there exists a g'\in \mathrm{Diff}(X) with p\circ g = g' \circ p.

It is important we are looking at the diffeomorphism group for \mathcal{C}, not J^r\mathcal{C}. In general diffeomorphisms of J^r\mathcal{C} will not preserve holonomy for sections of the form j^r\phi, a condition that is essential for solving PDEs. The condition that the symmetry operation “commutes with projections” is to ensure that g:\Gamma\mathcal{C}\to\Gamma\mathcal{C}, which in particular guarantees that g extends to a diffeomorphism of J^rC with itself that commutes with projections.

From this point of view, a (system of) partial differential equation(s) R^r is said to be \mathcal{G}-invariant if for every g\in\mathcal{G}, we have g(R^r) \subset R^r.

We give two examples showing that this description agrees with the classical notions.

Gauge theory. In classical gauged theories, the configuration space \mathcal{C} is a fibre bundle with structure group G which acts on the fibres. A section of G\times X \to X induces a diffeomorphism of \mathcal{C} by fibre-wise action. In fact, the gauge symmetry is a fibre bundle morphism (fixes the base points).

General relativity. In general relativity, the configuration space is the space of Lorentzian metrics. So the background manifold is the space-time X. And the configuration space is the open submanifold of S^2T^*X given by non-degenerate symmetric bilinear forms with signature (-+++). A diffeomorphism \Psi:X\to X induces T^*\Psi = (\Psi^{-1})^*: T^*X \to T^*X and hence a configuration space diffeomorphism that commutes with projection. It is in this sense that Einstein’s equations are diffeomorphism invariant.

Notice of course, this formulation does not contain the “physical” distinction between global and local gauge transformations. For example, for a linear PDE (so \mathcal{C} is a vector bundle and R^r is closed under linear operations), the trivial “global scaling” of a solution is considered in this frame work a gauge symmetry, though it is generally ignored in physics.

28/09/2011

Moving SVN servers

Filed under: Life of a mathematician — Willie Wong @ 10:57

This is mostly to document the process for my own benefit, so the information is easier to find the next time I need to do it. This is also a follow-up to this previous post of mine.

Suppose the repository on the old server is called “WorkSVN” as in my previous post. First we need to dump the contents from the old server by

svnadmin dump /home/XXYY/WorkSVN > WorkSVN.dmp

We then copy the file over to the new server. On the new server issue

svnadmin create ~/WorkSVN

as before, and followed by

svnadmin load ~/WorkSVN< WorkSVN.dmp

to load the database dump.

Then create the working copy as usual

svn co file:///home/XXYY/WorkSVN/PaperDrafts ~/PaperDrafts

For the client computers, update the server path via

svn switch –relocate OLDURL NEWURL

You can find the oldurl by reading “svn info” of the working copy.

30/08/2011

Extensions of (co)vector fields to tangent bundles

I am reading Sasaki’s original paper on the construction of the Sasaki metric (a canonical Riemannian metric on the tangent bundle of a Riemannian manifold), and the following took me way too long to understand. So I’ll write it down in case I forgot in the future.

In section two of the paper, Sasaki consider “extended transformations and extended tensors”. Basically he wanted to give a way to “lift” tensor fields from a manifold to tensor fields of the same rank on its tangent bundle. And he did so in the language of coordinate changes, which geometrical content is a bit hard to parse. I’ll discuss his construction in a bit. But first I’ll talk about something different.

The trivial lifts
Let M, N be smooth manifolds, and let f:M\to N a submersion. Then we can trivially lift covariant objects on N to equivalent objects on M by the pull-back operation. To define the pull-back, we start with a covariant tensor field \vartheta \in \Gamma T^0_kN, and set f^*\vartheta \in \Gamma T^0_kM by the formula:

\displaystyle f^*\vartheta(X_1,\ldots,X_k) = \vartheta(df\circ X_1, \ldots, df\circ X_k)

where the X_1, \ldots, X_k \in T_pM, and we use that df(p): T_pM \to T_{f(p)}N. Observe that for a function g: N \to \mathbb{R}, the pull-back is simply f^*g = g\circ f :M\to N\to\mathbb{R}.

On the other hand, for contravariant tensor fields, the pull-back is not uniquely defined: using that f is a submersion, we have that TM / \ker(df) = TN, so while, given a vector field v on N, we can always find a vector field w on M such that df(w) = v, the vector field w is only unique up to an addition of a vector field that lies in the kernel of df. If, however, that M is Riemannian, then we can take the orthogonal decomposition of TM into the kernel and its complement, thereby getting a well-defined lift of the vector field (in other words, by exploiting the identification between the tangent and cotangent spaces).

Remarkably, the extensions defined by Sasaki is not this one.

(Let me just add a remark here: given two manifolds, once one obtain a well defined way of lifting vectors, covectors, and functions from one to the other, such that they are compatible (\vartheta^*(v^*) = [\vartheta(v)]^*), one can extend this mapping to arbitrary tensor fields.)

The extensions defined by Sasaki
As seen above, if we just rely on the canonical submersion \pi:TM\to M, we cannot generally extend vector fields. Sasaki’s construction, however, strongly exploits the fact that TM is the tangent bundle of M.

We start by looking at the vector field extension defined by equation (2.6) of the linked paper. We first observe that a vector field v on a manifold M is a section of the tangent bundle. That is, v is a map M\to TM such that the composition with the canonical projection \pi\circ v:M\to M is the identity map. This implies, using the chain rule, that the map d(\pi\circ v)= d\pi \circ dv: TM\to TM is also the identity map. Now, d\pi: T(TM) \to TM is the projection induced by the projection map \pi, which is different from the canonical projection \pi_2: T(TM) \to TM from the tangent bundle of a manifold to the manifold itself. However, a Proposition of Kobayashi (see “Theory of Connections” (1957), Proposition 1.4), shows that there exists an automorphism \alpha:T(TM) \to T(TM) such that d\pi \circ \alpha = \pi_2 and \pi_2\circ\alpha = d\pi. So v as a differential mapping induces a map \alpha\circ dv: TM \to T(TM), which is a map from the tangent bundle TM to the double tangent bundle T(TM), which when composed with the canonical projection \pi_2 is the identity. In other words, \alpha\circ dv is a vector field on TM.

Next we look at the definition (2.7) for one-forms. Give \vartheta a one-form on M, it induces naturally a scalar function on TM: for p\in M, v\in T_pM, we have \vartheta: TM\to \mathbb{R} taking value \vartheta(p)\cdot v. Hence its differential d\vartheta is a one-form over TM.

Now, what about scalar functions? Let \vartheta be a one-form and v be a vector field on M, we consider the pairing of their extensions to TM. It is not too hard to check that the corresponding scalar field to \vartheta(v), when evaluated at (p,w)\in TM, is in fact d(\vartheta(v))|_{p,w}, the derivative of the scalar function \vartheta(v) in the direction of w at point p. In general, the compatible lift of scalar fields g:M\to \mathbb{R} to TM is the function \tilde{g}(p,v) = dg(p)[v].

Using this we can extend the construction to arbitrary tensor fields, and a simple computation yields that this construction is in fact identical, for rank-2 tensors, to the expressions given in (2.8), (2.9), and (2.10) in the paper.

The second extension
The above extension is not the only map sending vectors on M to vectors on TM. In the statement of Lemmas 3 there is also another construction. Given a vector field v, it induces a one parameter family of diffeomorphisms on TM via that maps \psi_t(p,w) = (p, w+vt). Its differential \frac{d}{dt}\psi_t|_{t=0} is a vector field over TM.

The construction in the statement of Lemma 4 is the trivial one mentioned at the start of this post.

09/06/2011

Decay of waves IIIb: tails for homogeneous linear equation on curved background

Now we will actually show that the specific decay properties of the linear wave equation on Minkowski space–in particular the strong Huygens’ principle–is very strongly tied to the global geometry of that space-time. In particular, we’ll build, by hand, an example of a space-time where geometry itself induces back-scattering, and even linear, homogeneous waves will exhibit a tail.

For convenience, the space-time we construct will be spherically symmetric, and we will only consider spherically symmetric solutions of the wave equation on it. We will also focus on the 1+3 dimensional case. (more…)

16/05/2011

Decay of waves IIIa: nonlinear tails in Minkowski space redux

Before we move on to the geometric case, I want to flesh out the nonlinear case mentioned in the end of the last post a bit more. Recall that it was shown for generic nonlinear (actually semilinear; for quasilinear and worse equations we cannot use Duhamel’s principle) wave equations, if we put in compact support for the initial data, we expect the first iterate to exhibit a tail. One may ask whether it is possible that, in fact, this is an artifact of the successive approximation scheme; that in fact somehow it always transpires that a conspiracy happens, and all the higher order iterates cancel out the tail coming from the first iterate. This is rather unlikely, owing to the fact that the convergence to \phi_\infty is dominated by a geometric series. But to just make double sure, here we give a nonlinear system of wave equations such that the successive approximation scheme converges after finitely many steps (in fact, after the first iterate), and so we can also explicitly compute the rate of decay for the nonlinear tail. While the decay rate is not claimed to be generic (though it is), the existence of one such example with a fixed decay rate shows that for a statement quantifying over all nonlinear wave equations, it would be impossible to demonstrate better decay rate than the one exhibited. (more…)

14/05/2011

Decay of waves IIb: Minkowski space, with right-hand side

In the first half of this second part of the series, we considered solutions to the linear, homogeneous wave equation on flat Minkowski space, and showed that for compactly supported initial data, we have strong Huygens’ principle. We further made references to the fact that this behaviour is expected to be unstable. In this post, we will further illustrate this instability by looking at Equation 1 first with a fixed source F = F(t,x), and then with a nonlinearity F = F(t,x, \phi, \partial\phi).

Duhamel’s Principle

To study how one can incorporate inhomogeneous terms into a linear equation, and to get a qualitative grasp of how the source term contributes to the solution, we need to discuss the abstract method known as Duhamel’s Principle. We start by illustrating this for a very simple ordinary differential equation.

Consider the ODE satisfied by a scalar function \alpha:

Equation 13
\displaystyle \frac{d}{ds}\alpha(s) = k(s)\alpha(s) + \beta(s)

when \beta\equiv 0, we can easily solve the equation with integration factors

\displaystyle \alpha(s) = \alpha(0) e^{\int_0^s k(t) dt}

Using this as a sort of an ansatz, we can solve the inhomogeneous equation as follows. For convenience we denote by K(s) = \int_0^s k(t) dt the anti-derivative of k. Then multiplying Equation 13 through by \exp -K(s), we have that

Equation 14
\displaystyle \frac{d}{ds} \left( e^{-K(s)}\alpha(s)\right) = e^{-K(s)}\beta(s)

which we solve by integrating

Equation 15
\displaystyle \alpha(s) = e^{K(s)}\alpha(0) + e^{K(s)} \int_0^s e^{-K(t)}\beta(t) dt

If we write K(s;t) = \int_t^s k(u) du, then we can rewrite Equation 15 as given by an integral operator

Equation 15′
\displaystyle \alpha(s) = e^{K(s)}\alpha(0) + \int_0^s e^{K(s;t)}\beta(t) dt

(more…)

12/05/2011

Decay of waves IIa: Minkowski background, homogeneous case

Now let us get into the mathematics. The wave equations that we will consider take the form

Equation 1
-\partial_t^2 \phi + \triangle \phi = F

where \phi:\mathbb{R}^{1+n}\to\mathbb{R} is a real valued function defined on (1+n)-dimensional Minkowski space that describes our solution, and F represents a “source” term. When F vanishes identically, we say that we are looking at the linear, homogeneous wave equation. When F is itself a function of \phi and its first derivatives, we say that the equation is a semilinear wave equation.

We first start with the homogeneous, linear case.

Homogeneous wave equation in one spatial dimension

One interesting aspect of the wave equation is that it only possesses the second, multidimensional, dispersive mechanism as described in my previous post. In physical parlance, the “phase velocity” and the “group velocity” of the wave equation are the same. And therefore, a solution of the wave equation, quite unlike a solution of the Schroedinger equation, will not exhibit decay when there is only one spatial dimension (mathematically this is one significant difference between relativistic and quantum mechanics). In this section we make a computation to demonstrate this, a fact that would also be useful later on when we look at higher (in particular, three) dimensions.

Use x\in\mathbb{R} for the variable representing spatial position. The wave equation can be written as

-\partial_t^2 \phi + \partial_x^2\phi = 0

Now we perform a change of variables: let u = \frac{1}{2}(t-x) and v = \frac{1}{2}(t+x) be the canonical null variables. The change of variable formula replaces

Equation 2
\displaystyle \partial_t \to \frac{\partial u}{\partial t} \partial_u + \frac{\partial v}{\partial t} \partial v = \frac{1}{2}\partial_u + \frac{1}{2}\partial_v
\displaystyle \partial_x \to \frac{\partial u}{\partial x} \partial_u + \frac{\partial v}{\partial x} \partial v = -\frac{1}{2}\partial_u + \frac{1}{2}\partial_v

and we get that in the (u,v) coordinate system,

Equation 3
-\partial_u \partial_v \phi = 0

(more…)

10/05/2011

Getting Google Buzz to work

Filed under: Uncategorized — Willie Wong @ 12:59
Tags:

Allegedly one can link a blog to Google Buzz. And I am pretty sure it works, since Terry’s Buzz gets also his What’s New feed. But for the life of me I cannot get Google to realise that this blog is affiliated to my google profile.

And I hate to admit defeat to technology.

05/05/2011

Decay of waves I: Introduction

In the next week or so, I will compose a series of posts on the heuristics for the decay of the solutions of the wave equation on curved (and flat) backgrounds. (I have my fingers crossed that this does not end up aborted like my series of posts on compactness.) In this first post I will give some physical intuition of why waves decay. In the next post I will write about the case of linear and nonlinear waves on flat space-time, which will be used to motivate the construction, in post number three, of an example space-time which gives an upper bound on the best decay that can be generally expected for linear waves on non-flat backgrounds. This last argument, due to Mihalis Dafermos, shows that why the heuristics known as Price’s Law is as good as one can reasonably hope for in the linear case. (In the nonlinear case, things immediately get much much worse as we will see already in the next post.)

This first post will not be too heavily mathematical, indeed, the only realy foray into mathematics will be in the appendix; the next ones, however, requires some basic familiarity with partial differential equations and pseudo-Riemannian geometry. (more…)

Next Page »

Theme: Rubric. Blog at WordPress.com.

Follow

Get every new post delivered to your Inbox.