Isotropic Harmonic Oscillator Research Articles

Large-cage occupation and quantum dynamics of hydrogen molecules in sII clathrate hydrates.

Hydrogen clathrate hydrates are ice-like crystalline substances in which hydrogen molecules are trapped inside polyhedral cages formed by the water molecules. Small cages can host only a single H2 molecule, while each large cage can be occupied by up to four H2 molecules. Here, we present a neutron scattering study on the structure of the sII hydrogen clathrate hydrate and on the low-temperature dynamics of the hydrogen molecules trapped in its large cages, as a function of the gas content in the samples. We observe spectral features at low energy transfer (between 1 and 3meV), and we show that they can be successfully assigned to the rattling motion of a single hydrogen molecule occupying a large water cage. These inelastic bands remarkably lose their intensity with increasing the hydrogen filling, consistently with the fact that the probability of single occupation (as opposed to multiple occupation) increases as the hydrogen content in the sample gets lower. The spectral intensity of the H2 rattling bands is studied as a function of the momentum transfer for partially emptied samples and compared with three distinct quantum models for a single H2 molecule in a large cage: (i) the exact solution of the Schrödinger equation for a well-assessed semiempirical force field, (ii) a particle trapped in a rigid sphere, and (iii) an isotropic three-dimensional harmonic oscillator. The first model provides good agreement between calculations and experimental data, while the last two only reproduce their qualitative trend. Finally, the radial wavefunctions of the three aforementioned models, as well as their potential surfaces, are presented and discussed.

Approaching certain central-field Schrödinger problems through the harmonic oscillator green's function

We propose a method to analytically solve the Schrödinger equation with central potentials. By means of Sturm–Liouville expansion, the Green's function (GF) of the Schrödinger differential operator can be represented in the GF of the isotropic harmonic oscillator. Applying this straightforward approach, we solve the Schrödinger problem for the Coulomb potential, for the Morse potential, and for the respective inverted potentials. Using the obtained results, we formulate a method to analytically solve Schrodinger equations for multi-well potentials, demonstrating that classically forbidden regions can be treated as inverted harmonic oscillators, whereas harmonic oscillators or Morse potentials can approximate the wells.

Landau Levels versus Hydrogen Atom

The Landau problem and harmonic oscillator in the plane share a Hilbert space that carries the structure of Dirac’s remarkable so(2,3) representation. We show that the orthosymplectic algebra osp(1|4) is the spectrum generating algebra for the Landau problem and, hence, for the 2D isotropic harmonic oscillator. The 2D harmonic oscillator is in duality with the 2D quantum Coulomb–Kepler systems, with the osp(1|4) symmetry broken down to the conformal symmetry so(2,3). The even so(2,3) submodule (coined Rac) generated from the ground state of zero angular momentum is identified with the Hilbert space of a 2D hydrogen atom. An odd element of the superalgebra osp(1|4) creates a pseudo-vacuum with intrinsic angular momentum 1/2 from the vacuum. The odd so(2,3)-submodule (coined Di) built upon the pseudo-vacuum is the Hilbert space of a magnetized 2D hydrogen atom: a quantum system of a dyon and an electron. Thus, the Hilbert space of the Landau problem is a direct sum of two massless unitary so(2,3) representations, namely, the Di and Rac singletons introduced by Flato and Fronsdal.

Employing an operator form of the Rodrigues formula to calculate wavefunctions without differential equations

The factorization method of Schrödinger shows us how to determine the energy eigenstates without needing to determine the wavefunctions in position or momentum space. A strategy to convert the energy eigenstates to wavefunctions is well known for the one-dimensional simple harmonic oscillator by employing the Rodrigues formula for the Hermite polynomials in position or momentum space. In this work, we illustrate how to generalize this approach in a representation-independent fashion to find the wavefunctions of other problems in quantum mechanics that can be solved by the factorization method. We examine three problems in detail: (i) the one-dimensional simple harmonic oscillator; (ii) the three-dimensional isotropic harmonic oscillator; and (iii) the three-dimensional Coulomb problem. This approach can be used in either undergraduate or graduate classes in quantum mechanics.

Classical trajectories from the zeros of the quantum potential: the 2D isotropic harmonic oscillator

In the first part of this work, using the quantum potential approach, we show that a solution to the time-independent Schrödinger equation determines a subset of classical solutions, only if the region corresponding to the zeroes of the quantum potential is tangent to the caustic region determined by the classical trajectories. Thus, the solutions of the time-independent Schrödinger equation, according to their caustic and the zeros of the quantum potential, can be classified in three different cases given by the following conditions: the two regions coincide, they are tangent at certain subset of points, and the two regions are not tangent at any point. In the second part, as examples of the first type of wave functions, we present the solutions of the Schrödinger equation for the 2D isotropic harmonic oscillator, which are eigenfunctions of both the Hamiltonian and the angular momentum operators. That is, we show that for this family of solutions, the zeroes of the quantum potential coincide with the caustic region. Furthermore, we find that the classical trajectories, determined from the quantum ones and the zeroes of the quantum potential, conform to a family of elliptical curves for a particle with energy, (2n + l + 1)ℏ ω, and orbital angular momentum l ℏ.

