Numerical work with sage 7 – Eigenvalues of the volume operator in Loop quantum Gravity

This week I have been looking at how to actually  calculate the eigenvalues of the volume operator in Loop Quantum gravity. I’ve still got a lot more work to do on this but so far I can calculate the volume of a 4-valent spinnetwork node when the dimension d of the Hilbert space is 2. In order to generalize this I’ll need to generate appropriate Jacobi Matrices when the dimension d is different from 2, calculate the matrix elements and then I’ll need to find their eigenvalues. All this looks relatively straightforward given the work that has been done so far.  I’ll also look to either interactively enter j1, j2, j3, j4 or use a data file for data entry.

volume spectrum sage 7 fig 1

volume spectrum sage 7 fig 2

 

So here we see that the case (j1,j2,j3,j4) = (1/2, 1/2, 3/2, 3/2) has eigenvalue 0.463, which is known to be correct in the Loop Quantum Gravity case. I have also checked other 2 dimensional cases including (1/2,1/2,1/2,1/2) where the eigenvalue is 0.310.

Advertisements

9 thoughts on “Numerical work with sage 7 – Eigenvalues of the volume operator in Loop quantum Gravity”

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s