The Kepler problem on the lattice

Diego Sanjinés,Evaristo Mamani,Javier Velasco
2024-06-27
Abstract:We study the motion of a particle in a 3-dimensional lattice in the presence of a Coulomb potential, but we demonstrate semiclassicaly that the trajectories will always remain in a plane which can be taken as a rectangular lattice. The Hamiltonian model for this problem is the conservative tight-binding one with lattice constants a, b and hopping elements A, B in the XY axes, respectively. We use the semiclassical and quantum formalisms; for the latter we apply the pseudo-spectral algorithm to integrate the Schroedinger equation. Since the lattice discrete subspace is not isotropic, the angular momentum is not conserved, which has interesting consequences as chaotic trajectories and precession trajectories, similar to the astronomical precession trajectories due to non-central gravitational forces, notably, the non-relativistic Mercury's perihelion precession. Although the elements of the mass tensor are naturally different in a rectangular lattice, these can be chosen to be still different in the continuum, which permits to study the motion with the usual Newtonian kinetic energies. We calculate also the contour plots of an initial Gaussian wavepacket as it moves in the lattice and we propose an "intrinsec angular momentum" associated to its asymmetrical deformation, such that the quantum and semiclassical angular momenta could be simply related.
Quantum Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the motion of charged particles (such as electrons) on the lattice in the Coulomb electric field generated by another fixed charged particle (such as donor impurity). Specifically, although the electric field is three - dimensional, research shows that the motion of particles will be restricted to a plane. The paper explores this problem by using the semiclassical model (Hamiltonian equations) and quantum formalism (tight - binding model and pseudo - spectral algorithm to solve the Schrödinger equation). The main concerns include: 1. **Energy Conservation and Non - Conservation of Angular Momentum**: Due to the lack of rotational symmetry in the discrete subspace of the lattice, the total energy is conserved, but the angular momentum is not. 2. **Educational Significance**: The author believes that these topics are very suitable for teaching in undergraduate courses, such as classical mechanics, quantum mechanics and solid - state physics. 3. **Continuous Limit of Rectangular Lattice**: Through the continuous limit of the rectangular lattice, the mass differences in different directions can be explained, so as to understand the anisotropy of particle trajectories. 4. **Orbital Precession Effect**: In the square lattice, due to the non - conservation of angular momentum, the orbital precession effect can be observed, which is similar to the precession of the perihelion of Mercury in astronomy due to the non - central gravitational field. ### Specific Results - **Reproducing Cortés' Results**: The paper reproduces Cortés' research results on the trajectories of particles in the attractive central field. Cortés assumed that the masses of particles on the X - axis and Y - axis are different, and this paper explains this through the continuous limit of the rectangular lattice. - **Orbital Precession Effect**: In the square lattice, due to the non - conservation of angular momentum, the orbital precession effect is observed. This is similar to the orbital precession caused by the non - central gravitational field in astronomy. ### Mathematical Formulas - **Hamiltonian Function**: \[ H(r, k)=2A(1 - \cos(ak_x))+2B(1 - \cos(bk_y))+2C(1 - \cos(ck_z))-\frac{V_1}{\sqrt{x^2 + y^2+z^2}} \] - **Hamiltonian Equations**: \[ \dot{x}=\frac{\partial H}{\partial p_x}=\frac{2Aa}{\hbar}\sin(ak_x) \] \[ \dot{y}=\frac{\partial H}{\partial p_y}=\frac{2Bb}{\hbar}\sin(bk_y) \] \[ \dot{z}=\frac{\partial H}{\partial p_z}=\frac{2Cc}{\hbar}\sin(ck_z) \] \[ \hbar\dot{k}_x =-\frac{\partial H}{\partial x}=\frac{V_1x}{(x^2 + y^2+z^2)^{3/2}} \] \[ \hbar\dot{k}_y =-\frac{\partial H}{\partial y}=\frac{V_1y}{(x^2 + y^2+z^2)^{3/2}} \] \[ \hbar\dot{k}_z =-\frac{\partial H}{\partial z}=\frac{V_1z}{(x^2 + y^2+z^2)^{3/2}} \] - **Quantum Angular Momentum**: \[ L_q = |z\times k|=z_xk_y - z_yk_x \] - **Semiclassical Angular Momentum**: \[ L_c = |r\times k|=xk_y - yk_x \] - **Intrinsic Angular Momentum**: \[ S = |s\times k| \] where \(s=(s_x, s_y)\), and \(s_x=\sum_{n,m}(an - z_x)\)