{ "cells": [ { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "# 5D Kerr-AdS spacetime: general solution with 2 angular momenta" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "*NB:* a version of SageMath at least equal to 7.5 is required to run this worksheet:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/plain": [ "'SageMath version 8.0.beta12, Release Date: 2017-06-22'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "version()" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "First we set up the notebook to display mathematical objects using LaTeX rendering:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "%display latex" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "We also define a viewer for 3D plots (use `'threejs'` or `'jmol'` for interactive 3D graphics):" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "viewer3D = 'threejs' # must be 'threejs', jmol', 'tachyon' or None (default)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Since some computations are quite long, we ask for running them in parallel on 8 cores:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "Parallelism().set(nproc=8)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "## Spacetime manifold\n", "\n", "We declare the Kerr-AdS spacetime as a 5-dimensional diffentiable manifold:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "5-dimensional differentiable manifold M\n" ] } ], "source": [ "M = Manifold(5, 'M', r'\\mathcal{M}')\n", "print(M)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Let us define **Boyer-Lindquist-type coordinates (rational polynomial version)** on $\\mathcal{M}$, via the method `chart()`, the argument of which is a string expressing the coordinates names, their ranges (the default is $(-\\infty,+\\infty)$) and their LaTeX symbols:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(\\mathcal{M},(t, r, {\\mu}, {\\phi}, {\\psi})\\right)$ " ], "text/plain": [ "Chart (M, (t, r, mu, ph, ps))" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "BL. = M.chart(r't r:(0,+oo) mu:(-1,1):\\mu ph:(0,2*pi):\\phi ps:(0,2*pi):\\psi')\n", "BL" ] }, { "cell_type": "markdown", "metadata": { "collapsed": false, "deletable": true, "editable": true }, "source": [ "Note that $\\mu$ is related to the standard Boyer-Lindquist coordinate $\\theta$ by\n", "$$ \\mu = \\cos\\theta$$" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "## Metric tensor\n", "\n", "The 4 parameters $m$, $a$, $b$ and $\\ell$ of the Kerr-AdS spacetime are declared as symbolic variables, $a$ and $b$ being the two angular momentum parameters and $\\ell$ being related to the cosmological constant by $\\Lambda = - 6 \\ell^2$:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\left(m, a, b\\right)$ " ], "text/plain": [ "(m, a, b)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "var('m a b', domain='real')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}{\\ell}$ " ], "text/plain": [ "l" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "var('l', domain='real', latex_name=r'\\ell')" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "# Particular cases\n", "# a = 0\n", "# m = 0\n", "# b = a" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Some auxiliary functions:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "sig = (1+r^2*l^2)/r^2\n", "Delta = (r^2+a^2)*(r^2+b^2)*sig - 2*m\n", "sinth2 = 1-mu^2\n", "Delta_th = 1 - a^2*l^2*mu^2 - b^2*l^2*sinth2\n", "rho2 = r^2 + a^2*mu^2 + b^2*sinth2\n", "Xi_a = 1 - a^2*l^2\n", "Xi_b = 1 - b^2*l^2" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "The metric is set by its components in the coordinate frame associated with the Boyer-Lindquist-type coordinates, which is the current manifold's default frame. These components are given by \n", "Eq. (5.22) of the article [S.W. Hawking, C.J. Hunter & M.M. Taylor-Robinson, Phys. Rev. D **59**, 064005 (1999)](https://doi.org/10.1103/PhysRevD.59.064005):" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}g = \\left( \\frac{{\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} {\\mu}^{4} - a^{2} b^{2} {\\ell}^{2} - {\\ell}^{2} r^{4} - {\\left({\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} - {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} - b^{2} + 2 \\, m}{{\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2} + r^{2}} \\right) \\mathrm{d} t\\otimes \\mathrm{d} t + \\left( -\\frac{a^{3} b^{2} {\\ell}^{2} - {\\left(a^{5} - a^{3} b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a {\\ell}^{2} {\\mu}^{2} - a {\\ell}^{2}\\right)} r^{4} + {\\left({\\left(a^{5} - 2 \\, a^{3} b^{2}\\right)} {\\ell}^{2} + 2 \\, a m\\right)} {\\mu}^{2} - {\\left(2 \\, a b^{2} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{3} - a b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a^{3} + a b^{2}\\right)} {\\ell}^{2}\\right)} r^{2} - 2 \\, a m}{a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\right) \\mathrm{d} t\\otimes \\mathrm{d} {\\phi} + \\left( -\\frac{b {\\ell}^{2} {\\mu}^{2} r^{4} + {\\left(a^{2} b^{3} - b^{5}\\right)} {\\ell}^{2} {\\mu}^{4} + {\\left(b^{5} {\\ell}^{2} - 2 \\, b m\\right)} {\\mu}^{2} + {\\left(2 \\, b^{3} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{2} b - b^{3}\\right)} {\\ell}^{2} {\\mu}^{4}\\right)} r^{2}}{b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(b^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\right) \\mathrm{d} t\\otimes \\mathrm{d} {\\psi} + \\left( \\frac{r^{4} + {\\left({\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2}\\right)} r^{2}}{{\\ell}^{2} r^{6} + {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{4} + a^{2} b^{2} + {\\left(a^{2} b^{2} {\\ell}^{2} + a^{2} + b^{2} - 2 \\, m\\right)} r^{2}} \\right) \\mathrm{d} r\\otimes \\mathrm{d} r + \\left( \\frac{{\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2} + r^{2}}{{\\left(a^{2} - b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - b^{2} {\\ell}^{2} - {\\left({\\left(a^{2} - 2 \\, b^{2}\\right)} {\\ell}^{2} + 1\\right)} {\\mu}^{2} + 1} \\right) \\mathrm{d} {\\mu}\\otimes \\mathrm{d} {\\mu} + \\left( -\\frac{a^{3} b^{2} {\\ell}^{2} - {\\left(a^{5} - a^{3} b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a {\\ell}^{2} {\\mu}^{2} - a {\\ell}^{2}\\right)} r^{4} + {\\left({\\left(a^{5} - 2 \\, a^{3} b^{2}\\right)} {\\ell}^{2} + 2 \\, a m\\right)} {\\mu}^{2} - {\\left(2 \\, a b^{2} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{3} - a b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a^{3} + a b^{2}\\right)} {\\ell}^{2}\\right)} r^{2} - 2 \\, a m}{a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\right) \\mathrm{d} {\\phi}\\otimes \\mathrm{d} t + \\left( -\\frac{a^{4} b^{2} {\\ell}^{2} + {\\left(a^{4} - a^{2} b^{2} - {\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{2} - 2 \\, a^{2} m\\right)} {\\mu}^{4} + {\\left(a^{2} {\\ell}^{2} - {\\left(a^{2} {\\ell}^{2} - 1\\right)} {\\mu}^{2} - 1\\right)} r^{4} - a^{2} b^{2} - 2 \\, a^{2} m - {\\left(a^{4} - 2 \\, a^{2} b^{2} - {\\left(a^{6} - 2 \\, a^{4} b^{2}\\right)} {\\ell}^{2} - 4 \\, a^{2} m\\right)} {\\mu}^{2} - {\\left({\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{4} - {\\left(a^{4} + a^{2} b^{2}\\right)} {\\ell}^{2} + 2 \\, {\\left(a^{2} b^{2} {\\ell}^{2} - b^{2}\\right)} {\\mu}^{2} + a^{2} + b^{2}\\right)} r^{2}}{a^{4} b^{2} {\\ell}^{4} - 2 \\, a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{4} {\\ell}^{4} - 2 \\, a^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\right) \\mathrm{d} {\\phi}\\otimes \\mathrm{d} {\\phi} + \\left( -\\frac{2 \\, {\\left(a b m {\\mu}^{4} - a b m {\\mu}^{2}\\right)}}{a^{2} b^{4} {\\ell}^{4} - {\\left(a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + {\\left({\\left(a^{4} b^{2} - a^{2} b^{4}\\right)} {\\ell}^{4} - {\\left(a^{4} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} b^{2} {\\ell}^{4} - {\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\right) \\mathrm{d} {\\phi}\\otimes \\mathrm{d} {\\psi} + \\left( -\\frac{b {\\ell}^{2} {\\mu}^{2} r^{4} + {\\left(a^{2} b^{3} - b^{5}\\right)} {\\ell}^{2} {\\mu}^{4} + {\\left(b^{5} {\\ell}^{2} - 2 \\, b m\\right)} {\\mu}^{2} + {\\left(2 \\, b^{3} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{2} b - b^{3}\\right)} {\\ell}^{2} {\\mu}^{4}\\right)} r^{2}}{b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(b^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\right) \\mathrm{d} {\\psi}\\otimes \\mathrm{d} t + \\left( -\\frac{2 \\, {\\left(a b m {\\mu}^{4} - a b m {\\mu}^{2}\\right)}}{a^{2} b^{4} {\\ell}^{4} - {\\left(a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + {\\left({\\left(a^{4} b^{2} - a^{2} b^{4}\\right)} {\\ell}^{4} - {\\left(a^{4} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} b^{2} {\\ell}^{4} - {\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\right) \\mathrm{d} {\\psi}\\otimes \\mathrm{d} {\\phi} + \\left( -\\frac{{\\left(b^{2} {\\ell}^{2} - 1\\right)} {\\mu}^{2} r^{4} - {\\left(a^{2} b^{2} - b^{4} - {\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{2} + 2 \\, b^{2} m\\right)} {\\mu}^{4} + {\\left(b^{6} {\\ell}^{2} - b^{4}\\right)} {\\mu}^{2} + {\\left({\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{4} + 2 \\, {\\left(b^{4} {\\ell}^{2} - b^{2}\\right)} {\\mu}^{2}\\right)} r^{2}}{b^{6} {\\ell}^{4} - 2 \\, b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(b^{4} {\\ell}^{4} - 2 \\, b^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\right) \\mathrm{d} {\\psi}\\otimes \\mathrm{d} {\\psi}$ " ], "text/plain": [ "g = ((a^4 - 2*a^2*b^2 + b^4)*l^2*mu^4 - a^2*b^2*l^2 - l^2*r^4 - ((a^4 - 2*a^2*b^2 + b^4)*l^2 + a^2 - b^2)*mu^2 - ((a^2 + b^2)*l^2 + 1)*r^2 - b^2 + 2*m)/((a^2 - b^2)*mu^2 + b^2 + r^2) dt*dt - (a^3*b^2*l^2 - (a^5 - a^3*b^2)*l^2*mu^4 - (a*l^2*mu^2 - a*l^2)*r^4 + ((a^5 - 2*a^3*b^2)*l^2 + 2*a*m)*mu^2 - (2*a*b^2*l^2*mu^2 + (a^3 - a*b^2)*l^2*mu^4 - (a^3 + a*b^2)*l^2)*r^2 - 2*a*m)/(a^2*b^2*l^2 + ((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^2 + (a^2*l^2 - 1)*r^2 - b^2) dt*dph - (b*l^2*mu^2*r^4 + (a^2*b^3 - b^5)*l^2*mu^4 + (b^5*l^2 - 2*b*m)*mu^2 + (2*b^3*l^2*mu^2 + (a^2*b - b^3)*l^2*mu^4)*r^2)/(b^4*l^2 + ((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^2 + (b^2*l^2 - 1)*r^2 - b^2) dt*dps + (r^4 + ((a^2 - b^2)*mu^2 + b^2)*r^2)/(l^2*r^6 + ((a^2 + b^2)*l^2 + 1)*r^4 + a^2*b^2 + (a^2*b^2*l^2 + a^2 + b^2 - 2*m)*r^2) dr*dr + ((a^2 - b^2)*mu^2 + b^2 + r^2)/((a^2 - b^2)*l^2*mu^4 - b^2*l^2 - ((a^2 - 2*b^2)*l^2 + 1)*mu^2 + 1) dmu*dmu - (a^3*b^2*l^2 - (a^5 - a^3*b^2)*l^2*mu^4 - (a*l^2*mu^2 - a*l^2)*r^4 + ((a^5 - 2*a^3*b^2)*l^2 + 2*a*m)*mu^2 - (2*a*b^2*l^2*mu^2 + (a^3 - a*b^2)*l^2*mu^4 - (a^3 + a*b^2)*l^2)*r^2 - 2*a*m)/(a^2*b^2*l^2 + ((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^2 + (a^2*l^2 - 1)*r^2 - b^2) dph*dt - (a^4*b^2*l^2 + (a^4 - a^2*b^2 - (a^6 - a^4*b^2)*l^2 - 2*a^2*m)*mu^4 + (a^2*l^2 - (a^2*l^2 - 1)*mu^2 - 1)*r^4 - a^2*b^2 - 2*a^2*m - (a^4 - 2*a^2*b^2 - (a^6 - 2*a^4*b^2)*l^2 - 4*a^2*m)*mu^2 - (((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^4 - (a^4 + a^2*b^2)*l^2 + 2*(a^2*b^2*l^2 - b^2)*mu^2 + a^2 + b^2)*r^2)/(a^4*b^2*l^4 - 2*a^2*b^2*l^2 + ((a^6 - a^4*b^2)*l^4 - 2*(a^4 - a^2*b^2)*l^2 + a^2 - b^2)*mu^2 + (a^4*l^4 - 2*a^2*l^2 + 1)*r^2 + b^2) dph*dph - 2*(a*b*m*mu^4 - a*b*m*mu^2)/(a^2*b^4*l^4 - (a^2*b^2 + b^4)*l^2 + ((a^4*b^2 - a^2*b^4)*l^4 - (a^4 - b^4)*l^2 + a^2 - b^2)*mu^2 + (a^2*b^2*l^4 - (a^2 + b^2)*l^2 + 1)*r^2 + b^2) dph*dps - (b*l^2*mu^2*r^4 + (a^2*b^3 - b^5)*l^2*mu^4 + (b^5*l^2 - 2*b*m)*mu^2 + (2*b^3*l^2*mu^2 + (a^2*b - b^3)*l^2*mu^4)*r^2)/(b^4*l^2 + ((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^2 + (b^2*l^2 - 1)*r^2 - b^2) dps*dt - 2*(a*b*m*mu^4 - a*b*m*mu^2)/(a^2*b^4*l^4 - (a^2*b^2 + b^4)*l^2 + ((a^4*b^2 - a^2*b^4)*l^4 - (a^4 - b^4)*l^2 + a^2 - b^2)*mu^2 + (a^2*b^2*l^4 - (a^2 + b^2)*l^2 + 1)*r^2 + b^2) dps*dph - ((b^2*l^2 - 1)*mu^2*r^4 - (a^2*b^2 - b^4 - (a^2*b^4 - b^6)*l^2 + 2*b^2*m)*mu^4 + (b^6*l^2 - b^4)*mu^2 + (((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^4 + 2*(b^4*l^2 - b^2)*mu^2)*r^2)/(b^6*l^4 - 2*b^4*l^2 + ((a^2*b^4 - b^6)*l^4 - 2*(a^2*b^2 - b^4)*l^2 + a^2 - b^2)*mu^2 + (b^4*l^4 - 2*b^2*l^2 + 1)*r^2 + b^2) dps*dps" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "g = M.lorentzian_metric('g')\n", "tmp = 1/rho2*( -Delta + Delta_th*(a^2*sinth2 + b^2*mu^2) + a^2*b^2*sig )\n", "g[0,0] = tmp.simplify_full()\n", "tmp = a*sinth2/(rho2*Xi_a)*( Delta - (r^2+a^2)*(Delta_th + b^2*sig) )\n", "g[0,3] = tmp.simplify_full()\n", "tmp = b*mu^2/(rho2*Xi_b)*( Delta - (r^2+b^2)*(Delta_th + a^2*sig) )\n", "g[0,4] = tmp.simplify_full()\n", "g[1,1] = (rho2/Delta).simplify_full()\n", "g[2,2] = (rho2/Delta_th/(1-mu^2)).simplify_full()\n", "tmp = sinth2/(rho2*Xi_a^2)*( - Delta*a^2*sinth2 + (r^2+a^2)^2*(Delta_th + sig*b^2*sinth2) ) \n", "g[3,3] = tmp.simplify_full()\n", "tmp = a*b*sinth2*mu^2/(rho2*Xi_a*Xi_b)*( - Delta + sig*(r^2+a^2)*(r^2+b^2) )\n", "g[3,4] = tmp.simplify_full()\n", "tmp = mu^2/(rho2*Xi_b^2)*( - Delta*b^2*mu^2 + (r^2+b^2)^2*(Delta_th + sig*a^2*mu^2) )\n", "g[4,4] = tmp.simplify_full()\n", "g.display()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\begin{array}{lcl} g_{ \\, t \\, t }^{ \\phantom{\\, t}\\phantom{\\, t} } & = & \\frac{{\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} {\\mu}^{4} - a^{2} b^{2} {\\ell}^{2} - {\\ell}^{2} r^{4} - {\\left({\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} - {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} - b^{2} + 2 \\, m}{{\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2} + r^{2}} \\\\ g_{ \\, t \\, {\\phi} }^{ \\phantom{\\, t}\\phantom{\\, {\\phi}} } & = & -\\frac{a^{3} b^{2} {\\ell}^{2} - {\\left(a^{5} - a^{3} b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a {\\ell}^{2} {\\mu}^{2} - a {\\ell}^{2}\\right)} r^{4} + {\\left({\\left(a^{5} - 2 \\, a^{3} b^{2}\\right)} {\\ell}^{2} + 2 \\, a m\\right)} {\\mu}^{2} - {\\left(2 \\, a b^{2} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{3} - a b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - {\\left(a^{3} + a b^{2}\\right)} {\\ell}^{2}\\right)} r^{2} - 2 \\, a m}{a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\\\ g_{ \\, t \\, {\\psi} }^{ \\phantom{\\, t}\\phantom{\\, {\\psi}} } & = & -\\frac{b {\\ell}^{2} {\\mu}^{2} r^{4} + {\\left(a^{2} b^{3} - b^{5}\\right)} {\\ell}^{2} {\\mu}^{4} + {\\left(b^{5} {\\ell}^{2} - 2 \\, b m\\right)} {\\mu}^{2} + {\\left(2 \\, b^{3} {\\ell}^{2} {\\mu}^{2} + {\\left(a^{2} b - b^{3}\\right)} {\\ell}^{2} {\\mu}^{4}\\right)} r^{2}}{b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(b^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\\\ g_{ \\, r \\, r }^{ \\phantom{\\, r}\\phantom{\\, r} } & = & \\frac{r^{4} + {\\left({\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2}\\right)} r^{2}}{{\\ell}^{2} r^{6} + {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{4} + a^{2} b^{2} + {\\left(a^{2} b^{2} {\\ell}^{2} + a^{2} + b^{2} - 2 \\, m\\right)} r^{2}} \\\\ g_{ \\, {\\mu} \\, {\\mu} }^{ \\phantom{\\, {\\mu}}\\phantom{\\, {\\mu}} } & = & \\frac{{\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2} + r^{2}}{{\\left(a^{2} - b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - b^{2} {\\ell}^{2} - {\\left({\\left(a^{2} - 2 \\, b^{2}\\right)} {\\ell}^{2} + 1\\right)} {\\mu}^{2} + 1} \\\\ g_{ \\, {\\phi} \\, {\\phi} }^{ \\phantom{\\, {\\phi}}\\phantom{\\, {\\phi}} } & = & -\\frac{a^{4} b^{2} {\\ell}^{2} + {\\left(a^{4} - a^{2} b^{2} - {\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{2} - 2 \\, a^{2} m\\right)} {\\mu}^{4} + {\\left(a^{2} {\\ell}^{2} - {\\left(a^{2} {\\ell}^{2} - 1\\right)} {\\mu}^{2} - 1\\right)} r^{4} - a^{2} b^{2} - 2 \\, a^{2} m - {\\left(a^{4} - 2 \\, a^{2} b^{2} - {\\left(a^{6} - 2 \\, a^{4} b^{2}\\right)} {\\ell}^{2} - 4 \\, a^{2} m\\right)} {\\mu}^{2} - {\\left({\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{4} - {\\left(a^{4} + a^{2} b^{2}\\right)} {\\ell}^{2} + 2 \\, {\\left(a^{2} b^{2} {\\ell}^{2} - b^{2}\\right)} {\\mu}^{2} + a^{2} + b^{2}\\right)} r^{2}}{a^{4} b^{2} {\\ell}^{4} - 2 \\, a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{4} {\\ell}^{4} - 2 \\, a^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\\\ g_{ \\, {\\phi} \\, {\\psi} }^{ \\phantom{\\, {\\phi}}\\phantom{\\, {\\psi}} } & = & -\\frac{2 \\, {\\left(a b m {\\mu}^{4} - a b m {\\mu}^{2}\\right)}}{a^{2} b^{4} {\\ell}^{4} - {\\left(a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + {\\left({\\left(a^{4} b^{2} - a^{2} b^{4}\\right)} {\\ell}^{4} - {\\left(a^{4} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} b^{2} {\\ell}^{4} - {\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\\\ g_{ \\, {\\psi} \\, {\\psi} }^{ \\phantom{\\, {\\psi}}\\phantom{\\, {\\psi}} } & = & -\\frac{{\\left(b^{2} {\\ell}^{2} - 1\\right)} {\\mu}^{2} r^{4} - {\\left(a^{2} b^{2} - b^{4} - {\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{2} + 2 \\, b^{2} m\\right)} {\\mu}^{4} + {\\left(b^{6} {\\ell}^{2} - b^{4}\\right)} {\\mu}^{2} + {\\left({\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{4} + 2 \\, {\\left(b^{4} {\\ell}^{2} - b^{2}\\right)} {\\mu}^{2}\\right)} r^{2}}{b^{6} {\\ell}^{4} - 2 \\, b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(b^{4} {\\ell}^{4} - 2 \\, b^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\end{array}$ " ], "text/plain": [ "g_t,t = ((a^4 - 2*a^2*b^2 + b^4)*l^2*mu^4 - a^2*b^2*l^2 - l^2*r^4 - ((a^4 - 2*a^2*b^2 + b^4)*l^2 + a^2 - b^2)*mu^2 - ((a^2 + b^2)*l^2 + 1)*r^2 - b^2 + 2*m)/((a^2 - b^2)*mu^2 + b^2 + r^2) \n", "g_t,ph = -(a^3*b^2*l^2 - (a^5 - a^3*b^2)*l^2*mu^4 - (a*l^2*mu^2 - a*l^2)*r^4 + ((a^5 - 2*a^3*b^2)*l^2 + 2*a*m)*mu^2 - (2*a*b^2*l^2*mu^2 + (a^3 - a*b^2)*l^2*mu^4 - (a^3 + a*b^2)*l^2)*r^2 - 2*a*m)/(a^2*b^2*l^2 + ((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^2 + (a^2*l^2 - 1)*r^2 - b^2) \n", "g_t,ps = -(b*l^2*mu^2*r^4 + (a^2*b^3 - b^5)*l^2*mu^4 + (b^5*l^2 - 2*b*m)*mu^2 + (2*b^3*l^2*mu^2 + (a^2*b - b^3)*l^2*mu^4)*r^2)/(b^4*l^2 + ((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^2 + (b^2*l^2 - 1)*r^2 - b^2) \n", "g_r,r = (r^4 + ((a^2 - b^2)*mu^2 + b^2)*r^2)/(l^2*r^6 + ((a^2 + b^2)*l^2 + 1)*r^4 + a^2*b^2 + (a^2*b^2*l^2 + a^2 + b^2 - 2*m)*r^2) \n", "g_mu,mu = ((a^2 - b^2)*mu^2 + b^2 + r^2)/((a^2 - b^2)*l^2*mu^4 - b^2*l^2 - ((a^2 - 2*b^2)*l^2 + 1)*mu^2 + 1) \n", "g_ph,ph = -(a^4*b^2*l^2 + (a^4 - a^2*b^2 - (a^6 - a^4*b^2)*l^2 - 2*a^2*m)*mu^4 + (a^2*l^2 - (a^2*l^2 - 1)*mu^2 - 1)*r^4 - a^2*b^2 - 2*a^2*m - (a^4 - 2*a^2*b^2 - (a^6 - 2*a^4*b^2)*l^2 - 4*a^2*m)*mu^2 - (((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^4 - (a^4 + a^2*b^2)*l^2 + 2*(a^2*b^2*l^2 - b^2)*mu^2 + a^2 + b^2)*r^2)/(a^4*b^2*l^4 - 2*a^2*b^2*l^2 + ((a^6 - a^4*b^2)*l^4 - 2*(a^4 - a^2*b^2)*l^2 + a^2 - b^2)*mu^2 + (a^4*l^4 - 2*a^2*l^2 + 1)*r^2 + b^2) \n", "g_ph,ps = -2*(a*b*m*mu^4 - a*b*m*mu^2)/(a^2*b^4*l^4 - (a^2*b^2 + b^4)*l^2 + ((a^4*b^2 - a^2*b^4)*l^4 - (a^4 - b^4)*l^2 + a^2 - b^2)*mu^2 + (a^2*b^2*l^4 - (a^2 + b^2)*l^2 + 1)*r^2 + b^2) \n", "g_ps,ps = -((b^2*l^2 - 1)*mu^2*r^4 - (a^2*b^2 - b^4 - (a^2*b^4 - b^6)*l^2 + 2*b^2*m)*mu^4 + (b^6*l^2 - b^4)*mu^2 + (((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^4 + 2*(b^4*l^2 - b^2)*mu^2)*r^2)/(b^6*l^4 - 2*b^4*l^2 + ((a^2*b^4 - b^6)*l^4 - 2*(a^2*b^2 - b^4)*l^2 + a^2 - b^2)*mu^2 + (b^4*l^4 - 2*b^2*l^2 + 1)*r^2 + b^2) " ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "g.display_comp(only_nonredundant=True)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "## Einstein equation\n", "\n", "The Ricci tensor of $g$ is" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Field of symmetric bilinear forms Ric(g) on the 5-dimensional differentiable manifold M\n" ] } ], "source": [ "Ric = g.ricci()\n", "print(Ric)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\begin{array}{lcl} \\mathrm{Ric}\\left(g\\right)_{ \\, t \\, t }^{ \\phantom{\\, t}\\phantom{\\, t} } & = & -\\frac{4 \\, {\\left({\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{4} {\\mu}^{4} - a^{2} b^{2} {\\ell}^{4} - {\\ell}^{4} r^{4} - b^{2} {\\ell}^{2} + 2 \\, {\\ell}^{2} m - {\\left({\\left(a^{4} - 2 \\, a^{2} b^{2} + b^{4}\\right)} {\\ell}^{4} + {\\left(a^{2} - b^{2}\\right)} {\\ell}^{2}\\right)} {\\mu}^{2} - {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{4} + {\\ell}^{2}\\right)} r^{2}\\right)}}{{\\left(a^{2} - b^{2}\\right)} {\\mu}^{2} + b^{2} + r^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, t \\, {\\phi} }^{ \\phantom{\\, t}\\phantom{\\, {\\phi}} } & = & \\frac{4 \\, {\\left(a^{3} b^{2} {\\ell}^{4} - {\\left(a^{5} - a^{3} b^{2}\\right)} {\\ell}^{4} {\\mu}^{4} - {\\left(a {\\ell}^{4} {\\mu}^{2} - a {\\ell}^{4}\\right)} r^{4} - 2 \\, a {\\ell}^{2} m + {\\left({\\left(a^{5} - 2 \\, a^{3} b^{2}\\right)} {\\ell}^{4} + 2 \\, a {\\ell}^{2} m\\right)} {\\mu}^{2} - {\\left(2 \\, a b^{2} {\\ell}^{4} {\\mu}^{2} + {\\left(a^{3} - a b^{2}\\right)} {\\ell}^{4} {\\mu}^{4} - {\\left(a^{3} + a b^{2}\\right)} {\\ell}^{4}\\right)} r^{2}\\right)}}{a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, t \\, {\\psi} }^{ \\phantom{\\, t}\\phantom{\\, {\\psi}} } & = & \\frac{4 \\, {\\left(b {\\ell}^{4} {\\mu}^{2} r^{4} + {\\left(a^{2} b^{3} - b^{5}\\right)} {\\ell}^{4} {\\mu}^{4} + {\\left(b^{5} {\\ell}^{4} - 2 \\, b {\\ell}^{2} m\\right)} {\\mu}^{2} + {\\left(2 \\, b^{3} {\\ell}^{4} {\\mu}^{2} + {\\left(a^{2} b - b^{3}\\right)} {\\ell}^{4} {\\mu}^{4}\\right)} r^{2}\\right)}}{b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} - a^{2} + b^{2}\\right)} {\\mu}^{2} + {\\left(b^{2} {\\ell}^{2} - 1\\right)} r^{2} - b^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, r \\, r }^{ \\phantom{\\, r}\\phantom{\\, r} } & = & -\\frac{4 \\, {\\left({\\ell}^{2} r^{4} + {\\left({\\left(a^{2} - b^{2}\\right)} {\\ell}^{2} {\\mu}^{2} + b^{2} {\\ell}^{2}\\right)} r^{2}\\right)}}{{\\ell}^{2} r^{6} + {\\left({\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{4} + a^{2} b^{2} + {\\left(a^{2} b^{2} {\\ell}^{2} + a^{2} + b^{2} - 2 \\, m\\right)} r^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, {\\mu} \\, {\\mu} }^{ \\phantom{\\, {\\mu}}\\phantom{\\, {\\mu}} } & = & -\\frac{4 \\, {\\left({\\left(a^{2} - b^{2}\\right)} {\\ell}^{2} {\\mu}^{2} + b^{2} {\\ell}^{2} + {\\ell}^{2} r^{2}\\right)}}{{\\left(a^{2} - b^{2}\\right)} {\\ell}^{2} {\\mu}^{4} - b^{2} {\\ell}^{2} - {\\left({\\left(a^{2} - 2 \\, b^{2}\\right)} {\\ell}^{2} + 1\\right)} {\\mu}^{2} + 1} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, {\\phi} \\, {\\phi} }^{ \\phantom{\\, {\\phi}}\\phantom{\\, {\\phi}} } & = & \\frac{4 \\, {\\left(a^{4} b^{2} {\\ell}^{4} - a^{2} b^{2} {\\ell}^{2} - 2 \\, a^{2} {\\ell}^{2} m - {\\left({\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{4} + 2 \\, a^{2} {\\ell}^{2} m - {\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2}\\right)} {\\mu}^{4} + {\\left(a^{2} {\\ell}^{4} - {\\left(a^{2} {\\ell}^{4} - {\\ell}^{2}\\right)} {\\mu}^{2} - {\\ell}^{2}\\right)} r^{4} + {\\left({\\left(a^{6} - 2 \\, a^{4} b^{2}\\right)} {\\ell}^{4} + 4 \\, a^{2} {\\ell}^{2} m - {\\left(a^{4} - 2 \\, a^{2} b^{2}\\right)} {\\ell}^{2}\\right)} {\\mu}^{2} + {\\left({\\left(a^{4} + a^{2} b^{2}\\right)} {\\ell}^{4} - {\\left({\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{4} - {\\left(a^{2} - b^{2}\\right)} {\\ell}^{2}\\right)} {\\mu}^{4} - {\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} - 2 \\, {\\left(a^{2} b^{2} {\\ell}^{4} - b^{2} {\\ell}^{2}\\right)} {\\mu}^{2}\\right)} r^{2}\\right)}}{a^{4} b^{2} {\\ell}^{4} - 2 \\, a^{2} b^{2} {\\ell}^{2} + {\\left({\\left(a^{6} - a^{4} b^{2}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{4} - a^{2} b^{2}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{4} {\\ell}^{4} - 2 \\, a^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, {\\phi} \\, {\\psi} }^{ \\phantom{\\, {\\phi}}\\phantom{\\, {\\psi}} } & = & \\frac{8 \\, {\\left(a b {\\ell}^{2} m {\\mu}^{4} - a b {\\ell}^{2} m {\\mu}^{2}\\right)}}{a^{2} b^{4} {\\ell}^{4} - {\\left(a^{2} b^{2} + b^{4}\\right)} {\\ell}^{2} + {\\left({\\left(a^{4} b^{2} - a^{2} b^{4}\\right)} {\\ell}^{4} - {\\left(a^{4} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(a^{2} b^{2} {\\ell}^{4} - {\\left(a^{2} + b^{2}\\right)} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\\\ \\mathrm{Ric}\\left(g\\right)_{ \\, {\\psi} \\, {\\psi} }^{ \\phantom{\\, {\\psi}}\\phantom{\\, {\\psi}} } & = & \\frac{4 \\, {\\left({\\left(b^{2} {\\ell}^{4} - {\\ell}^{2}\\right)} {\\mu}^{2} r^{4} + {\\left({\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{4} - 2 \\, b^{2} {\\ell}^{2} m - {\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2}\\right)} {\\mu}^{4} + {\\left(b^{6} {\\ell}^{4} - b^{4} {\\ell}^{2}\\right)} {\\mu}^{2} + {\\left({\\left({\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{4} - {\\left(a^{2} - b^{2}\\right)} {\\ell}^{2}\\right)} {\\mu}^{4} + 2 \\, {\\left(b^{4} {\\ell}^{4} - b^{2} {\\ell}^{2}\\right)} {\\mu}^{2}\\right)} r^{2}\\right)}}{b^{6} {\\ell}^{4} - 2 \\, b^{4} {\\ell}^{2} + {\\left({\\left(a^{2} b^{4} - b^{6}\\right)} {\\ell}^{4} - 2 \\, {\\left(a^{2} b^{2} - b^{4}\\right)} {\\ell}^{2} + a^{2} - b^{2}\\right)} {\\mu}^{2} + {\\left(b^{4} {\\ell}^{4} - 2 \\, b^{2} {\\ell}^{2} + 1\\right)} r^{2} + b^{2}} \\end{array}$ " ], "text/plain": [ "Ric(g)_t,t = -4*((a^4 - 2*a^2*b^2 + b^4)*l^4*mu^4 - a^2*b^2*l^4 - l^4*r^4 - b^2*l^2 + 2*l^2*m - ((a^4 - 2*a^2*b^2 + b^4)*l^4 + (a^2 - b^2)*l^2)*mu^2 - ((a^2 + b^2)*l^4 + l^2)*r^2)/((a^2 - b^2)*mu^2 + b^2 + r^2) \n", "Ric(g)_t,ph = 4*(a^3*b^2*l^4 - (a^5 - a^3*b^2)*l^4*mu^4 - (a*l^4*mu^2 - a*l^4)*r^4 - 2*a*l^2*m + ((a^5 - 2*a^3*b^2)*l^4 + 2*a*l^2*m)*mu^2 - (2*a*b^2*l^4*mu^2 + (a^3 - a*b^2)*l^4*mu^4 - (a^3 + a*b^2)*l^4)*r^2)/(a^2*b^2*l^2 + ((a^4 - a^2*b^2)*l^2 - a^2 + b^2)*mu^2 + (a^2*l^2 - 1)*r^2 - b^2) \n", "Ric(g)_t,ps = 4*(b*l^4*mu^2*r^4 + (a^2*b^3 - b^5)*l^4*mu^4 + (b^5*l^4 - 2*b*l^2*m)*mu^2 + (2*b^3*l^4*mu^2 + (a^2*b - b^3)*l^4*mu^4)*r^2)/(b^4*l^2 + ((a^2*b^2 - b^4)*l^2 - a^2 + b^2)*mu^2 + (b^2*l^2 - 1)*r^2 - b^2) \n", "Ric(g)_r,r = -4*(l^2*r^4 + ((a^2 - b^2)*l^2*mu^2 + b^2*l^2)*r^2)/(l^2*r^6 + ((a^2 + b^2)*l^2 + 1)*r^4 + a^2*b^2 + (a^2*b^2*l^2 + a^2 + b^2 - 2*m)*r^2) \n", "Ric(g)_mu,mu = -4*((a^2 - b^2)*l^2*mu^2 + b^2*l^2 + l^2*r^2)/((a^2 - b^2)*l^2*mu^4 - b^2*l^2 - ((a^2 - 2*b^2)*l^2 + 1)*mu^2 + 1) \n", "Ric(g)_ph,ph = 4*(a^4*b^2*l^4 - a^2*b^2*l^2 - 2*a^2*l^2*m - ((a^6 - a^4*b^2)*l^4 + 2*a^2*l^2*m - (a^4 - a^2*b^2)*l^2)*mu^4 + (a^2*l^4 - (a^2*l^4 - l^2)*mu^2 - l^2)*r^4 + ((a^6 - 2*a^4*b^2)*l^4 + 4*a^2*l^2*m - (a^4 - 2*a^2*b^2)*l^2)*mu^2 + ((a^4 + a^2*b^2)*l^4 - ((a^4 - a^2*b^2)*l^4 - (a^2 - b^2)*l^2)*mu^4 - (a^2 + b^2)*l^2 - 2*(a^2*b^2*l^4 - b^2*l^2)*mu^2)*r^2)/(a^4*b^2*l^4 - 2*a^2*b^2*l^2 + ((a^6 - a^4*b^2)*l^4 - 2*(a^4 - a^2*b^2)*l^2 + a^2 - b^2)*mu^2 + (a^4*l^4 - 2*a^2*l^2 + 1)*r^2 + b^2) \n", "Ric(g)_ph,ps = 8*(a*b*l^2*m*mu^4 - a*b*l^2*m*mu^2)/(a^2*b^4*l^4 - (a^2*b^2 + b^4)*l^2 + ((a^4*b^2 - a^2*b^4)*l^4 - (a^4 - b^4)*l^2 + a^2 - b^2)*mu^2 + (a^2*b^2*l^4 - (a^2 + b^2)*l^2 + 1)*r^2 + b^2) \n", "Ric(g)_ps,ps = 4*((b^2*l^4 - l^2)*mu^2*r^4 + ((a^2*b^4 - b^6)*l^4 - 2*b^2*l^2*m - (a^2*b^2 - b^4)*l^2)*mu^4 + (b^6*l^4 - b^4*l^2)*mu^2 + (((a^2*b^2 - b^4)*l^4 - (a^2 - b^2)*l^2)*mu^4 + 2*(b^4*l^4 - b^2*l^2)*mu^2)*r^2)/(b^6*l^4 - 2*b^4*l^2 + ((a^2*b^4 - b^6)*l^4 - 2*(a^2*b^2 - b^4)*l^2 + a^2 - b^2)*mu^2 + (b^4*l^4 - 2*b^2*l^2 + 1)*r^2 + b^2) " ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Ric.display_comp(only_nonredundant=True)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Let us check that $g$ is a solution of the vacuum Einstein equation with the cosmological constant $\\Lambda = - 6 \\ell^2$:" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "text/html": [ " $\\newcommand{\\Bold}[1]{\\mathbf{#1}}\\mathrm{True}$ " ], "text/plain": [ "True" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Lambda = -6*l^2\n", "Ric == 2/3*Lambda*g" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "SageMath 8.0.beta12", "language": "", "name": "sagemath" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 0 }