Total space - M
 | (2.1.1) |
Constants:
Constants(Lambda);
 | (2.1.2) |
Vector fields:
Vectors(E[i],X,Y,Z);
![{X, Y, Z, E[i]}](Maple/atlas/examples/images/ricciflat_26.gif) | (2.1.3) |
Differential p-forms:
Forms(e[j]=1);
![{e[j]}](Maple/atlas/examples/images/ricciflat_27.gif) | (2.1.4) |
p-Sphere dimension (change it here):
p:=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])]](Maple/atlas/examples/images/ricciflat_29.gif) | (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])]](Maple/atlas/examples/images/ricciflat_30.gif) | (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[...](Maple/atlas/examples/images/ricciflat_31.gif)
![g = `+`(`/`(`*`(`&.`(e[1], e[1])), `*`(`+`(1, `-`(`/`(1, `*`(rho)))))), `*`(4, `*`(`+`(1, `-`(`/`(1, `*`(rho)))), `*`(`&.`(e[2], e[2])))), `/`(`*`(4, `*`(`^`(rho, 2), `*`(`+`(`&.`(e[3], e[3]), `&.`(e[...](Maple/atlas/examples/images/ricciflat_32.gif) | (2.1.8) |
Connection 1-forms:
Connection(omega);
![omega[i, j]](Maple/atlas/examples/images/ricciflat_33.gif) | (2.1.9) |
Curvature 2-forms:
Curvature(Omega);
![Omega[i, j]](Maple/atlas/examples/images/ricciflat_34.gif) | (2.1.10) |
Curvature tensor field:
Riemann(R);
Verify that total space is Ricci flat:
 | (2.1.12) |
Base space - B
Declare base domain:
Domain(B);
 | (2.2.1) |
Declare forms:
Forms(u[i]=1);
![{e[j], u[i]}](Maple/atlas/examples/images/ricciflat_43.gif) | (2.2.2) |
Declare vectors:
Vectors(U[k]);
![{U[k]}](Maple/atlas/examples/images/ricciflat_44.gif) | (2.2.3) |
Declare coframe
Coframe(u[1]=d(zeta),u[2]=d(xi));
![[u[1] = d(zeta), u[2] = d(xi)]](Maple/atlas/examples/images/ricciflat_45.gif) | (2.2.4) |
Frame declaration:
Frame(U[k]);
![[U[1] = Diff(``, zeta), U[2] = Diff(``, xi)]](Maple/atlas/examples/images/ricciflat_46.gif) | (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])))))](Maple/atlas/examples/images/ricciflat_47.gif) | (2.2.6) |
Submersion definition
Let us define submersion
:
such that
:
| > | Mapping(pi,M,B, zeta=rho, xi=theta); |
 |
 | (2.3.1) |
Projectors of the submersion
Now we can calculate vertical projector V and horizontal projector H:
![table( [( horizontal ) = `+`(`&.`(e[1], E[1]), `&.`(e[2], E[2])), ( vertical ) = `+`(`&.`(e[3], E[3]), `&.`(e[4], E[4])) ] )](Maple/atlas/examples/images/ricciflat_56.gif)
![table( [( horizontal ) = `+`(`&.`(e[1], E[1]), `&.`(e[2], E[2])), ( vertical ) = `+`(`&.`(e[3], E[3]), `&.`(e[4], E[4])) ] )](Maple/atlas/examples/images/ricciflat_57.gif) | (2.4.1) |
![`+`(`&.`(e[3], E[3]), `&.`(e[4], E[4]))](Maple/atlas/examples/images/ricciflat_58.gif) | (2.4.2) |
![`+`(`&.`(e[1], E[1]), `&.`(e[2], E[2]))](Maple/atlas/examples/images/ricciflat_59.gif) | (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](e[3]), `*`(E[3])), `*`(iota[X](e[4]), `*`(E[4])))](Maple/atlas/examples/images/ricciflat_60.gif) |
 = `+`(`*`(iota[X](e[1]), `*`(E[1])), `*`(iota[X](e[2]), `*`(E[2])))](Maple/atlas/examples/images/ricciflat_61.gif) | (2.4.4) |
Invariants T and A of the submersion:
Let us calculate invariants of the submersion:
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 `&/`:
![`+`(`/`(`*`(rho, `*`(`&.`(e[1], e[1]))), `*`(`+`(rho, `-`(1)))), `/`(`*`(4, `*`(`+`(rho, `-`(1)), `*`(`&.`(e[2], e[2])))), `*`(rho)))](Maple/atlas/examples/images/ricciflat_67.gif) | (2.5.2) |
We obtain the horizontal part of g metric.
For tensor field
we have:
![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]...](Maple/atlas/examples/images/ricciflat_69.gif)
![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]...](Maple/atlas/examples/images/ricciflat_70.gif) | (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]...](Maple/atlas/examples/images/ricciflat_71.gif)
![`+`(`/`(`*`(`&.`(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]...](Maple/atlas/examples/images/ricciflat_72.gif) | (2.5.4) |
For vectors X and Y:
| > | 'T[X](Y)'=iota[X,Y](TT); |
 = `+`(`/`(`*`(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...](Maple/atlas/examples/images/ricciflat_73.gif)
 = `+`(`/`(`*`(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...](Maple/atlas/examples/images/ricciflat_74.gif) | (2.5.5) |
Jump to total manifold:
 | (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...](Maple/atlas/examples/images/ricciflat_76.gif)
![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...](Maple/atlas/examples/images/ricciflat_77.gif) | (2.5.7) |
For mean curvature vector field
where
is vertical projection of the metric tensor
we obtain:
| > | N:=eval(Inv[meanCurvature]); |
![`+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2)))))](Maple/atlas/examples/images/ricciflat_81.gif) | (2.5.8) |
But for warped product we have
, where
is warped function. Let us verify that:
![N = `+`(`-`(`/`(`*`(2, `*`(`+`(rho, `-`(1)), `*`(E[1]))), `*`(`^`(rho, 2)))))](Maple/atlas/examples/images/ricciflat_84.gif) | (2.5.9) |