Ricci - flat warped product 

Description: 

  • Einstein manifolds (manifolds with constant Ricci curvature) are riemannian manifolds with metric tensor field g and Ricci tensor field r = `*`(lambda, `*`(g)) where lambda = const.(see Arthur L. Besse. "Einstein Manifolds" Springer-Verlag). Thus for Ricci flat manifolds we have lambda = 0 and  r = 0.
 

  • Warped products is a simple class of riemannian submersions which are defined as follows.
    Let {B, g[B]} be riemannian manifold B (base space) with metric g[B] and {F, g[F]} be riemannian manifold F (fiber space) with metric g[F].
    The warped product is the riemannian manifold {`+`(g[B], `*`(f, `*`(g[F]))), Typesetting:-delayCrossProduct(B, F)}  where f = f(b) is positive function (warped function) on
 

  • In this worksheet we deal with warped product with 2-dimensional base B = `*`(`^`(R, 2)) and metric:
    g = `+`(`/`(`*`(`⊗`(d(r), d(r))), `*`(`+`(1, `-`(`^`(r, `+`(1, `-`(p))))))), `/`(`*`(`*`(4, `+`(1, `-`(`^`(r, `+`(1, `-`(p)))))), `*`(d(theta), `*`(`⊗`(d(theta))))), `*`(`^`(`+`(p, `-`(1...  where g[F] is complete metric on p-dimensional Einstein manifold F with Ricci constant lambda[F] = `+`(p, `-`(1)) and p = dim(F). We take p-dimensional sphere {g[can], `^`(S, p)} as the fiber space {F, g[F]}.
 

`*`(`^`(R, 2), `*`(cross, `*`(`^`(S, p)))) 

> restart:
with(atlas):
 

Total space - M 

> Domain(M);
 

M(2.1.1)
 

Constants:
Constants(Lambda);
 

{`+`(`-`(I)), I, Pi, _Z, Catalan, Lambda}(2.1.2)
 

Vector fields:
Vectors(E[i],X,Y,Z);
 

{X, Y, Z, E[i]}(2.1.3)
 

Differential p-forms:
Forms(e[j]=1);
 

{e[j]}(2.1.4)
 

p-Sphere dimension (change it here):
p:=2;
 

2(2.1.5)
 

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])](2.1.6)
 

Frame vector fields:
Frame(E[i]);
 

[E[1] = Diff(``, rho), E[2] = Diff(``, theta), E[3] = Diff(``, x[1]), E[4] = Diff(``, x[2])](2.1.7)
 

Metric tensor field:
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 = `+`(`/`(`*`(`&.`(e[1], e[1])), `*`(`+`(1, `-`(`/`(1, `*`(rho)))))), `*`(4, `*`(`+`(1, `-`(`/`(1, `*`(rho)))), `*`(`&.`(e[2], e[2])))), `/`(`*`(4, `*`(`^`(rho, 2), `*`(`+`(`&.`(e[3], e[3]), `&.`(e[...
g = `+`(`/`(`*`(`&.`(e[1], e[1])), `*`(`+`(1, `-`(`/`(1, `*`(rho)))))), `*`(4, `*`(`+`(1, `-`(`/`(1, `*`(rho)))), `*`(`&.`(e[2], e[2])))), `/`(`*`(4, `*`(`^`(rho, 2), `*`(`+`(`&.`(e[3], e[3]), `&.`(e[...
(2.1.8)
 

Connection 1-forms:
Connection(omega); 

omega[i, j](2.1.9)
 

Curvature 2-forms:
Curvature(Omega); 

Omega[i, j](2.1.10)
 

Curvature tensor field:
Riemann(R); 

R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
R = `+`(`/`(`*`(4, `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2])))), `*`(`^`(rho, 3))), `-`(`/`(`*`(2, `*`(`&.`(`&^`(e[1], e[3]), `&^`(e[1], e[3])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))...
(2.1.11)
 

Verify that total space is Ricci flat: 

> Ricci(r);
 

r = 0(2.1.12)
 

Base space - B 

Declare base domain:
Domain(B); 

B(2.2.1)
 

Declare forms:
Forms(u[i]=1); 

{e[j], u[i]}(2.2.2)
 

Declare vectors:
Vectors(U[k]); 

{U[k]}(2.2.3)
 

Declare coframe
Coframe(u[1]=d(zeta),u[2]=d(xi)); 

[u[1] = d(zeta), u[2] = d(xi)](2.2.4)
 

Frame declaration:
Frame(U[k]); 

[U[1] = Diff(``, zeta), U[2] = Diff(``, xi)](2.2.5)
 

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 = `+`(`/`(`*`(`&.`(u[1], u[1])), `*`(`+`(1, `-`(`/`(1, `*`(zeta)))))), `*`(4, `*`(`+`(1, `-`(`/`(1, `*`(zeta)))), `*`(`&.`(u[2], u[2])))))(2.2.6)
 

Submersion definition 

Let us define submersion  Pi: proc (M) options operator, arrow; B end proc such that  Zeta = rho, xi = theta: 

> Mapping(pi,M,B,
       zeta=rho,
       xi=theta);
 

 

pi
`.`(M, `--->`, B)(2.3.1)
 

> Who(pi);
 

 

pi: mapping
table( [( equations ) = [zeta = rho, xi = theta], ( manifolds ) = [M, B], ( natural ) = {Diff(``, rho) = Diff(``, zeta), Diff(``, theta) = Diff(``, xi), Diff(``, x[1]) = 0, Diff(``, x[2]) = 0}, ( fram...
table( [( equations ) = [zeta = rho, xi = theta], ( manifolds ) = [M, B], ( natural ) = {Diff(``, rho) = Diff(``, zeta), Diff(``, theta) = Diff(``, xi), Diff(``, x[1]) = 0, Diff(``, x[2]) = 0}, ( fram...
table( [( equations ) = [zeta = rho, xi = theta], ( manifolds ) = [M, B], ( natural ) = {Diff(``, rho) = Diff(``, zeta), Diff(``, theta) = Diff(``, xi), Diff(``, x[1]) = 0, Diff(``, x[2]) = 0}, ( fram...
(2.3.2)
 

Projectors of the submersion 

Now we can calculate vertical projector V and horizontal projector H: 

> P:=Projectors(pi);
 

table( [( horizontal ) = `+`(`&.`(e[1], E[1]), `&.`(e[2], E[2])), ( vertical ) = `+`(`&.`(e[3], E[3]), `&.`(e[4], E[4])) ] )
table( [( horizontal ) = `+`(`&.`(e[1], E[1]), `&.`(e[2], E[2])), ( vertical ) = `+`(`&.`(e[3], E[3]), `&.`(e[4], E[4])) ] )
(2.4.1)
 

> V:=P[vertical];
 

`+`(`&.`(e[3], E[3]), `&.`(e[4], E[4]))(2.4.2)
 

> H:=P[horizontal];
 

`+`(`&.`(e[1], E[1]), `&.`(e[2], E[2]))(2.4.3)
 

Thus vertical and horizontal projections of arbitrary vector X are: 

> 'iota[X](V)'=iota[X](V);
'iota[X](H)'=iota[X](H);
 

 

iota[X](V) = `+`(`*`(iota[X](e[3]), `*`(E[3])), `*`(iota[X](e[4]), `*`(E[4])))
iota[X](H) = `+`(`*`(iota[X](e[1]), `*`(E[1])), `*`(iota[X](e[2]), `*`(E[2])))(2.4.4)
 

Invariants T and A of the submersion: 

Let us calculate invariants of the submersion: 

> Inv:=Invariants(pi);
 

table( [( integrabilityObstruction ) = 0, ( A ) = 0, ( meanCurvature ) = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2))))), ( riemannianObstruction ) = 0, ( T ) = table( [( 4, ...
table( [( integrabilityObstruction ) = 0, ( A ) = 0, ( meanCurvature ) = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2))))), ( riemannianObstruction ) = 0, ( T ) = table( [( 4, ...
table( [( integrabilityObstruction ) = 0, ( A ) = 0, ( meanCurvature ) = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2))))), ( riemannianObstruction ) = 0, ( T ) = table( [( 4, ...
table( [( integrabilityObstruction ) = 0, ( A ) = 0, ( meanCurvature ) = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2))))), ( riemannianObstruction ) = 0, ( T ) = table( [( 4, ...
(2.5.1)
 

So, submersion invariant A[X](Y) = `+`(`*`(H, `*`(D[`*`(H, `*`(X))](`*`(V, `*`(Y))))), `*`(V, `*`(D[`*`(H, `*`(X))](`*`(H, `*`(Y)))))) 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;
 

`+`(`/`(`*`(rho, `*`(`&.`(e[1], e[1]))), `*`(`+`(rho, `-`(1)))), `/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(`&.`(e[2], e[2])))), `*`(rho)))(2.5.2)
 

We obtain the horizontal part of g metric. 

For tensor field T[X](Y) = `+`(`*`(H, `*`(D[`*`(V, `*`(X))](`*`(V, `*`(Y))))), `*`(V, `*`(D[`*`(V, `*`(X))](`*`(H, `*`(Y)))))) we have: 

> T:=eval(Inv['T']);
 

table( [( 4, 1 ) = [`/`(`*`(E[4]), `*`(rho))], ( 3, 3 ) = [`+`(`-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))), 2)))))], ( 3, 1 ) = [`/`(`*`(E[3]...
table( [( 4, 1 ) = [`/`(`*`(E[4]), `*`(rho))], ( 3, 3 ) = [`+`(`-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))), 2)))))], ( 3, 1 ) = [`/`(`*`(E[3]...
(2.5.3)
 

To construct the T - tensor:  

> TT:=add(add(`&.`(e[i],e[j],T[i,j]),j=1..p+2),i=1..p+2);
 

`+`(`/`(`*`(`&.`(e[3], e[1], E[3])), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(`&.`(e[3], e[3], E[1])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))), 2)))), `/`(`*`(`&.`(e[4], e[1]...
`+`(`/`(`*`(`&.`(e[3], e[1], E[3])), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(`&.`(e[3], e[3], E[1])))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2)), `*`(`^`(x[2], 2))), 2)))), `/`(`*`(`&.`(e[4], e[1]...
(2.5.4)
 

For vectors X and Y: 

> 'T[X](Y)'=iota[X,Y](TT);
 

T[X](Y) = `+`(`/`(`*`(iota[Y](e[3]), `*`(iota[X](e[1]), `*`(E[3]))), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(iota[Y](e[3]), `*`(iota[X](e[3]), `*`(E[1]))))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2...
T[X](Y) = `+`(`/`(`*`(iota[Y](e[3]), `*`(iota[X](e[1]), `*`(E[3]))), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(iota[Y](e[3]), `*`(iota[X](e[3]), `*`(E[1]))))), `*`(`^`(`+`(1, `*`(`^`(x[1], 2...
(2.5.5)
 

Jump to total manifold: 

> Domain(M);
 

M(2.5.6)
 

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) = `+`(`/`(`*`(delta[3, j], `*`(delta[1, i], `*`(E[3]))), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(delta[3, j], `*`(delta[3, i], `*`(E[1]))))), `*`(`^`(`+`(1, `*`(`^`(x[1...
iota[E[i], E[j]](T) = `+`(`/`(`*`(delta[3, j], `*`(delta[1, i], `*`(E[3]))), `*`(rho)), `-`(`/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(delta[3, j], `*`(delta[3, i], `*`(E[1]))))), `*`(`^`(`+`(1, `*`(`^`(x[1...
(2.5.7)
 

For mean curvature vector field N = trace[`*`(V, `*`(g))](T) where `*`(V, `*`(g)) is vertical projection of the metric tensor g we obtain: 

> N:=eval(Inv[meanCurvature]);
 

`+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2)))))(2.5.8)
 

But for warped product we have N = `+`(`-`(`/`(`*`(p, `*`(grad(sqrt(f)))), `*`(sqrt(f))))) , where f = `*`(`^`(rho, 2)) is warped function. Let us verify that: 

> 'N'=-p/rho*grad(rho);
 

N = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2)))))(2.5.9)
 

>