Schwarzschild black hole with cosmological constant 

Problem: 

Schwarzschild black hole with cosmological constant is 4-dimentional Lorentz manifold with constant Ricci curvature, timelike Killing vector field and group SO(3) as a subgroup of the manifold isometry group (with spacelike orbits).
For  Schwarzschild metric calculate the following: 

  • connetion 1-forms
 

  • curvature 2-forms
 

  • Riemannin tensor field
 

  • Ricci tensor field
 

Verify that are Killing vector vields. 

Schwarzschild metric 

> restart:
with(atlas):
 

Constants:
Constants(r[g],Lambda);
 

{`+`(`-`(I)), I, Pi, _Z, dim, Catalan, Lambda, r[g]}(2.1)
 

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

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

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

{e[j]}(2.3)
 

Coframe 1-forms:
Coframe(e[1]=d(t),e[2]=d(rho),e[3]=d(theta),e[4]=d(phi));
 

[e[1] = d(t), e[2] = d(rho), e[3] = d(theta), e[4] = d(phi)](2.4)
 

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

[E[1] = Diff(``, t), E[2] = Diff(``, rho), E[3] = Diff(``, theta), E[4] = Diff(``, phi)](2.5)
 

Metric tensor field:
Metric( g=(1-r[g]/rho+Lambda/3*rho^2)*d(t)&.d(t)-1/(1-r[g]/rho+Lambda/3*rho^2)*d(rho)&.d(rho)-rho^2*(d(theta)&.d(theta)+sin(theta)^2*d(phi)&.d(phi)) );
 

g = `+`(`*`(`+`(1, `-`(`/`(`*`(r[g]), `*`(rho))), `*`(`/`(1, 3), `*`(Lambda, `*`(`^`(rho, 2))))), `*`(`&.`(e[1], e[1]))), `-`(`/`(`*`(`&.`(e[2], e[2])), `*`(`+`(1, `-`(`/`(`*`(r[g]), `*`(rho))), `*`(`...
g = `+`(`*`(`+`(1, `-`(`/`(`*`(r[g]), `*`(rho))), `*`(`/`(1, 3), `*`(Lambda, `*`(`^`(rho, 2))))), `*`(`&.`(e[1], e[1]))), `-`(`/`(`*`(`&.`(e[2], e[2])), `*`(`+`(1, `-`(`/`(`*`(r[g]), `*`(rho))), `*`(`...
(2.6)
 

Connection 1-forms:
Connection(omega); 

omega[i, j](2.7)
 

> eval(omega);
 

table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
table( [( 4, 3 ) = `/`(`*`(cos(theta), `*`(e[4])), `*`(sin(theta))), ( 2, 2 ) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*`(Lambda, `*`(`^`(rho, 3)))), `*`(3, `*`(r[g]))), `*`(e[2]))), `*`(rho, `*`(...
(2.8)
 

Curvature 2-forms:
Curvature(Omega); 

Omega[i, j](2.9)
 

> eval(Omega);
 

table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
table( [( 4, 3 ) = `+`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&^`(e[3], e[4])))), `*`(rho))), ( 2, 2 ) = 0, ( 3, 2 ) = `+`(`/`(`*`(`/`(1, 2), `*`(`+`(`...
(2.10)
 

Curvature tensor field:
Riemann(R); 

R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
R = `+`(`-`(`/`(`*`(`/`(1, 3), `*`(`+`(`*`(Lambda, `*`(`^`(rho, 3))), `-`(`*`(3, `*`(r[g])))), `*`(`&.`(`&^`(e[1], e[2]), `&^`(e[1], e[2]))))), `*`(`^`(rho, 3)))), `/`(`*`(`/`(1, 2), `*`(`+`(`*`(2, `*...
(2.11)
 

Ricci tensor field calculation: 

> Ricci(ric);
 

ric = `+`(`/`(`*`(`/`(1, 3), `*`(Lambda, `*`(`+`(`*`(3, `*`(rho)), `-`(`*`(3, `*`(r[g]))), `*`(Lambda, `*`(`^`(rho, 3)))), `*`(`&.`(e[1], e[1]))))), `*`(rho)), `-`(`/`(`*`(3, `*`(rho, `*`(Lambda, `*`(...
ric = `+`(`/`(`*`(`/`(1, 3), `*`(Lambda, `*`(`+`(`*`(3, `*`(rho)), `-`(`*`(3, `*`(r[g]))), `*`(Lambda, `*`(`^`(rho, 3)))), `*`(`&.`(e[1], e[1]))))), `*`(rho)), `-`(`/`(`*`(3, `*`(rho, `*`(Lambda, `*`(...
(2.12)
 

Verify that metric g is Einstein one: 

> 'ric'-Lambda*g=simplify(ric-Lambda*ToBasis(g));
 

`+`(ric, `-`(`*`(Lambda, `*`(g)))) = 0(2.13)
 

Verify that and are Killing vector fields: 

> 'L[E[1]](g)'=L[E[1]](g);
 

L[E[1]](g) = 0(2.14)
 

> 'L[E[4]](g)'=L[E[4]](g);
 

L[E[4]](g) = 0(2.15)
 

>