Loading [MathJax]/jax/output/HTML-CSS/jax.js

Sunday, January 14, 2018

New Location

I have a new blog at github pages. All of the content of this blog has been migrated there. This blog is no longer being maintained.

Tuesday, September 15, 2015

Santalo Formula

Let M be a simple Riemannian manifold with boundary M. For (x,v)SM, let τ(x,v) denote the exit time of the geodesic starting at x with tangent vector v, i.e. τ(x,v) is the (necessarily unique, and finite) time at which expx(tv)M.

We let +(SM) denote the set
+(SM)={(x,v)SM | xM,v,ν>0}
where ν denotes the inward unit normal to M in M. The exponential map identifies SM with the set
Ω={(x,v,t)+(SM)×R | 0tτ(x,v)},
via (x,v,t)expx(tv). Let Φ:ΩM denote this diffeomorphism. Then we have, for all fC(SM)
SMfdvol(SM)=Ω(Φf)(Φdvol(SM))=+(SM)τ(x,v)0f(ϕt(x,v))Φdvol(SM).
Therefore, we can compute integrals of functions over SM by integrating along geodesics, provided that we can cmopute Φdvol(SM). This is the content of the Santalo formula.

Theorem (Santalo formula). For all fC(SM), we have
SMfdvol(SM)=+(SM)τ(x,v)0f(ϕt(x,v))v,νdtdvol((SM))

Proof. Necessarily, we must have
Φ(dvol(SM))=a(x,v)dtdvol((SM))),
for some function a(x,v). The reason we can assume that a is independent of t is that Φ is defined via geodesic flow, and geodesic flow preserves the volume form on SM. To compute the factor a(x,v), we just need to compute
i/tΦ(dvol(SM))=Φ(iΦ(/t)dvol(SM))
From the definition of Φ, we have that Φ(/t) is the Reeb vector field on SM, i.e. the vector field generating geodesic flow. Therefore, Φ(/t) is equal, at a point (x,v) to the horizontal lift of the vector v. Therefore, using the definition of the induced volume form on a hypersurface of a Riemannian manifold, we find
i/tΦ(dvol(SM))=v,νdvol((SM))
where ν is the inward pointing unit normal to (SM) in SM. This shows that a(x,v)=v,ν and completes the proof.

Friday, September 11, 2015

The Index Form

Let f:[0,T]×(ϵ,ϵ)M be a family of parametrized curves in a Riemannian manifold (M,g). To simplify this calculation, we assume that f(0,s)=p,f(T,s)=q for some p,qM and all s(ϵ,ϵ). (This assumption is not necessary, but without it our variational formulae will have additional boundary terms.)

For convenience, set ˙f=f/t and f=f/s. For each s(ϵ,ϵ) we define the energy functional E=E(s) to be
E(s)=12T0|˙f|2dt.
The first variation is
dEds=T0f˙f,˙fdt =T0˙ff,˙fdt =T0f,˙f˙fdt

Set γ(t):=f(t,0) and X(t)=f(t) (thought of as a vector field supported on γ). Evaluating the above at s=0 we obtain
dEds|s=0=T0X,˙γ˙γdt,
which shows immediately that

Theorem. γ is a critical point of the energy functional if and only if ˙γ˙γ=0.


The second variation is
d2Eds2=T0ff,˙f˙f+f,f˙f˙fdt =T0ff,˙f˙f+f,˙ff˙fdt+f,R(f,˙f)˙fdt =T0ff,˙f˙f˙ff,f˙fdt+f,R(f,˙f)˙fdt =T0ff,˙f˙f˙ff,˙ffdt+f,R(f,˙f)˙fdt

Assume now that γ is a geodesic, i.e. ˙γ˙γ=0. Then evaluating the above at s=0, we obtain
d2Eds2=T0|˙γX|2X,R(X,˙γ)˙γdt.

Definition. Let γ be a geodesic. The index form associated to variations X,Y of γ is
I(X,Y)=T0˙γX,˙γYdtY,R(X,˙γ)˙γ =T0Y,2˙γX+R(X,˙γ)˙γ
It follows from symmetries of the Riemann tensor that I(X,Y)=I(Y,X) and also I(X,X)=E as above.

