CoCalc Shared Filesfix-tex / fixed.texOpen in CoCalc with one click!
Authors: J , Andy Huchala
Views : 13
1
2
\documentclass{article}
3
\usepackage{fullpage}
4
\usepackage{amsmath}
5
\usepackage[utf8x]{inputenc}
6
\usepackage{amssymb}
7
\usepackage{graphicx}
8
\usepackage{etoolbox}
9
\usepackage{url}
10
\usepackage{hyperref}
11
\usepackage[T1]{fontenc}
12
\makeatletter
13
\preto{\@verbatim}{\topsep=0pt \partopsep=0pt }
14
\makeatother
15
\usepackage{listings}
16
\lstdefinelanguage{Sage}[]{Python}
17
{morekeywords={True,False,sage,singular},
18
sensitive=true}
19
\lstset{
20
showtabs=False,
21
showspaces=False,
22
showstringspaces=False,
23
commentstyle={\ttfamily\color{dbrowncolor}},
24
keywordstyle={\ttfamily\color{dbluecolor}\bfseries},
25
stringstyle ={\ttfamily\color{dgraycolor}\bfseries},
26
backgroundcolor=\color{lightyellow},
27
language = Sage,
28
basicstyle={\ttfamily},
29
aboveskip=1em,
30
belowskip=0.1em,
31
breaklines=true,
32
prebreak = \raisebox{0ex}[0ex][0ex]{\ensuremath{\backslash}},
33
%frame=single
34
}
35
\usepackage{color}
36
\definecolor{lightyellow}{rgb}{1,1,.92}
37
\definecolor{dblackcolor}{rgb}{0.0,0.0,0.0}
38
\definecolor{dbluecolor}{rgb}{.01,.02,0.7}
39
\definecolor{dredcolor}{rgb}{1,0,0}
40
\definecolor{dbrowncolor}{rgb}{0.625,0.3125,0}
41
\definecolor{dgraycolor}{rgb}{0.30,0.3,0.30}
42
\definecolor{graycolor}{rgb}{0.35,0.35,0.35}
43
\title{SM\_de\_Sitter.sagews}
44
\author{Andy Huchala}
45
\begin{document}
46
\maketitle
47
48
{
49
\section{de Sitter spacetime}
50
{This worksheet demonstrates a few capabilities of \url{SageManifolds} (version 0.8) in computations regarding de Sitter spacetime.}\\\\
51
{It is released under the GNU General Public License version 3.}\\\\
52
{(c) Eric Gourgoulhon, Michal Bejger (2015)}\\\\
53
{The corresponding worksheet file can be downloaded from }{\url{here}}
54
\subsection{Spacetime manifold}
55
{We declare the de Sitter spacetime as a 4-dimensional differentiable manifold:}
56
57
\begin{lstlisting}
58
M = Manifold(4, 'M', r'\mathcal{M}')
59
print M ; M
60
\end{lstlisting}\begin{verbatim}4-dimensional manifold 'M'\end{verbatim}
61
62
{$\mathcal{M}$}\\\\
63
\iffalse
64
\begin{lstlisting}
65
%html
66
<p>We consider hyperspherical coordinates $(\tau,\chi,\theta,\phi)$ on $\mathcal{M}$. Allowing for the standard coordinate singularities at $\chi=0$, $\chi=\pi$, $\theta=0$ or $\theta=\pi$, these coordinates cover the entire spacetime manifold (which is topologically $\mathbb{R}\times\mathbb{S}^3$). If we restrict ourselves to <em>regular</em> coordinates (i.e. to consider only mathematically well defined charts), the hyperspherical coordinates cover only an open part of $\mathcal{M}$, which we call $\mathcal{M}_0$, on which $\chi$ spans the open interval $(0,\pi)$, $\theta$ the open interval $(0,\pi)$ and $\phi$ the open interval $(0,2\pi)$. Therefore, we declare:</p>
67
68
\end{lstlisting}\fi{We consider hyperspherical coordinates $(\tau,\chi,\theta,\phi)$ on $\mathcal{M}$. Allowing for the standard coordinate singularities at $\chi=0$, $\chi=\pi$, $\theta=0$ or $\theta=\pi$, these coordinates cover the entire spacetime manifold (which is topologically $\mathbb{R}\times\mathbb{S}^3$). If we restrict ourselves to {regular} coordinates (i.e. to consider only mathematically well defined charts), the hyperspherical coordinates cover only an open part of $\mathcal{M}$, which we call $\mathcal{M}_0$, on which $\chi$ spans the open interval $(0,\pi)$, $\theta$ the open interval $(0,\pi)$ and $\phi$ the open interval $(0,2\pi)$. Therefore, we declare:}
69
\begin{lstlisting}
70
M0 = M.open_subset('M_0', r'\mathcal{M}_0')
71
X_hyp.<ta,ch,th,ph> = M0.chart(r'ta:\tau ch:(0,pi):\chi th:(0,pi):\theta ph:(0,2*pi):\phi')
72
print X_hyp ; X_hyp
73
\end{lstlisting}\begin{verbatim}chart (M_0, (ta, ch, th, ph))\end{verbatim}
74
75
{$\left(\mathcal{M}_0,({\tau}, {\chi}, {\theta}, {\phi})\right)$}
76
\begin{lstlisting}
77
\end{lstlisting}\subsection{$\mathbb{R}^5$ as an ambient space}
78
{The de Sitter metric can be defined as that induced by the embedding of $\mathcal{M}$ into a 5-dimensional Minkowski space, i.e. $\mathbb{R}^5$ equipped with a flat Lorentzian metric. We therefore introduce $\mathbb{R}^5$ as a 5-dimensional manifold covered by canonical coordinates:}
79
\begin{lstlisting}
80
R5 = Manifold(5, 'R5', r'\mathbb{R}^5')
81
X5.<T,W,X,Y,Z> = R5.chart()
82
print X5 ; X5
83
\end{lstlisting}\begin{verbatim}chart (R5, (T, W, X, Y, Z))\end{verbatim}
84
85
{$\left(\mathbb{R}^5,(T, W, X, Y, Z)\right)$}
86
\begin{lstlisting}
87
\end{lstlisting}{The embedding of $\mathcal{M}$ into $\mathbb{R}^5$ is defined as a differential mapping $\Phi$ from $\mathcal{M}$ to $\mathbb{R}^5$, by providing its expression in terms of $\mathcal{M}$'s default chart (which is X\_hyp = $(\mathcal{M}_0,(\tau,\chi,\theta,\phi))$ ) and $\mathbb{R}^5$'s default chart (which is X5 = $(\mathbb{R}^5,(T,W,X,Y,Z))$ ):}
88
\begin{lstlisting}
89
var('b')
90
Phi = M.diff_mapping(R5, [sinh(b*ta)/b,
91
cosh(b*ta)/b * cos(ch),
92
cosh(b*ta)/b * sin(ch)*sin(th)*cos(ph),
93
cosh(b*ta)/b * sin(ch)*sin(th)*sin(ph),
94
cosh(b*ta)/b * sin(ch)*cos(th)],
95
name='Phi', latex_name=r'\Phi')
96
print Phi ; Phi.display()
97
\end{lstlisting}
98
99
{$b$}\begin{verbatim}differentiable mapping 'Phi' from the 4-dimensional manifold 'M' to the 5-dimensional
100
manifold 'R5'\end{verbatim}
101
102
{$\begin{array}{llcl} \Phi:& \mathcal{M} & \longrightarrow & \mathbb{R}^5 \\ \mbox{on}\ \mathcal{M}_0 : & \left({\tau}, {\chi}, {\theta}, {\phi}\right) & \longmapsto & \left(T, W, X, Y, Z\right) = \left(\frac{\sinh\left(b {\tau}\right)}{b}, \frac{\cos\left({\chi}\right) \cosh\left(b {\tau}\right)}{b}, \frac{\cos\left({\phi}\right) \cosh\left(b {\tau}\right) \sin\left({\chi}\right) \sin\left({\theta}\right)}{b}, \frac{\cosh\left(b {\tau}\right) \sin\left({\chi}\right) \sin\left({\phi}\right) \sin\left({\theta}\right)}{b}, \frac{\cos\left({\theta}\right) \cosh\left(b {\tau}\right) \sin\left({\chi}\right)}{b}\right) \end{array}$}
103
\begin{lstlisting}
104
\end{lstlisting}{The constant $b$ is a scale parameter. Considering de Sitter metric as a solution of vacuum Einstein equation with positive cosmological constant $\Lambda$, one has $b = \sqrt{\Lambda/3}$.}
105
{Let us evaluate the image of a point via the mapping $\Phi$:}
106
\begin{lstlisting}
107
p = M.point((ta, ch, th, ph), name='p') ; print p
108
\end{lstlisting}\begin{verbatim}point 'p' on 4-dimensional manifold 'M'\end{verbatim}
109
\begin{lstlisting}
110
p.coord()
111
\end{lstlisting}
112
113
{(${\tau}$, ${\chi}$, ${\theta}$, ${\phi}$)}
114
\begin{lstlisting}
115
q = Phi(p) ; print q
116
\end{lstlisting}\begin{verbatim}point 'Phi(p)' on 5-dimensional manifold 'R5'\end{verbatim}\begin{verbatim}\end{verbatim}
117
\begin{lstlisting}
118
q.coord()
119
\end{lstlisting}
120
121
{($\frac{\sinh\left(b {\tau}\right)}{b}$, $\frac{\cos\left({\chi}\right) \cosh\left(b {\tau}\right)}{b}$, $\frac{\cos\left({\phi}\right) \cosh\left(b {\tau}\right) \sin\left({\chi}\right) \sin\left({\theta}\right)}{b}$, $\frac{\cosh\left(b {\tau}\right) \sin\left({\chi}\right) \sin\left({\phi}\right) \sin\left({\theta}\right)}{b}$, $\frac{\cos\left({\theta}\right) \cosh\left(b {\tau}\right) \sin\left({\chi}\right)}{b}$)}
122
\begin{lstlisting}
123
\end{lstlisting}{The image of $\mathcal{M}$ by $\Phi$ is a hyperboloid of one sheet, of equation $-T^2+W^2+X^2+Y^2+Z^2=b^{-2}$. Indeed:}
124
\begin{lstlisting}
125
(Tq,Wq,Xq,Yq,Zq) = q.coord()
126
s = -Tq^2 + Wq^2 + Xq^2 + Yq^2 + Zq^2
127
s.simplify_full()
128
\end{lstlisting}
129
130
{$\frac{1}{b^{2}}$}
131
\begin{lstlisting}
132
\end{lstlisting}{We may use the embedding $\Phi$ to draw the coordinate grid $(\tau,\chi)$ in terms of the coordinates $(W,X,T)$ for $\theta=\pi/2$ and $\phi=0$ (red) and $\theta=\pi/2$ and $\phi=\pi$ (green) (the brown lines are the lines $\tau={\rm const}$):}
133
\begin{lstlisting}
134
graph1 = X_hyp.plot(X5, mapping=Phi, ambient_coords=(W,X,T), fixed_coords={th:pi/2, ph:0}, nb_values=9, color={ta:'red', ch:'brown'}, thickness=2, max_value=2, parameters={b:1}, label_axes=False)
135
graph2 = X_hyp.plot(X5, mapping=Phi, ambient_coords=(W,X,T), fixed_coords={th:pi/2, ph:pi}, nb_values=9, color={ta:'green', ch:'brown'}, thickness=2, max_value=2, parameters={b:1}, label_axes=False)
136
show(set_axes_labels(graph1+graph2,'W','X','T'), aspect_ratio=1)
137
\end{lstlisting}\url{https://cloud.sagemath.com/blobs/f97f4212-d46d-4ad8-b8f7-4979d61774e6.sage3d?uuid=f97f4212-d46d-4ad8-b8f7-4979d61774e6}
138
139
{}
140
\begin{lstlisting}
141
\end{lstlisting}\subsection{Spacetime metric}
142
{First, we introduce on $\mathbb{R}^5$ the Minkowski metric $h$:}
143
\begin{lstlisting}
144
h = R5.lorentz_metric('h')
145
h[0,0], h[1,1], h[2,2], h[3,3], h[4,4] = -1, 1, 1, 1, 1
146
h.display()
147
\end{lstlisting}
148
149
{$h = -\mathrm{d} T\otimes \mathrm{d} T+\mathrm{d} W\otimes \mathrm{d} W+\mathrm{d} X\otimes \mathrm{d} X+\mathrm{d} Y\otimes \mathrm{d} Y+\mathrm{d} Z\otimes \mathrm{d} Z$}
150
\begin{lstlisting}
151
\end{lstlisting}{As mentionned above, the de Sitter metric $g$ on $\mathcal{M}$ is that induced by $h$, i.e.$g$ is the pullback of $h$ by the mapping $\Phi$:}
152
\begin{lstlisting}
153
g = M.metric('g')
154
g.set( Phi.pullback(h) )
155
\end{lstlisting}
156
\begin{lstlisting}
157
\end{lstlisting}{The expression of $g$ in terms of $\mathcal{M}$'s default frame is found to be}
158
\begin{lstlisting}
159
g.display()
160
\end{lstlisting}
161
162
{$g = -\mathrm{d} {\tau}\otimes \mathrm{d} {\tau} + \frac{\cosh\left(b {\tau}\right)^{2}}{b^{2}} \mathrm{d} {\chi}\otimes \mathrm{d} {\chi} + \frac{\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2}}{b^{2}} \mathrm{d} {\theta}\otimes \mathrm{d} {\theta} + \frac{\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2}}{b^{2}} \mathrm{d} {\phi}\otimes \mathrm{d} {\phi}$}
163
\begin{lstlisting}
164
g[:]
165
\end{lstlisting}
166
167
{$\left(\begin{array}{rrrr}
168
-1 & 0 & 0 & 0 \\
169
0 & \frac{\cosh\left(b {\tau}\right)^{2}}{b^{2}} & 0 & 0 \\
170
0 & 0 & \frac{\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2}}{b^{2}} & 0 \\
171
0 & 0 & 0 & \frac{\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2}}{b^{2}}
172
\end{array}\right)$}
173
\begin{lstlisting}
174
\end{lstlisting}\subsection{Curvature}
175
{The Riemann tensor of $g$ is}
176
\begin{lstlisting}
177
Riem = g.riemann()
178
print Riem
179
Riem.display()
180
\end{lstlisting}\begin{verbatim}tensor field 'Riem(g)' of type (1,3) on the 4-dimensional manifold 'M'\end{verbatim}
181
182
{$\mathrm{Riem}\left(g\right) = \cosh\left(b {\tau}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\chi} -\cosh\left(b {\tau}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\tau} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\theta} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\tau} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\phi} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\tau} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\tau} + b^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\chi} -b^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\tau} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\theta} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\chi} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\phi} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\chi} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\chi} + b^{2} \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\theta} -b^{2} \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\tau} + \left( -\frac{\sin\left({\chi}\right)^{2} \sinh\left(b {\tau}\right)^{2} - \cos\left({\chi}\right)^{2} + 1}{\sin\left({\chi}\right)^{2}} \right) \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\theta} + \cosh\left(b {\tau}\right)^{2} \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\chi} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\phi} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \frac{\partial}{\partial {\theta} }\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\theta} + b^{2} \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\phi} -b^{2} \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\tau}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\tau} + \left( -\frac{\sin\left({\chi}\right)^{2} \sinh\left(b {\tau}\right)^{2} - \cos\left({\chi}\right)^{2} + 1}{\sin\left({\chi}\right)^{2}} \right) \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\phi} + \cosh\left(b {\tau}\right)^{2} \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\chi}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\chi} -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\phi} + \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \frac{\partial}{\partial {\phi} }\otimes \mathrm{d} {\theta}\otimes \mathrm{d} {\phi}\otimes \mathrm{d} {\theta}$}
183
\begin{lstlisting}
184
Riem.display_comp(only_nonredundant=True)
185
\end{lstlisting}
186
187
{$\begin{array}{lcl} \mathrm{Riem}\left(g\right)_{ \phantom{\, {\tau} } \, {\chi} \, {\tau} \, {\chi} }^{ \, {\tau} \phantom{\, {\chi} } \phantom{\, {\tau} } \phantom{\, {\chi} } } & = & \cosh\left(b {\tau}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\tau} } \, {\theta} \, {\tau} \, {\theta} }^{ \, {\tau} \phantom{\, {\theta} } \phantom{\, {\tau} } \phantom{\, {\theta} } } & = & \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\tau} } \, {\phi} \, {\tau} \, {\phi} }^{ \, {\tau} \phantom{\, {\phi} } \phantom{\, {\tau} } \phantom{\, {\phi} } } & = & \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\chi} } \, {\tau} \, {\tau} \, {\chi} }^{ \, {\chi} \phantom{\, {\tau} } \phantom{\, {\tau} } \phantom{\, {\chi} } } & = & b^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\chi} } \, {\theta} \, {\chi} \, {\theta} }^{ \, {\chi} \phantom{\, {\theta} } \phantom{\, {\chi} } \phantom{\, {\theta} } } & = & \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\chi} } \, {\phi} \, {\chi} \, {\phi} }^{ \, {\chi} \phantom{\, {\phi} } \phantom{\, {\chi} } \phantom{\, {\phi} } } & = & \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\theta} } \, {\tau} \, {\tau} \, {\theta} }^{ \, {\theta} \phantom{\, {\tau} } \phantom{\, {\tau} } \phantom{\, {\theta} } } & = & b^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\theta} } \, {\chi} \, {\chi} \, {\theta} }^{ \, {\theta} \phantom{\, {\chi} } \phantom{\, {\chi} } \phantom{\, {\theta} } } & = & -\frac{\sin\left({\chi}\right)^{2} \sinh\left(b {\tau}\right)^{2} - \cos\left({\chi}\right)^{2} + 1}{\sin\left({\chi}\right)^{2}} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\theta} } \, {\phi} \, {\theta} \, {\phi} }^{ \, {\theta} \phantom{\, {\phi} } \phantom{\, {\theta} } \phantom{\, {\phi} } } & = & \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\phi} } \, {\tau} \, {\tau} \, {\phi} }^{ \, {\phi} \phantom{\, {\tau} } \phantom{\, {\tau} } \phantom{\, {\phi} } } & = & b^{2} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\phi} } \, {\chi} \, {\chi} \, {\phi} }^{ \, {\phi} \phantom{\, {\chi} } \phantom{\, {\chi} } \phantom{\, {\phi} } } & = & -\frac{\sin\left({\chi}\right)^{2} \sinh\left(b {\tau}\right)^{2} - \cos\left({\chi}\right)^{2} + 1}{\sin\left({\chi}\right)^{2}} \\ \mathrm{Riem}\left(g\right)_{ \phantom{\, {\phi} } \, {\theta} \, {\theta} \, {\phi} }^{ \, {\phi} \phantom{\, {\theta} } \phantom{\, {\theta} } \phantom{\, {\phi} } } & = & -\cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \end{array}$}
188
\begin{lstlisting}
189
\end{lstlisting}{The Ricci tensor:}
190
\begin{lstlisting}
191
Ric = g.ricci()
192
print Ric
193
Ric.display()
194
\end{lstlisting}\begin{verbatim}field of symmetric bilinear forms 'Ric(g)' on the 4-dimensional manifold 'M'\end{verbatim}
195
196
{$\mathrm{Ric}\left(g\right) = -3 \, b^{2} \mathrm{d} {\tau}\otimes \mathrm{d} {\tau} + 3 \, \cosh\left(b {\tau}\right)^{2} \mathrm{d} {\chi}\otimes \mathrm{d} {\chi} + 3 \, \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \mathrm{d} {\theta}\otimes \mathrm{d} {\theta} + 3 \, \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2} \mathrm{d} {\phi}\otimes \mathrm{d} {\phi}$}
197
\begin{lstlisting}
198
Ric[:]
199
\end{lstlisting}
200
201
{$\left(\begin{array}{rrrr}
202
-3 \, b^{2} & 0 & 0 & 0 \\
203
0 & 3 \, \cosh\left(b {\tau}\right)^{2} & 0 & 0 \\
204
0 & 0 & 3 \, \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} & 0 \\
205
0 & 0 & 0 & 3 \, \cosh\left(b {\tau}\right)^{2} \sin\left({\chi}\right)^{2} \sin\left({\theta}\right)^{2}
206
\end{array}\right)$}
207
\begin{lstlisting}
208
\end{lstlisting}{The Ricci scalar:}
209
\begin{lstlisting}
210
R = g.ricci_scalar()
211
print R
212
R.display()
213
\end{lstlisting}\begin{verbatim}scalar field 'r(g)' on the 4-dimensional manifold 'M'\end{verbatim}
214
215
{$\begin{array}{llcl} \mathrm{r}\left(g\right):& \mathcal{M} & \longrightarrow & \mathbb{R} \\ \mbox{on}\ \mathcal{M}_0 : & \left({\tau}, {\chi}, {\theta}, {\phi}\right) & \longmapsto & 12 \, b^{2} \end{array}$}
216
\begin{lstlisting}
217
\end{lstlisting}{We recover the fact that de Sitter spacetime has a constant curvature. It is indeed a {maximally symmetric space}. In particular, the Riemann tensor is expressible as}
218
{\[ R^i_{\ \, jlk} = \frac{R}{n(n-1)} \left( \delta^i_{\ \, k} g_{jl} - \delta^i_{\ \, l} g_{jk} \right) \]}
219
{where $n$ is the dimension of $\mathcal{M}$: $n=4$ in the present case. Let us check this formula here, under the form $R^i_{\ \, jlk} = -\frac{R}{6} g_{j[k} \delta^i_{\ \, l]}$:}
220
\begin{lstlisting}
221
delta = M.tangent_identity_field()
222
Riem == - (R/6)*(g*delta).antisymmetrize(2,3) # 2,3 = last positions of the type-(1,3) tensor g*delta
223
\end{lstlisting}
224
225
{$\mathrm{True}$}
226
\begin{lstlisting}
227
\end{lstlisting}{We may also check that de Sitter metric is a solution of the vacuum {Einstein equation} with (positive) cosmological constant:}
228
\begin{lstlisting}
229
Lambda = 3*b^2
230
Ric - 1/2*R*g + Lambda*g == 0
231
\end{lstlisting}
232
233
{$\mathrm{True}$}
234
235
\end{document}
236
%sagemathcloud={"latex_command":"pdflatex -synctex=1 -interact=nonstopmode 'fixed.tex'"}