Towards quantum mechanics on the curved cotangent bundle

The minimal-length paradigm is a cornerstone of quantum gravity phenomenology. Recently, it has been demonstrated that minimal-length quantum mechanics can alternatively be described as an undeformed theory set on a nontrivial momentum space. However, there is no fully consistent formulation of these theories beyond Cartesian coordinates in flat space and, in particular, no position representation. This paper is intended to take the first steps in bridging this gap. We find a natural position representation of the position and momentum operators on general curved cotangent bundles. In an expansion akin to Riemann normal coordinates with curvature in both position and momentum space, we apply the formalism perturbatively to the isotropic harmonic oscillator and the hydrogenic atom. Due to the symmetry of the harmonic oscillator under exchange of positions and momenta, we show that it is impossible to distinguish position- from momentum-space curvature with oscillators alone. Thus, we obtain an instantiation of Born reciprocity on the curved cotangent bundle, i.e. in precisely the way Born originally envisioned. It manifests itself as a symmetry mixing UV and IR physics, reminiscent of T-duality in string theory.

Quarkonium production in high energy pp collisions

We investigate the charmonium and bottomonium production in pp collisions using the Wigner densities formalism. The Wigner density of the quarkonia is approximated by analytical 3-D isotropic harmonic oscillator Wigner densities with the same root-mean-square radius given by the solution of the Schrödinger equation. This approach reproduces quite well the available experimental transverse momentum and rapidity distributions.

Warning! The negative divergence of the stress-tensor does not always yield the Ehrenfest force.

It has been assumed that the negative divergence of all stress tensors in common use yields the same force. This work finds that this is untrue, and, in fact, can vary wildly. We demonstrate this for the hydrogen atom, the one-particle isotropic harmonic oscillator, and a particle in an infinite spherical well where the exact density, pair-density, and the first order reduced density matrix are known for ground and excited states without any approximation. The Ehrenfest stress-tensor is introduced as any stress-tensor whose negative divergence will yield the corresponding Ehrenfest force for the same system when the exact wave-function is utilized. Stress-tensors within the literature are examined to show those that are Ehrenfest stress-tensors. Those that differ are demonstrated by how they differ within an exact formulation. The proof that the negative divergence of an Ehrenfest stress-tensor yields the Ehrenfest force is summarized.

Three-dimensional non-approximate Coulomb interaction between two trapped quantum particles

The two trapped quantum particles interacting problem is generalized to three dimensions, and the exact Coulomb potential is used. The system is solved by expanding the wavefunction in terms of the isotropic harmonic oscillator eigenfunctions and Hydrogen atom eigenfunctions independently, showing that each one results in a prime approximation for different domains of the normalized coupling constant of the relative interactions, suggesting that the combination of the basis is enough to build a well-suited base for the non-approximate problem. The results are compared to previous works that use a model of approximate finite-rage soft-core interaction model of the problem to give insights into the many-body states of strongly correlated ultracold bosons and fermions. We conclude that the proposed three-dimensional approach facilitates the distinction between bosons and fermions while the solutions given by the expansions better define the behavior of the particles for repulsive potentials. In addition, we discuss the substantial differences between our work and the previous approximate model.

Elliptic billiard with harmonic potential: Classical description.

The classical dynamics of the isotropic two-dimensional harmonic oscillator confined by an elliptic hard wall is discussed. The interplay between the harmonic potential with circular symmetry and the boundary with elliptical symmetry does not spoil the separability in elliptic coordinates; however, it generates nontrivial energy and momentum dependencies in the billiard. We analyze the equimomentum surfaces in the parameter space and classify the kinds of motion the particle can have in the billiard. The winding numbers and periods of the rotational and librational trajectories are analytically calculated and numerically verified. A remarkable finding is the possibility of having degenerate rotational trajectories with the same energy but different second constant of motion and different caustics and periods. The conditions to get these degenerate trajectories are analyzed. Similarly, we show that obtaining two different rotational trajectories with the same period and second constant of motion but different energy is possible.

Rational extensions of the Dunkl oscillator in the plane and exceptional orthogonal polynomials

