The problem of quantum charged particles in the presence of electric and magnetic fields has always remained on target of any theoretical physicist who wishes to predict the behavior of one or a swarm of such particles imbedded in the structure of a solid. The exact solution in three dimensions has eluded these attempts for years. However, when the particles remain in a two dimensional solid, the problem of the exact solution becomes both feasible and possible. Experimentally, the last situation has only been possible to be realized until very recently in the laboratory thanks to the technological means that allow the experimentalists to produce quasi dimensional samples of metals and semi- metals. One can induce large magnetic behavior combining them with ferromag- netic 2D materials such as and. Indeed we should also mention the recent discovery of Magnetic Graphene which falls in the class of these two dimensional samples when combined with the aforemen- tioned substrates. The propagator can also be used in three dimensions with some little arrangements. However, I restrict myself to the two dimensional case given the importance of two dimensional solids nowadays.
However, the stationary solution of a bounded Hamiltonian is only of academic interest as we would like to know the evolution of the particle inside of the two dimensional solid. This goal can only be achieved if we know the propagator of the particle(s) along the time and the properties that may or may not change in the sample: phase transitions from metal to semimetal or to an insulator phase, the behavior of the Anderson localization, the magneto-optics properties of the solid and the like. This wealth of physical effects can be studied if we know the Feynman propagator of the system and the help of computing algebraic calculus. The aim of this paper is to show that these results can be attained with little effort and a touch of elegance.
The problem of the exact solution of the Feynman propagator of the system simplifies a great deal with the help of the Stationary Phase Approximation that turns to be exact in the case of a quadratic Hamiltonians. In this paper, we use the most general two-dimensional Hamiltonian with magnetic effects provided by a magnetic field pointing in the perpendicular direction of the plane. The rest of the interaction is represented by anisotropic harmonic oscillators lying in the plane along perpendicular directions.
2. Quantum Mechanical Hamiltonian
The Hamiltonian described in the Introduction has the form:
where and are the first order frequencies on each direction of the plane. The set of variables are well known in electromagnetic theory. One can select the following gauge for the magnetic vector that produces a constant magnetic field of arbitrary strength and perpendicular to the plane:
and the Hamiltonian reads:
Obviously is the strength of the constant magnetic field. Let us now call:
and we scale the position and momentum operators in a scale free way in the form:
Again, with this change of variables the Hamiltonian takes the form:
that looks as a quadratic form that we can easily diagonalize:
where and are constants:
The non vanishing real roots of the quartic equation take the form:
and multiplying the Hamiltonian by, one finally obtains the following eigenvalues:
or in a more tractable manner:
If one diagonalizes the isotropic harmonic oscillator in both cartesian and angular coordinates the Landau eigenlevels are labeled respectively by the quantum numbers and. Both eigenlevels take the form:
3. Lagrangian Formalism and Construction of the Propagator
The Lagrangian of the system has the obvious form:
and now select a different gauge, but representing also a constant magnetic field in the perpendicular direction to the plane:
The change of gauge causes no problems in the final solution of the propagator as it was pointed out many years ago by Bialynicki-Birula  . The propagator and quantum amplitudes in general are invariant under a gauge transformation. In our case:
It is easy to find out the actual fom of, obviously independent of. One finds:
which vanishes for the isotropic limit. For a throughout and recent discussion on gauge independence of the propagator see reference  .
The Lagrangian looks now as:
and the Equations of Motion are:
where where and are the first order frequencies on each direction of the plane due to the potential of the net. The real solutions of these Equations are always harmonic functions that we can write of the form:
We shall name:, to the points at time equal zero and, , those points at an arbitrary evolution time. Thus, one can see trivially that the following relationships hold:
Next we shall be concerned with the evaluation of the classical action of the anisotropic net in the presence of a uniform and constant magnetic field. The various steps of the calculation are listed in Section 4 at the end of the paper. For the moment we shall be concerned formally with the action derived from this Lagrangian having the form:
where we have used (16). Next we integrate by parts the terms and leading easily to:
The integral of the right hand side is trivially zero as it contains only the equations of motion (17). This property has been used in other contexts in reference  for the general case of kinetic energy plus an arbitrary potential depending just upon the coordinates. One finally founds:
and substituting the values of the coordinates and their derivatives one finally obtains the following expression for the action:
The Amplitude of the Propagator is now calculated with the help of the VanVleck-DeWitt-Morette Determinant. As it is well known the form of the Determinant is:
where is the number of spatial dimensions of the problem. In our case. Furthermore the form of the Determinant is:
The result leads to:
Substituting we obtain (with):
and the propagator takes the final form:
The expressions for and can be found in Section 4 but now we have substituted everywhere as in the expression. As a final word one can say that these seemingly tedious calculations simplify greatly when one uses any of the available packages of algebraic calculations such as MATHEMATICAT or MAPLET. Also limits in the case of vanishing magnetic fields or isotropic harmonic oscillator (i.e.) have been discussed in  . One of the main goals of this work has been precisely to put at work these software applications to the non trivial field of path integrals and its use in material science and quantum optics: for instance propagating a gaussian wave packet in two dimensional magnetic solids when the gaussian represents charged particles or intense femtosecond light pulses. A particular case is the one of a “Schrödinger cat” propagating in the net. The wave function can be simulated as two widely separated gaussians charged with “-2e-charge” as in a Cooper pair  .
where in atomic units containing the electric charge and is the
variance of each packet. Although these projects for Master Thesis or Graduate Courses can be considered as a little cumbersome, they might attract the interest of the young physicists to the condensed matter field in which much of the research is being done nowadays.
4. The Main Result and Calculations with Useful Definitions of Interest
To obtain in Equations (23)-(25) as functions of. I have used for the calculation of the inverse matrix with the help of the MATHEMATICAT package:
where the matrix elements and the determinant of (24) are listed below:
and the determinant takes the form:
Likewise one can easily calculate
The derivatives then take the form:
Notice the interesting relationship:
The following constants have been defined in the body of the Paper as:
With these definitions the following expressions hold:
We have presented and fully solved the propagator of the anisotropic two dimensional harmonic oscillator in the presence of a constant magnetic field in the perpendicular direction of the plane. The solution is based on the Stationary Phase Approximation that is exact in this case as the Hamiltonian contains terms no more than quadratic terms. We believe this solution is most interesting today when a large amount of plane samples as grapheme with large mobility (large number of electrons per unit of area) are submitted to various both strong and constant magnetic fields.
This research has been carried out with the partial support of MINECO under the project MAT2013-46308-C2-1-R and JCyL Project Number SA226U13.