---
title: Position-Dependent Mass Systems
url: https://www.emergentmind.com/topics/position-dependent-mass-systems
type: topic
---

# Position-Dependent Mass Systems

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) \) [2002.02297][1404.2650][1208.2300]. 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 [1411.7152][1110.1582][1412.4232]. 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)) \) [1411.7152][2002.02297][1512.09220].

## 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 [1110.1582][2007.11184]. 

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 \) [1110.1582].

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 [2007.11184]. 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 [1512.09220][1606.05780].

## 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 [1605.02336][1412.4232][2403.01235][2407.20112][1011.0708]. 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 [2403.01235].
- 27 systems invariant under either dilatation or shift symmetries with explicit lists of mass and potential forms admitting multiple second-order commuting integrals [2407.20112].
- 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 [1605.02336].

All these families admit explicit analytic solutions via separation of variables, Liouville-type changes of coordinates, or via their shape invariance/SUSY structures [1412.4232][1512.09220].

## 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 [1110.1582][1809.04469]. 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 [1809.04469].

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 [2503.16748]. 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 [1812.07756]. 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 [2503.23043][1606.05780]. 

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 [1512.09220]. 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 [1507.05217]. 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) \) [1110.1582] | Rotational: \( H = -\nabla \cdot [f(r)\nabla] + V(r) \), \( f(r)=1/[2m(r)] \) [1412.4232] |
| Representative mass function(s)     | \( m(x) = m_0/(1+\gamma x)^2 \), \( m(x) = m_0/(1+\kappa^2 x^2) \) [1110.1582][2007.11184] | \( m(r) = r^{2n} \), \( m(r) = 1/(r^2) \), \( m(r) = (r^4-1)^2/(r^4-2kr^2-1) \) [1412.4232]       |
| Exactly solvable potential(s)       | Infinite well, Pöschl–Teller, Harmonic Oscillator [1110.1582][1512.09220] | Coulomb, Harmonic oscillator, double shape-invariant potentials [1412.4232]        |
| Coherent states / algebra           | Nonlinear f-coherent, SU(2) [1606.05780][1002.4590][2503.23043] | su(1,1), su(2), shape-invariant SUSY [1512.09220][1412.4232]      |

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:**  
[1110.1582]  
[2002.02297]  
[2007.11184]  
[1812.07756]  
[1411.7152]  
[2503.16748]  
[1809.04469]  
[1804.01074]  
[1404.2650]  
[1208.2300]  
[1606.05780]  
[2403.01235]  
[2407.20112]  
[2503.23043]  
[1512.09220]  
[1002.4590]  
[1605.02336]  
[1412.4232]  
[1011.0708]  
[1507.05217]

Source: https://www.emergentmind.com/topics/position-dependent-mass-systems