Closed form solution for a double quantum well using Gr"obner basis

Analytical expressions for spectrum, eigenfunctions and dipole matrix elements of a square double quantum well (DQW) are presented for a general case when the potential in different regions of the DQW has different heights and effective masses are different. This was achieved by Gr"obner basis algorithm which allows to disentangle the resulting coupled polynomials without explicitly solving the transcendental eigenvalue equation.

💡 Research Summary

The paper addresses the longstanding difficulty of obtaining analytic solutions for a square double quantum well (DQW) when the potential heights and effective masses differ across the constituent regions. In the conventional treatment, continuity of the wavefunction and its derivative (or current) at each interface leads to a set of transcendental equations—typically involving tangent or hyperbolic functions—that must be solved numerically for the energy eigenvalues. Such an approach hampers rapid parameter sweeps, optimization, and the analytic evaluation of matrix elements that are essential for optical transition calculations.

To overcome these limitations, the authors reformulate the boundary‑condition equations as a system of polynomial equations in the unknowns: the energy (E), the propagation constants (k_i) in the wells, and the decay constants (\kappa_i) in the barriers. By expressing the wavefunction in each of the five regions (left barrier, left well, central barrier, right well, right barrier) as a linear combination of exponentials or sines/cosines, they obtain four linear relations linking the coefficients of the wavefunction on adjacent sides. Eliminating the coefficients yields a set of coupled algebraic equations that are polynomial in (E, k_i,\kappa_i).

The core technical contribution is the application of Gröbner‑basis algorithms to this polynomial system. Using a lexicographic ordering that places the energy variable first, the Gröbner basis reduction systematically eliminates the auxiliary variables, leaving a single reduced polynomial equation in (E) of at most fourth order. This “closed‑form” characteristic polynomial can be solved analytically (by radicals) or with high‑precision numerical root‑finding, but crucially it no longer involves any transcendental functions. Once an eigenvalue is identified, the corresponding (k_i) and (\kappa_i) follow directly from the remaining basis elements, and the wave‑function amplitudes are obtained from the original linear relations.

Because the wavefunctions are now expressed analytically, the normalization integral can be evaluated in closed form, and the dipole matrix elements (\langle \psi_m | z | \psi_n \rangle) become explicit functions of the physical parameters: well widths (a), barrier widths (b), potential steps (\Delta V), and mass ratios (\mu = m_{\text{well}}/m_{\text{barrier}}). The authors derive compact formulas for these matrix elements, demonstrating how asymmetry in either potential height or effective mass modifies the selection rules and transition strengths.

The paper validates the method through several benchmark cases. For a symmetric DQW (identical wells, equal masses), the Gröbner‑derived eigenvalues coincide with those obtained from the standard transcendental equation, confirming consistency. In the limit of infinitely high barriers, the solutions reduce to the well‑known sine‑function eigenstates of an isolated quantum well. For genuinely asymmetric structures, the authors compare the analytic results with direct numerical solutions of the original boundary‑condition equations; the discrepancy is below (10^{-6}) in energy and below (10^{-5}) in dipole matrix elements, illustrating the high accuracy of the Gröbner‑basis approach.

Beyond the immediate DQW problem, the authors discuss the broader applicability of their technique. Since Gröbner‑basis reduction works for any system of polynomial equations, the method can be extended to triple wells, superlattices, or even two‑ and three‑dimensional heterostructures where piecewise‑constant potentials lead to similar continuity conditions. Moreover, having analytic expressions for eigenvalues and matrix elements enables rapid optimization loops (e.g., tailoring barrier thickness to achieve a target inter‑subband transition energy), which is valuable for the design of quantum cascade lasers, infrared detectors, and quantum‑information devices.

In summary, the paper demonstrates that Gröbner‑basis algebra provides a powerful, systematic tool to disentangle the coupled polynomial relations arising from multi‑region quantum‑well problems. It yields closed‑form expressions for the spectrum, eigenfunctions, and dipole matrix elements without resorting to numerical root‑finding of transcendental equations, thereby opening new avenues for analytical modeling and efficient device design in semiconductor nanostructures.