diracfdm.rar from Semiclassical approach for excitonic spectrum of Coulomb coupling between two Dirac particles

The properties of the energy spectrum of excitons in monolayer transition metal dichalcogenides are investigated using a multiband model. In the multiband model, we use the excitonic Hamiltonian in the product base of the Dirac single-particle states at the conduction and valence band edges. Following the separation of variables, we decouple the corresponding energy eigenvalue system of the first-order ODE radial equations rigorously and solve the resulting second-order ODE self-consistently, using the finite difference method, thus we determine the energy eigenvalues of the discrete excitonic spectrum and the corresponding wave functions. We also developed a WKB approach to solve the same spectral problem in semiclassical approximation for the resulting ODE. We compare the results for the energy spectrum and the corresponding eigen-function forms for WS₂ and WSe₂ obtained by means of both methods. We also compare our results for the energy spectrum with other theoretical works for excitons, and with available experimental data.