Position-Dependent Mass Systems

Updated 6 February 2026

Position-dependent mass systems are models in which a particle’s effective mass explicitly varies with its position, leading to modified kinetic terms and dynamic behavior.
These systems require reexamination of classical and quantum formulations, employing novel canonical transformations and operator ordering to maintain Hermiticity and conserve energy.
Their study enables exactly solvable models with practical applications in semiconductor heterostructures, curved quantum systems, and supersymmetric frameworks.

A position-dependent mass (PDM) system describes a classical or quantum mechanical model in which the effective mass of a particle depends explicitly on its position. PDM frameworks are motivated by diverse contexts, including charge carriers in semiconductor heterostructures, inhomogeneous media, curved spaces, Kaluza–Klein dimensional reduction, and nonclassical models in mathematical physics. They necessitate a re-examination of foundational structures: the kinetic term, the Hamiltonian and Lagrangian formulations, symmetries and superintegrability, as well as the quantization procedure and associated algebras.

1. Classical and Quantum Frameworks for PDM

In classical mechanics, the Lagrangian is generalized to $L(x, \dot{x}) = \frac{1}{2} m(x) \dot{x}^2 - V(x)$ with position-dependent mass $m(x)$ (Rosas-Ortiz, 2020, Lopez et al., 2014, Cruz et al., 2012). The corresponding equations of motion acquire extra terms beyond the constant-mass Newtonian form: $m(x) \ddot{x} + m'(x) \dot{x}^2 = F(x)$ Additional "thrust" terms quadratic in velocity appear. The canonical momentum, $p = m(x)\dot{x}$ , and the naive Hamiltonian, $H = p^2/(2m(x)) + V(x)$ , generally do not yield a conserved energy, as $dH/dt \neq 0$ . Constants of motion must be constructed to restore Hamiltonian structure, often via canonical transformations or by introducing effective invariants.