It is shown that rational extensions of the isotropic Dunkl oscillator in the plane can be obtained by adding some terms either to the radial equation or to the angular one obtained in the polar coordinates approach. In the former case, the isotropic harmonic oscillator is replaced by an isotropic anharmonic one, whose wave functions are expressed in terms of [Formula: see text]-Laguerre exceptional orthogonal polynomials. In the latter, it becomes an anisotropic potential, whose explicit form has been found in the simplest case associated with [Formula: see text]-Jacobi exceptional orthogonal polynomials.

Exactly solvable model behind Bose-Hubbard dimers, Ince-Gauss beams, and aberrated optical cavities

The authors present a widely applicable, solvable model (an isotropic two-dimensional harmonic oscillator subject to two simple perturbations) and show that it underlies a broad range of physical phenomena displaying a topological transition, including aberrated optical cavities that support Ince-Gauss beams as their modes and the Bose-Hubbard dimer describing two coupled superfluids.

Scaling asymptotics of Wigner distributions of harmonic oscillator orbital coherent states

The main result of this article gives scaling asymptotics of the Wigner distributions of isotropic harmonic oscillator orbital coherent states concentrating along Hamiltonian orbits γ in shrinking tubes around γ in phase space. In particular, these Wigner distributions exhibit a hybrid semi-classical scaling. That is, simultaneously, we have an Airy scaling when the tube has radius normal to the energy surface Σ E , and a Gaussian scaling when the tube has radius tangent to Σ E .

Uncertainty Relation and the Thermal Properties of an Isotropic Harmonic Oscillator (IHO) with the Inverse Quadratic (IQ) Potentials and the Pseudo-Harmonic (PH) with the Inverse Quadratic (IQ) Potentials

The solutions for a combination of the isotropic harmonic oscillator plus the inversely quadratic potentials and a combination of the pseudo-harmonic with inversely quadratic potentials has not been reported, though the individual potentials have been given attention. This study focuses on the solutions of the combination of the potentials, as stated above using the parametric Nikiforov–Uvarov (PNV) as the traditional technique to obtain the energy equations and their corresponding unnormalized radial wave functions. To deduce the application of these potentials, the expectation values, the uncertainty in the position and momentum, and the thermodynamic properties, such as the mean energy, entropy, heat capacity, and the free mean energy, are also calculated via the partition function. The result shows that the spectra for the PHIQ are higher than the spectra for the IHOIQ. It is also shown that the product of the uncertainties obeyed the Heisenberg uncertainty relation/principle. Finally, the thermal properties of the two potentials exhibit similar behaviours.

Electronic structure and terahertz intersubband absorption spectra of phosphorus donor in silicon using a variationally optimized diagonalization method

The electronic structure and terahertz (THz) intersubband absorption spectra of phosphorus (P) donor in silicon (Si) are theoretically investigated using a variationally optimized diagonalization method based on the three-dimensional (3D) nonorthogonal associated Laguerre basis set, the basis set of eigenfunctions of 3D hydrogen atom and the basis set of eigenfunctions of 3D isotropic harmonic oscillator. A fully analytical expression of the matrix elements of the Hamiltonian is derived. Due to the large band anisotropy of Si, the donor electron series appear in an anomalous energy-level order and interval of the orbital angular momentum. The donor electron energy levels with even a small number of basis functions are in excellent agreement with those reported in the previous literature. Compared to two different basis sets of eigenfunctions of 3D hydrogen atom and isotropic harmonic oscillator, my proposed variationally optimized diagonalization method based on the 3D nonorthogonal associated Laguerre basis set is proven to be more efficient and reliable in calculating highly accurate impurity states of anisotropic 3D materials. The intersubband absorption spectra of P donor electron involving transitions from the 1s to excited states are strongly dependent on the orientations of the valley and the polarization vector of the incident light, which results from the fact that the large band anisotropy of Si induces strong mixing of the donor electron states. Their absorption frequencies are found to fall into the THz part of the electromagnetic spectrum. These results show that P donor in Si provides a material platform that enables intersubband devices operating in the THz frequency domain. • A fully analytical expression of the matrix elements of the donor Hamiltonian is derived. • The donor electron energy levels are in excellent agreement with those reported in the previous literature. • The intersubband absorption spectra of P donor electron are strongly dependent on the orientations of the incident light polarization. • P donor in Si provides a material platform that enables optoelectronic devices operating in the THz frequency domain. • The 3D nonorthogonal Laguerre basis is more efficient and reliable in calculating impurity states of anisotropic 3D materials.

Fisher information of a modified trigonometric inversely quadratic potential