Theorem. Suppose that X is the infinitesimal variation of a family of affine geodesics about a fixed geodesic γ. Then
2˙γX+R(X,˙γ)˙γ=0.
In particular, I(X,)=0.

Proof. Let f(t,s) denote the family as above. By hypothesis, we have that ˙f˙f=0 for all s, so that
f˙f˙f=0.
Commuting the derivatives using the curvature tensor, we have
0=˙ff˙f+R(f,˙f)˙f.
Now use ˙ff=f˙f and evaluate at s=0 to obtain
0=2˙γX+R(X,˙γ)˙γ.

Thursday, September 3, 2015

Boundary Distance

Recently, I've been learning some topics related to machine learning, and especially manifold learning. These both fall under the general notion of inverse problems: given some mathematical object X (it could be a function f:AB, or a Riemannian manifold (M,g), or a probability measure dμ on a space X, etc.), can we effectively reconstruct X given only the information of some auxiliary measurements? What if we can only perform finitely many measurements? What if the measurements are noisy? Can we reconstruct X at least approximately? Can we measure in some precise way, how close our approximate reconstruction is to the unknown object X? And so on, and so forth.

Anyway, this post is about a cute observation, which I was reminded of while reading a paper on the inverse Gel'fand problem. Let M be a compact manifold with smooth boundary M. Then with no additional data required, we have a Banach space L(M) consisting of the essentially bounded measureable functions on the boundary. Since it is a Banach space, it comes with a complete metric d(f,g):=fgL(M).

Now, suppose that g is a Riemannian metric on M. Then we have the Riemannian distance function dg(x,y) which is defined to be the infimum of arclengths of all smooth paths connecting x and y. For any xM, we obtain a function rxL(M) defined by
rx(z)=dg(x,z),zM.
This gives a map ϕg:ML(M), defined by xrx.

Theorem. Suppose that for any two distinct x,yM, there is a unique length-minimizing geodesic connecting x and y. Then ϕg:ML(M) is an isometric embedding, i.e. dg(x,y)=d(rx,ry) for all x,yM.

Proof. Let x,y be distinct and let γ be the unique geodesic from x to y. For any point z on the boundary, we have
|dg(x,z)dg(y,z)|dg(x,y).
which is the triangle inequality. Now let γ be the unique geodesic from x to y, and extend γ until it hits some boundary point z. Then since x,y,z all lie on a length-minimizing geodesic, we have
dg(x,z)dg(y,z)=dg(x,y).
Therefore, the bound above is always saturated, and we find
supzM|dg(x,z)dg(y,z)|=dg(x,y).
But the expression on the left is nothing but the L(M)-norm of rxry, so the theorem is proved.

Monday, August 31, 2015

Hamilton-Jacobi equation and Riemannian distance

Consider the cotangent bundle TX as a symplectic manifold with canonical symplectic form ω. Consider the Hamilton-Jacobi equation
St+H(x,S)=0,
for the classical Hamilton function S(x,t). Setting x=x(t),p(t)=(S)(x(t),t) one sees immediately from the method of characteristics that this PDE is solved by the classical action
S(x,t)=t0(p˙xH)ds,
where the integral is taken over the solution (x(s),p(s)) of Hamilton's equations with x(0)=x0 and x(t)=x. The choice of basepoint x0 involves an overall additive constant of S, and really this solution is only valid in some neighbourhood U of x0. (Reason: S is in general multivalued, as the differential "dS" is closed but not necessarily exact.)

Now consider the case where X is Riemannian, with Hamiltonian H(x,p)=12|p|2. The solutions to Hamilton's equations are affinely parametrized geodesics, and by a simple Legendre transform we have
S(x,t)=12t0|˙x|2ds
where the integral is along the affine geodesic with x(0)=x0 and x(t)=x. Since x(s) is a geodesic, |˙x(s)| is a constant (in s) and therefore
S(x,t)=t2|˙x(0)|2.
Now consider the path γ(s)=x(|\dot x(0)|^{-1}s). This is an affine geodesic with γ(0)=x0, γ(|˙x(0)|t)=x and |˙γ|=1. Therefore, the Riemannian distance between x0 and x (provided x is sufficiently close to x0) is
d(x0,x)=|˙x(0)|t.
Combining this with the previous calculation, we see that
S(x,t)=12td(x0,x)2.
Now insert this back into the Hamilton-Jacobi equation above. With a bit of rearranging, we have the following.

