Show that the energy of a particle in a constant magnetic field oriented along the z axis can be written as

where and  are creation and anihilation operators.

Main reference: - Stephani, H., Kramer, D., MacCallum, M., Hoenselaers, C. Herlt, E.  Exact Solutions of Einstein's Field Equations, Cambridge Monographs on Mathematical Physics, second edition. Cambridge University Press, 2003.

The authors reviewed more than 4,000 papers containing solutions to Einstein’s equations in the literature and organized the material into chapters according to the physical properties of these solutions.

These solutions are now digitized within Maple 2016, so that it is now possible to actually compute with them.

Given the spacetime metric,

a) Compute the Ricci and Weyl scalars

b) Compute the trace of

where  is some function of the radial coordinate,  is the Ricci tensor,  is the covariant derivative operator and  is the stress-energy tensor

c) Compute the components of the traceless part of   of item b)

d) Compute an exact solution to the nonlinear system of differential equations conformed by the components of   obtained in c)

Background: paper from February/2013, "Withholding Potentials, Absence of Ghosts and Relationship between Minimal Dilatonic Gravity and f(R) Theories", by P. Fiziev.

From the "What is new in Physics in Maple 2016" page:

Generally speaking a canonical form is obtained using transformations that leave invariant the tetrad metric in a tetrad system of references, so that theWeyl scalars are fixed as much as possible (conventionally, either equal to 0 or to 1).

Bringing a tetrad in canonical form is a relevant step in the tackling of the equivalence problem between two spacetime metrics.

The implementation is as in "General Relativity, an Einstein century survey", edited by S.W. Hawking (Cambridge) and W. Israel (U. Alberta, Canada), specifically Chapter 7 written by S. Chandrasekhar, page 388:

The transformations (rotations of the tetrad system of references) used are of Class I, II and III as defined in Chandrasekar's chapter - equations (7.79) in page 384, (7.83) and (7.84) in page 385. Transformations of Class I can be performed with the command Physics:-Tetrads:-TransformTetrad using the optional argument nullrotationwithfixedl_, of Class II using nullrotationwithfixedn_ and of Class III by calling TransformTetrad(spatialrotationsm_mb_plan, boostsn_l_plane), so with the two optional arguments simultaneously.

The determination of appropriate transformation parameters to be used in these rotations, as well as the sequence of transformations happens all automatically by using the optional argument, canonicalform of TransformTetrad .