Kerr_sign_gpp.ipynb

Kernel: SageMath 7.1

Carter time machine

%display latex

var('a r')

var('th', latex_name=r'\theta')

f(r,a,th) = (r^2+a^2)*(r^2+a^2*cos(th)^2) + 2*a^2*r*sin(th)^2
f

g0 = plot(f(r,0.9,0), (r,-1.5,1.5), legend_label=r'$\theta=0$', 
          thickness=2, linestyle=':', color='red')
g1 = plot(f(r,0.9,pi/4), (r,-1.5,1.5), legend_label=r'$\theta=\pi/4$', 
          thickness=2, linestyle='-.', color='grey')
g2 = plot(f(r,0.9,pi/3), (r,-1.5,1.5), legend_label=r'$\theta=\pi/3$', 
          thickness=2, linestyle='--', color='blue')
g3 = plot(f(r,0.9,pi/2), (r,-1.5,1.5), legend_label=r'$\theta=\pi/2$', 
          thickness=2, color='violet')
graph = g0+g1+g2+g3
graph.axes_labels([r'$r/m$', r'$\rho^2 (r^2+a^2) + 2 a^2 m r \, \sin^2\theta$'])
graph.set_legend_options(loc='upper right')
graph

graph.save('ker_sign_gpp.pdf')

rp(a) = 1 + sqrt(1-a^2)
rm(a) = 1 - sqrt(1-a^2)

rp

rm

rp(0.9)

rm(0.9)

df = diff(f(r,a,th), r).simplify_full()
df

s = solve(df==0, r, solution_dict=True)
s

rmin = s[2][r]
rmin

df.subs(r=rmin).simplify_full()

plot(rmin.subs(th=pi/2), (a,0.01, 0.9))

fmin = f(rmin.subs(th=pi/2), a, pi/2).simplify_full()
fmin