Theorem. Let x0 denote a fixed basepoint of X. Then for all x in a sufficiently small neighborhood U of x0, the Riemannian distance function satisfies the Eikonal equation
|xd(x0,x)|2=1.

Now, for convenience set r(x)=d(x0,x). Then |r|2=1, from which we obtain (by differentiating twice and contracting)
gijgkl(lkirjr+kirljr)=0.
Quick calculation shows that
lkir=ilkrRlibkbr
Therefore, tracing over l and k we obtain
glklkir=i(Δr)+Rc(r,)
Plugging this back into the equation derived above, we have
r(Δr)+Rc(r,r)+|Hr|2=0,
where Hr denotes the Hessian of r regarded as a 2-tensor. Now, using r as a local coordinate, it is easy to see that r=r (as vector fields). So we can rewrite this identity as
r(Δr)+Rc(r,r)+|Hr|2=0.

Now, we can get a nice result out of this. First, note that the Hessian Hr always has at least one eigenvalue equal to zero, because the Eikonal equation implies that Hr(r,)=0. Let λ2,,λn denote the non-zero eigenvalues of Hr. We have
|Hr|2=λ22++λ2n,
while on the other hand
|Δr|2=(λ2++λn)2
By Cauchy-Schwarz, we have
|Δr|2(n1)|Hr|2

Proposition. Suppose that the Ricci curvature of X satisfies Rc(n1)κ, and let u=(n1)(Δr)1. Then
u1+κu2.

Proof. From preceding formulas, |Hr|2 can be expressed in terms of the Ricci curvature and the radial derivative of Δr. On the other hand, |Δ|2 is bounded above by (n1)|Hr|2. The claimed inequality then follows from simple rearrangement.

Now, the amazing thing is that this deceptively simple inequality is the main ingredient of the Bishop-Gromov comparison theorem. The Bishop-Gromov comparison theorem, in turn, is the main ingredient of the proof of Gromov(-Cheeger) precompactness. I hope to discuss these topics in a future post.

Tuesday, August 18, 2015

The Classical Partition Function

Let (M,ω) be a symplectic manifold of dimension 2n, and let H:MR be a classical Hamiltonian. The symplectic form ω allows us to define a measure on M, given by integration against the top form ωn/n!. We will denote this measure by dμ.

We imagine that (M,ω,H) represents some classical mechanical system. We suppose that the dynamics of this dynamical system are very complicated, e.g. some system of 1023 particles. The system is so complicated that not only can we not solve the equations of motion exactly, and even if we could, their solutions might be so complicated that we can't expect to learn very much from them.

So instead, we ask statistical questions. Imagine that we cannot measure the state of the system exactly (e.g. particles in a box), so we try to guess a probability distribution ρ(x,p,t) on M indicating that at time t the system has probability ρ(x,p,t)dμ of being in the state (x,p). Obviously, ρ should satisfy the constraint Mρdμ=1.

How does ρ evolve in time? We know that the system obeys Hamilton's equations,
(˙x,˙p)=XH=(H/p,H/x)
 in local Darboux coordinates. Therefore, a particle located at (x,p) in phase space at time t will be located at (x,p)+XHdt in phase space at time t+dt. Therefore, the probability that a particle is at point (x,p) at time t+dt, should be equal to the probability that the particle is at point (x,p)XHdt at time t. Therefore, we have
ρt=HxρpHpρx={H,ρ}

Given a probability distribution ρ, the entropy is defined to be
S[ρ]=Mρlogρdμ.

(A version of) the second law of thermodynamics. For a given average energy U, the system assumes a distribution of maximal possible entropy at thermodynamic equilibrium.

