Ricci - flat warped product Copyright © 2003-2010 DigiArea Group . All rights reserved. Description: - Einstein manifolds (manifolds with constant Ricci curvature) are riemannian manifolds with metric tensor field
and Ricci tensor field where .(see Arthur L. Besse. "Einstein Manifolds" Springer-Verlag). Thus for Ricci flat manifolds we have and . - Warped products is a simple class of riemannian submersions which is defined as follows. Let
be riemannian manifold (base space) with metric and be riemannian manifold (fiber space) with metric . Then riemannian manifold (total space) where is positive function on is warped product with warped function . - In this worksheet we deal with warped product with 2-dimensional base
and metric: where is complete metric on p-dimensional Einstein manifold with Ricci constant and . We take p-dimensional sphere as the fiber space .

restart: with(atlas): Total space - M Domain(M); 
Constants: Constants(Lambda); 
Vector fields: Vectors(E[i],X,Y,Z); ![{E[i], X, Y, Z}](prod/atlas/examples/images/ricciflat27.gif)
Differential p-forms: Forms(e[j]=1); ![{e[j]}](prod/atlas/examples/images/ricciflat28.gif)
p-Sphere dimension (change it here): p:=2; 
Coframe 1-forms: Coframe(e[1]=d(rho),e[2]=d(theta),seq(e[i]=d(x[i-2]),i=3..p+2)); ![[e[1] = d(rho), e[2] = d(theta), e[3] = d(x[1]), e[4] = d(x[2])]](prod/atlas/examples/images/ricciflat30.gif)
Frame vector fields: Frame(E[i]); ![[E[1] = Diff(``,rho), E[2] = Diff(``,theta), E[3] = Diff(``,x[1]), E[4] = Diff(``,x[2])]](prod/atlas/examples/images/ricciflat31.gif)
Metric tensor fie ld: Metric( g=1/(1-rho^(1-p))*d(rho)&.d(rho)+4*(1-rho^(1-p))/(p-1)^2*d(theta)&.d(theta)+rho^2*4*add(d(x[i])&.d(x[i]),i=1..p)/(1+add(x[i]*x[i],i=1..p))^2); ![g = 1/(1-1/rho)*`&.`(e[1],e[1])+4*(1-1/rho)*`&.`(e[2],e[2])+4*rho^2*(`&.`(e[3],e[3])+`&.`(e[4],e[4]))/(1+x[1]^2+x[2]^2)^2](prod/atlas/examples/images/ricciflat32.gif)
Connection 1-forms: Connection(omega); ![omega[i,j]](prod/atlas/examples/images/ricciflat33.gif)
Curvature 2-forms: Curvature(Omega); ![Omega[i,j]](prod/atlas/examples/images/ricciflat34.gif)
Curvature tensor field: Riemann(R); ![R = 4/rho^3*`&.`(`&^`(e[1],e[2]),`&^`(e[1],e[2]))-2/(1+x[1]^2+x[2]^2)^2/(rho-1)*`&.`(`&^`(e[1],e[3]),`&^`(e[1],e[3]))-2/(1+x[1]^2+x[2]^2)^2/(rho-1)*`&.`(`&^`(e[1],e[4]),`&^`(e[1],e[4]))-8/rho^2/(1+x[1]...](prod/atlas/examples/images/ricciflat37.gif)
Verify that total space is Ricci flat: Ricci(r); 
Base space - B Declare base domain: Domain(B); 
Declare forms: Forms(u[i]=1); ![{e[j], u[i]}](prod/atlas/examples/images/ricciflat40.gif)
Declare vectors: Vectors(U[k]); ![{U[k]}](prod/atlas/examples/images/ricciflat41.gif)
Declare coframe Coframe(u[1]=d(zeta),u[2]=d(xi)); ![[u[1] = d(zeta), u[2] = d(xi)]](prod/atlas/examples/images/ricciflat42.gif)
Frame declaration: Frame(U[k]); ![[U[1] = Diff(``,zeta), U[2] = Diff(``,xi)]](prod/atlas/examples/images/ricciflat43.gif)
Let us define metric on the base: Metric(G=1/(1-zeta^(1-p))*d(zeta)&.d(zeta)+4*(1-zeta^(1-p))/(p-1)^2*d(xi)&.d(xi)); ![G = 1/(1-1/zeta)*`&.`(u[1],u[1])+4*(1-1/zeta)*`&.`(u[2],u[2])](prod/atlas/examples/images/ricciflat44.gif)
Submersion definition Let us define submersion : such that : Mapping(pi,M,B, zeta=rho, xi=theta); 

Who(pi);
pi: mapping
![TABLE([equations = [zeta = rho, xi = theta], frame = {E[4] = 0, E[3] = 0, E[1] = U[1], E[2] = U[2]}, coframe = {u[1] = e[1], u[2] = e[2]}, natural = {Diff(``,rho) = Diff(``,zeta), Diff(``,theta) = Diff...](prod/atlas/examples/images/ricciflat53.gif)
Projectors of the submersion Now we can calculate vertical projector V and horizontal projector H : P:=Projectors(pi); ![P := TABLE([horizontal = `&.`(e[1],E[1])+`&.`(e[2],E[2]), vertical = `&.`(e[3],E[3])+`&.`(e[4],E[4])])](prod/atlas/examples/images/ricciflat54.gif)
V:=P[vertical]; ![V := `&.`(e[3],E[3])+`&.`(e[4],E[4])](prod/atlas/examples/images/ricciflat55.gif)
H:=P[horizontal]; ![H := `&.`(e[1],E[1])+`&.`(e[2],E[2])](prod/atlas/examples/images/ricciflat56.gif)
Thus vertical and horizontal projections of arbitrary vector X are: 'iota[X](V)'=iota[X](V); 'iota[X](H)'=iota[X](H);  = iota[X](e[3])*E[3]+iota[X](e[4])*E[4]](prod/atlas/examples/images/ricciflat57.gif)
 = iota[X](e[1])*E[1]+iota[X](e[2])*E[2]](prod/atlas/examples/images/ricciflat58.gif)
Invariants T and A of the submersion: Let us calculate invariants of the submersion: Inv:=Invariants(pi); ![Inv := TABLE([integrabilityObstruction = 0, A = 0, T = TABLE(zero,[(3, 1) = [1/rho*E[3]], (4, 4) = [-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*E[1]], (4, 1) = [1/rho*E[4]], (3, 3) = [-4/(1+x[1]^2+x[2]^2)^2*(rho-1)...](prod/atlas/examples/images/ricciflat64.gif)
So, submersion invariant is equal to zero. Thus obstruction against integrability of the horizontal distribution is equal to zero. It is obvious that the submersion is a riemannian one but we can verify it directly. To do this we "rise" G metric into total space using restriction operator `&/`: G &/ pi; ![4*(rho-1)/rho*`&.`(e[2],e[2])+rho/(rho-1)*`&.`(e[1],e[1])](prod/atlas/examples/images/ricciflat66.gif)
We obtain the horizontal part of g metric. For tensor field we have: T:=eval(Inv['T']); ![T := TABLE(zero,[(3, 1) = [1/rho*E[3]], (4, 4) = [-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*E[1]], (4, 1) = [1/rho*E[4]], (3, 3) = [-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*E[1]]])](prod/atlas/examples/images/ricciflat70.gif)
To construct the T - tensor: TT:=add(add(`&.`(e[i],e[j],T[i,j]),j=1..p+2),i=1..p+2); ![TT := 1/rho*`&.`(e[3],e[1],E[3])-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*`&.`(e[3],e[3],E[1])+1/rho*`&.`(e[4],e[1],E[4])-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*`&.`(e[4],e[4],E[1])](prod/atlas/examples/images/ricciflat71.gif)
For vectors X and Y: 'T[X](Y)'=iota[X,Y](TT);  = 1/rho*iota[Y](e[3])*iota[X](e[1])*E[3]-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*iota[Y](e[3])*iota[X](e[3])*E[1]+1/rho*iota[Y](e[4])*iota[X](e[1])*E[4]-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*iota[Y](e[4])*iota[X...](prod/atlas/examples/images/ricciflat73.gif)
Jump to total manifold: Domain(M); 
Now for coordinate representation of T we obtain: 'iota[E[i],E[j]](T)'=iota[E[i],E[j]](TT); ![iota[E[i],E[j]](T) = 1/rho*delta[3,j]*delta[1,i]*E[3]-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*delta[3,j]*delta[3,i]*E[1]+1/rho*delta[4,j]*delta[1,i]*E[4]-4/(1+x[1]^2+x[2]^2)^2*(rho-1)*delta[4,j]*delta[4,i]*E[1]](prod/atlas/examples/images/ricciflat75.gif)
For mean curvature vector field where is vertical projection of the metric tensor we obtain: N:=eval(Inv[meanCurvature]); ![N := -2/rho^2*(rho-1)*E[1]](prod/atlas/examples/images/ricciflat79.gif)
But for warped product we have , where is warped function. Let us verify that: 'N'=-p/rho*grad(rho); ![N = -2/rho^2*(rho-1)*E[1]](prod/atlas/examples/images/ricciflat82.gif)
|