A Classical model for diamagnetism

I'm trying to derive a classical model to describe diamagnetism, and I'd like to understand if it is possible to understand the basic properties of diamagnetic materials with it i.e.:

''A diamagnetic material placed in an external magnetic field $$B_{ext}$$ develops a magnetic dipole moment directed opposite $$B_{ext}$$. If the field is nonuniform, the diamagnetic material is repelled from a region of greater magnetic field toward a region of lesser field.''

or as much information as possible.

I started considering only the orbital motion of the electron in the atom. Consider an electron in orbit with angular velocity $$\omega$$. Its magnetic moment would be

\begin{align} \vec{\mu} = \frac{-e}{2m_e}\vec{L}_{orb} \end{align}

Subject to a field $${\bf H}$$, with induction, $$\mathbf{B} = \mu_0\mathbf{H}$$, this magnetic moment will experience a torque

\begin{align} \vec{\tau} &= \frac{d\vec{L}}{dt} \\ &= \vec{\mu} \times \vec{B} \\ &= \mu_0 \frac{-e}{2m_e} \vec{L} \times \vec{H} \\ &= \mu_0 \frac{e}{2m_e} \vec{H}\times\vec{L} \\ &= \vec{\omega} \times \vec{L} \end{align}

Now, it seems from the last relation that this magnetic moment is in precession. It should be this magnetic moment to cause the effects of diamagnetism, but I'm not sure how to treat this preceding magnetic moment and show that it is antiparallel to the external magnetic field and its repulsion to it.

The Hamiltonian $$H,$$ for a system of charged particles interacting via a potential energy $$U$$ is: $$H\left(\boldsymbol{r}_{k}, \boldsymbol{p}_{k}\right)=\sum_{j=1}^{N} \frac{\boldsymbol{p}_{j}^{2}}{2 m_{j}}+U\left(\boldsymbol{r}_{k}\right)$$ Particles defined with, mass $$m_{j},$$ position vector $$r_{j},$$ and momenta $$p_{j}$$. Under external magnetic field with vector potential $$A(r)$$ changes Hamiltonian to: $$H\left(\boldsymbol{r}_{k}, \boldsymbol{p}_{k}\right)=\sum_{j=1}^{N} \frac{\left[\boldsymbol{p}_{j}-\frac{e_{j}}{c} \boldsymbol{A}\left(\boldsymbol{r}_{j}\right)\right]^{2}}{2 m_{j}}+U\left(\boldsymbol{r}_{k}\right)$$
where $$e_{j}$$ are the charges of masses $$m_{j}$$. Using the above equation, statistical mechanics predicts that the energy-the thermal average of the Hamiltonian-does not depend on the external field. So, the system exhibits neither a paramagnetic nor a diamagnetic response.
• Note that Feynman's argument is patently false because an individual particle's energy absolutely does depend on $B$ (the relevant term in the energy is $p_i\cdot A(r_i)$). Once you integrate over momenta, however, that dependence becomes irrelevant. A correct argument is given in this answer physics.stackexchange.com/a/122529/8562. If Feynman's statement were true, the Lorentz force wouldn't exist at all. May 3 at 21:37