Typically when solving a Hamiltonian of ye olde form $ \frac{1}{2m} (\bar P - \frac{q}{c} \bar A)^2 $, you do separation of variables.
For simplicity say that $ A = - B Y \hat x $.
You can rewrite it as $\frac{1}{2m} ( (\hbar k_x + \frac{q}{c} B y )^2 + P_y^2 + \hbar^2 k_z^2 ) $.
Now, we're putting everything in the position basis by applying $ \langle x | \otimes \langle y | \otimes \langle z |$ on the left side, and then theoretically some sort of $ |\psi\rangle $ in 3 dimensions on the right side, and then computing the results.
However, then I'm slightly confused about what happens formally with the equation for the y, $ (\hbar k_x + \frac{qB}{c} y)^2 + P_y^2 $. You can solve this as a harmonic oscillator and you'll get solutions for $ \psi(y)$.
Now say, for instance, say you wanted to compute $ \langle Y \otimes I \otimes I\rangle $. In theory, you'd miss it altogether if you apply $ | \psi(x) \rangle \otimes | \psi(y) \rangle \otimes |\psi(z) \rangle $ since Y will only operate on $ \psi(x)$ which is just a plane wave. Frankly I'm not even sure how this computation would work/if it'll just pass through or something. You'd have to apply something like $ | \psi(y) \rangle \otimes | \psi(x) \rangle... $ and so on.
But this shouldn't occur because you should be able to somehow measure the average "Y" value to be $ \frac{\hbar k_x c}{qB} $. This also leads to bad outcomes because if you decide to compute $ \langle dX/dt \rangle $, you'd compute the commutator of $ X \otimes I \otimes I $ with the Hamiltonian; it'll commute with anything which is the identity in the X coordinate; $ D_t \langle X \rangle $ will end up being proportional to the expectation value of Y, and I know this quantity shouldn't end up being trivial or something.
Can someone point out the misstep here in how I propose to find the expectation value for Y?