• A new form of potential was proposed • The solution of the Schrödinger equation with the potential was obtained. • Using expectation values, the Fisher information was calculated. • The Cramer-Rao inequality was verified • The results of some special cases of the potential were obtained. The solution of the non-relativistic wave equation for modified trigonometric inversely quadratic potential (MTIAP) is obtained. The energy equation obtained, and the Hamiltonian which comprises the inverse squared term and the potential function were used to obtain the various radial expectation values. These expectation values were used to estimate the Fisher information for position configuration and momentum space. Numerical results were computed for Fisher information-based uncertainty relation with m = 0 and m = 1 . Some useful potentials like the isotropic harmonic oscillator and hydrogen atom were deduced from the modified trigonometric inversely quadratic potential and the Fisher information for the two potentials was also studied.

Islands of shape coexistence from single-particle spectra in covariant density functional theory

Using covariant density functional theory with the DDME2 functional and labeling single-particle energy orbitals by Nilsson quantum numbers, a search for particle-hole (p-h) excitations connected to the appearance of shape coexistence is performed for $Z=38$ to 84. Islands of shape coexistence are found near the magic numbers $Z=82$ and $Z=50$, restricted in regions around the relevant neutron midshells $N=104$ and $N=66$, respectively, in accordance with the well-accepted p-h interpretation of shape coexistence in these regions, which we call neutron-induced shape coexistence, since the neutrons act as elevators creating holes in the proton orbitals. Similar but smaller islands of shape coexistence are found near $N=90$ and $N=60$, restricted in regions around the relevant proton midshells $Z=66$ and $Z=39$, respectively, related to p-h excitations across the three-dimensional isotropic harmonic-oscillator magic numbers $N=112$ and $N=70$, which correspond to the beginning of the participation of the opposite parity orbitals $1{i}_{13/2}$ and $1{h}_{11/2}$, respectively, to the onset of deformation. We call this case proton-induced shape coexistence, since the protons act as elevators creating holes in the neutron orbitals, thus offering a possible microscopic mechanism for the appearance of shape coexistence in these regions. In the region around $N=40$, $Z=40$, an island is located on which both neutron p-h excitations and proton p-h excitations are present.

Scaling asymptotics of spectral Wigner functions* *SZ was partially supported by NSF Grant DMS-1810747. BH was funded by NSF CAREER Grant DMS-2143754 as well as NSF Grants DMS-1855684, DMS-2133806 and an ONR MURI on Foundations of Deep Learning.

We prove that smooth Wigner–Weyl spectral sums at an energy level E exhibit Airy scaling asymptotics across the classical energy surface Σ E . This was proved earlier by the authors for the isotropic harmonic oscillator and the proof is extended in this article to all quantum Hamiltonians −ℏ 2Δ + V where V is a confining potential with at most quadratic growth at infinity. The main tools are the Herman–Kluk initial value parametrix for the propagator and the Chester–Friedman–Ursell normal form for complex phases with a one-dimensional cubic degeneracy. This gives a rigorous account of Airy scaling asymptotics of spectral Wigner distributions of Berry, Ozorio de Almeida and other physicists.

More on Superintegrable Models on Spaces of Constant Curvature

A known general class of superintegrable systems on 2D spaces of constant curvature can be defined by potentials separating in (geodesic) polar coordinates. The radial parts of these potentials correspond either to an isotropic harmonic oscillator or a generalised Kepler potential. The angular components, on the contrary, are given implicitly by a transcendental, in general, equation. In the present note, devoted to the previously less studied models with the radial potential of the generalised Kepler type, a new two-parameter family of relevant angular potentials is constructed in terms of elementary functions. For an appropriate choice of parameters, the family reduces to an asymmetric spherical Higgs oscillator.

Interface asymptotics of Wigner—Weyl distributions for the Harmonic Oscillator

We prove several types of scaling results for Wigner distributions of spectral projections of the isotropic Harmonic oscillator on $\mathbb R^d$. In prior work, we studied Wigner distributions $W_{\hbar, E_N(\hbar)}(x, \xi)$ of individual eigenspace projections. In this continuation, we study Weyl sums of such Wigner distributions as the eigenvalue $E_N(\hbar)$ ranges over spectral intervals $[E - \delta(\hbar), E + \delta(\hbar)]$ of various widths $\delta(\hbar)$ and as $(x, \xi) \in T^*\mathbb R^d$ ranges over tubes of various widths around the classical energy surface $\Sigma_E \subset T^*\mathbb R^d$. The main results pertain to interface Airy scaling asymptotics around $\Sigma_E$, which divides phase space into an allowed and a forbidden region. The first result pertains to $\delta(\hbar) = \hbar$ widths and generalizes our earlier results on Wigner distributions of individual eigenspace projections. Our second result pertains to $\delta(\hbar) = \hbar^{2/3}$ spectral widths and Airy asymptotics of the Wigner distributions in $\hbar^{2/3}$-tubes around $\Sigma_E$. Our third result pertains to bulk spectral intervals of fixed width and the behavior of the Wigner distributions inside the energy surface, outside the energy surface and in a thin neighborhood of the energy surface.