The goal now, is to determine what distribution ρ will maximize the entropy, subject to the constraints (for fixed U)
MHρdμ=U Mρdμ=1

Setting aside technical issues of convergence, etc., this variational problem is easily solved using the method of Lagrange multipliers. Introducing parameters λ1,λ2, we consider the modified functional
S[ρ,λ1,λ2,U]=M(ρlogρ+λ1ρ+λ2(Hρ))dμλ1λ2U.

Note that S/U=λ2, and this is conventionally identified with (minus) the inverse temperature.

Taking the variation with respect to ρ, we find
0=δSδρ=logρ1+λ1+Hλ2
Therefore, rho is proportional to eβH where we have set β=λ2. Define the partition function Z to be
Z=MeβHdμ.
We therefore have proved (formally and heuristically only!):

Theorem. The probability distribution ρ assumed by the system at thermodynamic equilibrium is given by
ρ=eβHZ
where β>0 is a real parameter, called the inverse temperature.

Corollary. At thermodynamic equilibrium, the average energy is given by
U=logZβ,
and the entropy is given by
S=βU+logZ.

Thursday, January 29, 2015

What is generalized geometry?

The following are my notes for a short introductory talk. References below are not intended to be comprehensive!

Math references:

 
Physics references:

 What is geometry?

Before trying to define generalized geometry, we should first decide what we mean by ordinary geometry. Of course, this question doesn't have a unique answer, so there are many ways to generalize the classical notions of manifolds and varieties. The viewpoint taken in generalized geometry is the following: the distinguishing feature of smooth manifolds is the existence of a tangent bundle
TMM
which satisfies some nice axioms. The basic idea of generalized geometry is to replace the tangent bundle with some other vector bundle LM, again satisfying some nice axioms. Different generalized geometries on M will correspond to different choices of bundle LM, as well as auxiliary data compatible with L in some appropriate sense.

Definition. A Lie algebroid over M is a smooth vector bundle LM together with a vector bundle map a:LTM called the anchor map and a bracket [,]:H0(M,L)H0(M,L)H0(M,L) satisfying the following axioms:
  • [,] is  a Lie bracket on H0(M,L)
  • [X,fY]=f[X,Y]+a(X)fY for X,YH0(M,L) and fH0(M,OM)
Note that we can take L to be either a real or complex vector bundle. In the latter case the anchor map should map to the complexified tangent bundle.

Example 1. We can take L to be TM with anchor map the identity.

Example 2. Let σ be a Poisson tensor on M. Then we define a bracket by [X,Y]=σ(X,Y) and an anchor by Xσ(X,). This makes TM into a Lie algebroid.

Example 3. Let M be a complex manifold of and let LTMC be the sub-bundle of vectors spanned by {/z1,,/zn} in local holomorphic coordinates. Then LM is a (complex) Lie algebroid.


Courant Bracket

We'd like to try to fit the preceding examples into a common framework. Let TM=TMTM. This bundle has a natural symmetric bilinear pairing given by
Xα,Yβ=12α(Y)+12β(X)
Note that this bilinear form is of split signature (n,n). We define a bracket on sections of TM by
[Xα,Yβ]=[X,Y](LXβ+12(dα(Y))LYα12d(β(X)))
Note that this bracket is not a Lie bracket. We also have an anchor map a:TMTM which is just the projection.

 Let B be a 2-form on M. Define an action of B on sections of TM by
X+αX+α+iXB

Proposition. This action preserves the Courant bracket if and only if B is closed.

This shows that the diffeomorphisms of M as a generalized manifold are large than the ordinary diffeomorphisms of M. In fact is is the semidirect product of the diffeomorphism group of M with the vector space of closed 2-forms.

Dirac Structures

Definition. A Dirac structure on M is an Lagrangian sub-bundle LTM which is closed under the Courant bracket.

Theorem (Courant). A Lagrangian sub-bundle LTM is a Dirac structure if and only if LM is a Lie algebroid over M, with bracket induced by the Courant bracket and anchor given by projection.

Example 1. TMTM.