In quantum theory, the quantization of the kinetic term is nontrivial due to the operator-valued mass, leading to operator-ordering ambiguities. The general Hermitian kinetic operator can be written in the von Roos form: $\hat{T}_{\alpha,\beta,\gamma} = \frac{1}{4}\left[m^{\alpha}(\vec{x})\, \hat{p}\, m^{\beta}(\vec{x})\, \hat{p}\, m^{\gamma}(\vec{x}) + \text{H.c.}\right], \quad \alpha + \beta + \gamma = -1$ Canonical choices (e.g., BenDaniel–Duke, Gora–Williams, Zhu–Kroemer) reflect different physical constraints and symmetries (Ruby et al., 2014, Filho et al., 2011, Nikitin, 2014). Exact solvability and the removal of ordering ambiguities are often achieved by employing point-canonical transformations (PCTs), with canonical coordinates $s(x) = \int^x \sqrt{m(x')} dx'$ and wavefunction scalings $\psi(x) = m^{1/4}(x)\Phi(s(x))$ (Ruby et al., 2014, Rosas-Ortiz, 2020, Bravo et al., 2015).

2. Algebraic Structures and Operator Ordering

The PDM setting induces deformed algebraic structures. A deformed translation operator, $\mathcal{T}_\gamma$ , generates a generalized momentum $\hat{p}_\gamma = -i\hbar(1 + \gamma x)\frac{d}{dx}$ , which obeys a modified commutation relation $[\hat{x}, \hat{p}_\gamma] = i\hbar(1 + \gamma x)$ and leads to position-dependent uncertainty relations (Filho et al., 2011, Costa et al., 2020).

For a kinetic energy operator in the BenDaniel–Duke ordering, the Schrödinger equation takes the divergence form: $-\frac{\hbar^2}{2} \frac{d}{dx} \left[ \frac{1}{m(x)} \frac{d}{dx} \right] \psi(x) + V(x)\psi(x) = E\psi(x)$ Hermiticity and physical compatibility in boundaries/interfaces are ensured by matching the continuity of wavefunctions and the generalized flux $[1/m(x)]d\psi/dx$ (Filho et al., 2011).

When the mass function is tied to deformed algebras, such as κ-algebras (motivated by Kaniadakis statistics), classical and quantum operators are mapped via PCTs to constant-mass systems in deformed spaces: $m(x) = m_0/(1 + \kappa^2 x^2)$ , with generalized momentum and position operators reflecting the algebra (Costa et al., 2020). Supersymmetry (SUSY) and shape-invariance properties are extended to PDM, generating exactly solvable models and spectrum-generating algebras (e.g., su(1,1), su(2)) for reflectionless and finite-gap settings (Bravo et al., 2015, Amir et al., 2016).

3. Superintegrability, Symmetries, and Exactly Solvable Families

Extensive algebraic classification has established large families of integrable, superintegrable, and maximally superintegrable quantum systems with PDM in one, two, and three dimensions, invariant under various Lie group symmetries (Rañada, 2016, Nikitin, 2014, Nikitin, 2024, Nikitin, 2024, Ballesteros et al., 2010). In 2D and 3D, these are constructed by enforcing the existence of quadratic or higher-order integrals of motion (e.g., Killing tensors compatible with the kinetic energy metric) and exploiting the underlying symmetry algebra (e.g., so(3), so(4), conformal algebras).

Key structural results include:

Complete classification of 68 inequivalent 3D systems with cylindrical symmetry, among which 27 are superintegrable and 12 maximally superintegrable (Nikitin, 2024).
27 systems invariant under either dilatation or shift symmetries with explicit lists of mass and potential forms admitting multiple second-order commuting integrals (Nikitin, 2024).
General families of 2D superintegrable models with radial PDM of the form $m(r) = r^{2n}$ and potential terms deforming the oscillator and Kepler problems, with quadratically conserved integrals arising from deformed Fradkin or Laplace–Runge–Lenz constructs (Rañada, 2016).

All these families admit explicit analytic solutions via separation of variables, Liouville-type changes of coordinates, or via their shape invariance/SUSY structures (Nikitin, 2014, Bravo et al., 2015).

4. Applications: Quantum Wells, Thermodynamics, and Relativistic Systems

PDM frameworks are central in modeling graded semiconductor heterostructures, quantum wells, and quantum dots, where the effective electron mass varies with position (Filho et al., 2011, Serafim et al., 2018). In curved quantum systems, the position dependence of the mass fundamentally alters the geometric potential, localization, and transport properties. For instance, in deformed carbon nanotubes, the radial mass variation substantially renormalizes conductance resonances and band-gap structures; the kinetic operator includes additional derivative terms and mass-weighted geometric curvatures (Serafim et al., 2018).

In the context of equilibrium statistical mechanics, PDM alters all canonical thermodynamic potentials via $Z_1 = \int \sqrt{m(x)} e^{-\beta_0 V(x)} dx$ , and its inhomogeneity is naturally interpreted as a specific spatial distribution of local inverse temperature under the lens of superstatistics (Gomez et al., 20 Mar 2025). Quadratic and exponential mass profiles generate either typical 3D oscillator thermodynamics or anomalous specific heat behavior in agreement with the third law.

Relativistic extensions include Dirac systems with radially singular mass functions—e.g., $m(\rho) = m_0 + \kappa/\rho$ in the Aharonov–Bohm–Coulomb (ABC) background—leading to discrete spectra and bound states in regimes otherwise continuous for constant mass (Oliveira et al., 2018). Here, the PDM acts as a Lorentz-scalar coupling and modifies the degeneracy and density of states.

5. Coherent States, Supersymmetry, and Special Models

Exactly solvable PDM models support the construction of coherent states with generalized or nonlinear properties, sub-Poissonian photon statistics, and nonclassical phase-space signatures. For the harmonic PDM oscillator with $m(x) = m_0/(1 + \alpha x^2)^2$ , Gazeau–Klauder coherent states have been built explicitly, satisfying all completeness, overcompleteness, and temporal stability criteria; the corresponding Wigner functions display negative regions indicative of strong quantum character (Takou et al., 29 Mar 2025, Amir et al., 2016).

Finite-gap PDM systems, including Lamé and Darboux–Treibich–Verdier types, can be generated via coordinate reductions from integrable geometric backgrounds, with SUSY and nonlinear algebraic symmetries (Bravo et al., 2015). In particular, position-dependent mass models provide a unifying framework for analogues of the Higgs oscillator, Mathews–Lakshmanan oscillator, and reflectionless systems relevant in cosmological scenarios.

6. Physical Origins and Geometric/Field-Theoretic Realizations

PDM arises in effective-mass treatments in condensed matter, but field-theoretic contexts—e.g., Kaluza–Klein dimensional reduction of Einstein gravity with inhomogeneous compactification—yield PDM via dilaton couplings to 4D particles (Morris, 2015). The reduced effective mass is $m(x) = m_0 \exp(-a\phi(x)/2)$ with the dilaton field $\phi(x)$ controlling both geometry and mass variation. In such scenarios, both the classical and quantum equations show explicit modifications in energy levels and eigenfunction structure; e.g., a nonlinear oscillator with dilaton-induced PDM possesses a spectrum isospectral to the canonical oscillator, but with distinctly deformed wavefunctions and density profiles.

7. Tables of Representative Exact PDM Systems in Quantum Mechanics

Property	1D/2D Example	3D Example (Symmetry Type)
BenDaniel–Duke ordering	$H = -\frac{\hbar^2}{2} d/dx (m^{-1} d/dx) + V(x)$ (Filho et al., 2011)	Rotational: $H = -\nabla \cdot [f(r)\nabla] + V(r)$ , $f(r)=1/[2m(r)]$ (Nikitin, 2014)
Representative mass function(s)	$m(x) = m_0/(1+\gamma x)^2$ , $m(x) = m_0/(1+\kappa^2 x^2)$ (Filho et al., 2011, Costa et al., 2020)	$m(r) = r^{2n}$ , $m(r) = 1/(r^2)$ , $m(r) = (r^4-1)^2/(r^4-2kr^2-1)$ (Nikitin, 2014)
Exactly solvable potential(s)	Infinite well, Pöschl–Teller, Harmonic Oscillator (Filho et al., 2011, Bravo et al., 2015)	Coulomb, Harmonic oscillator, double shape-invariant potentials (Nikitin, 2014)
Coherent states / algebra	Nonlinear f-coherent, SU(2) (Amir et al., 2016, Dutra et al., 2010, Takou et al., 29 Mar 2025)	su(1,1), su(2), shape-invariant SUSY (Bravo et al., 2015, Nikitin, 2014)

In each class, operator ordering is fixed by physical or algebraic requirements (e.g., Hermiticity, boundary matching). Each system supports a full range of model-specific properties: ladder operators, dynamical algebras, superintegrable or maximally superintegrable structure, and, in special cases, dual shape-invariant factorizations that connect families of exactly solvable spectral problems.

References:

(Filho et al., 2011, Rosas-Ortiz, 2020, Costa et al., 2020, Oliveira et al., 2018, Ruby et al., 2014, Gomez et al., 20 Mar 2025, Serafim et al., 2018, Ghose-Choudhury et al., 2018, Lopez et al., 2014, Cruz et al., 2012, Amir et al., 2016, Nikitin, 2024, Nikitin, 2024, Takou et al., 29 Mar 2025, Bravo et al., 2015, Dutra et al., 2010, Rañada, 2016, Nikitin, 2014, Ballesteros et al., 2010, Morris, 2015)