Example 2. Take L to be the graph of a Poisson tensor.

Example 3. Take L to be the graph of a closed 2-form.

Admissible Functions

We now let LM be a Dirac structure on M.

Definition. A smooth function f on M is called admissible if there exists a vector field Xf such that (Xf,df) is a section of L.

The Poisson bracket is defined as follows. If f,g are admissible, then define
{f,g}=Xfg.
It is easy to check from the definitions that the bracket on admissible functions is well-defined (independent of choice of Xf) and skew-symmetric. With a little bit of calculation, we find the following.

Proposition. The vector space of admissible functions is naturally a Poisson algebra, and moreover the natural bracket satisfies the Leibniz rule.


Generalized Complex Structures

Definition. A generalized complex structure is a skew endomorphism J of TM such that J2=1 and such that the +i-eigenbundle is involutive under the Courant bracket.

Equivalently: A generalized complex structure is a (complex) Dirac structure LTM satisfying the condition L¯L=0.

Example 1. Let J be an ordinary complex structure on M. Then the endomorphism
[J00J]
defines a generalized complex structure on M.

Example 2. Let ω be a symplectic form on M. Then the endomorphism
[0ω1ω0]
defines a generalized complex structure on M.

Thus, generalized geometry gives a common framework for both complex geometry and symplectic geometry. Such a connection is exactly what is conjectured by mirror symmetry.

Example 3. Let J be a complex structure on M and let σ be a holomorphic Poisson tensor. Consider the subbundle LTM defined as the span of
ˉz1,,ˉzn,dz1σ(dz1),,dznσ(dzn)
Then L defines a generalized complex structure on M.

The last example shows that deformations of M as a generalized  complex manifold contain non-commutative deformations of the structure sheaf.  We also have the following theorem, which shows that there is an intimate relation between generalized complex geometry and holomorphic Poisson geometry.

Theorem (Bailey). Near any point of a generalized complex manifold, M is locally isomorphic to the product of a holomorphic Poisson manifold with a symplectic manifold.


Generalized Kähler Manifolds

Let (g,J,ω) be a Kähler triple. The Kähler property requires that
ω=gJ.
Let I1 denote the generalized complex structure induced by J, and let I1 denote the generalized complex structure induced by the symplectic form ω. We have
I1I2=[J00J][0ω1ω0]=[0g1g0]=I2I1

 Definition. A generalized Kähler manifold is a manifold with two commuting generalized complex structure I1,I2 such that the bilinear pairing (I1I2u,v) is positive definite.

Theorem (Gualtieri). A generalized Kähler structure on M induces a Riemannian metric g, two integrable almost complex structures J± Hermitian with respect to g, and two affine connections ± with skew-torsion ±H which preserve the metric and complex structure J±. Conversely, these data determine a generalized Kähler structure which is unique up to a B-field transformation.

Thus the notion of generalized Kähler manifold recovers the bihermitian geometry investigated by physicists in the context of susy non-linear σ-models.



Generalized Calabi-Yau Manifolds

Definition. A generalized Calabi-Yau manifold is a manifold M together with a complex-valued differential form ϕ, which is either purely even or purely odd, which is a pure spinor for the action of Cl(TM) and satisfies the non-degeneracy condition (ϕ,ˉϕ)0.

Note that (by definition) ϕ is pure if its annihilator is a maximal isotropic subspace. Let LTM be its annihilator. Then it is not hard to see that L defines a generalized complex structure on M, so indeed a generalized Calabi-Yau manifold is in particular a generalized complex manifold.

Example. If M is a complex manifold with a nowhere vanishing holomorphic (n,0) form, then it is generalized Calabi-Yau.

Example. If M is symplectic with symplectic form ω, then ϕ=exp(iω) gives M the structure of a generalized Calabi-Yau manifold.

If (M,ϕ) is generalized Calabi-Yau, then so is (M,exp(B)ϕ) for any closed real 2-form B. In the symplectic case, we obtain
ϕ=exp(B+iω)
This explains the appearance of the B-field (or "complexified Kähler form") in discussions of mirror symmetry.