Central delta barrier in a symmetric infinite well Created 2026-10-14 Updated 2026-10-18
Inside walls at , a central delta potential of strength leaves odd parity operator eigenstates unchanged because their wavefunctions vanish at zero. For even states, continuity and the derivative jump give the displayed quantization equation. For , one root lies between and for every .
For two identical fermions of mass in the same one-dimensional quantum harmonic oscillator, use center of mass and separation . The corresponding oscillator masses are and reduced mass , with the same frequency . If are their occupation numbers, the unperturbed energy is . The fermionic exchange symmetry requires even with a spin-one-half singlet state and odd with a spin-one-half triplet state.
For a delta potential the first-order shift is . Every odd relative wavefunction vanishes at contact, so these states remain exact unaffected eigenstates of the spatial problem. For the even relative ground state, the oscillator length is , where , and the shift is , independent of . This separates interaction-sensitive singlets from contact-insensitive triplets.
Finite rectangular barrier matching Created 2026-10-14 Updated 2026-10-20
Take a nonrelativistic particle of constant mass with energy , and choose the same zero of potential energy on both sides of a rectangular potential barrier:
The equality of the exterior potential energies matters: it gives the same exterior wavenumber and makes . Here is the reduced Planck constant. This model has a finite jump in potential energy, no delta potential, and no absorbing term.
For , the time-independent Schrödinger equation gives the spatial wavefunction
The full stationary state is . We set the incident probability amplitude to one and impose no incoming wave from the right; this scattering wavefunction is not a normalizable wavefunction. The tunnelling decay constant describes exponential decay, rather than an oscillatory wavenumber in the barrier. Both exponential functions are needed on the finite interval: the growing term is not excluded by a condition at infinity.
The wavefunction and its derivative are continuous at each face, giving four linear matching conditions:
Their solution supplies the transmission amplitude
A short derivation is to use the two right-face equations to express and in terms of , then substitute into the left-face equations. Taking the squared modulus gives the rectangular barrier transmission probability.
A hard wall at the origin and an attractive delta potential at distance form a leaky quantum cavity. For unit incident amplitude from outside, the interior standing wave has coefficient , where and . Its modulus squared is the title expression. A resonance pole produces an enhancement with a locally Lorentzian line shape, proportional to . At fixed low and small , and .
The coefficient can exceed one because it measures a standing-wave amplitude, while the net interior probability current is zero. Conservation of probability current gives unit reflected flux: no current passes through the hard wall. Near zero, when ; the threshold case requires a separate limit.
Write the trap angular frequency as and the single-particle oscillator length as . The center of mass coordinate and relative coordinate give independent quantum harmonic oscillators, both of frequency , with masses and reduced mass . Label their excitation numbers by and . The unperturbed energy eigenvalue is .
Particle exchange leaves unchanged and sends to . The exchange symmetry of two identical spin-one-half fermions therefore requires even to accompany the spin-one-half singlet state, and odd to accompany a spin-one-half triplet state. For the contact shifts of two spin-one-half fermions in a harmonic trap, the delta potential depends only on . For odd , the relative wavefunction vanishes at , so its contact shift is zero. Within the first-excited manifold, both perturbations are diagonal in the labels : this is a perturbation diagonal within a degenerate eigenspace. For , the relative oscillator length is , giving the first-order nondegenerate perturbation theory shift
The Zeeman splitting is for total spin projection . Thus the five lowest states, counted with multiplicity, have
under the assumed conditions . The three triplet states have slopes , while both singlet states are field independent at this order. The three triplet state branches coincide at zero magnetic field. Their crossing with the excited singlet state is allowed because the total-spin sectors do not mix. The odd relative wavefunctions are in fact unaffected by this delta potential at any coupling; the displayed singlet state shifts retain only first order.
Figure 1.
Five lowest two-fermion states in a weakly interacting harmonic trap
. Original plot for repulsive contact shift , with . The two horizontal singlet state branches and three triplet state branches count five states; coincident energies count as fewer distinct energy levels. An attractive coupling reverses the sign of and moves both singlet state branches downward. Details.
For a real potential energy and a stationary state, the probability current is
The time-independent Schrödinger equation gives , so . In an exterior plane wave superposition, interference terms contribute no current and
Consequently for every choice of incoming amplitudes. If , this says for every complex vector . Since is Hermitian, vanishing of its quadratic form implies that the matrix itself vanishes. Therefore
The second identity follows by choosing . Since is a square matrix, the first also implies , proving unitarity. A real delta potential preserves this proof: continuity of and the real derivative jump make the current continuous across the interaction point.
Assume and put , , and . At a delta potential located at , integration of the Schrödinger equation gives the matching conditions
Writing the adjacent plane wave coefficients as and gives
Solving these two linear equations yields the transfer matrix of a delta potential:
For the symmetric double-delta potential, use separate plane wave pairs in , , and . Successive matching gives between the two exterior pairs. For real , multiplication gives
The determinant identity reads . Left incidence means , whence
Here the first relation also follows from . Therefore the requested transmission probability and quantum reflection probability, as functions of the printed momentum , are
They obey . As checks, gives , and the coincident limit gives , the answer for a single delta potential of strength .
Take fixed and define . For sufficiently large , the barrier region has and solves . Solving this constant-coefficient equation gives the exact transfer relation
In the specified limit, and , while
The transfer matrix therefore tends to . As the two boundaries collapse to the origin, this proves continuity and the derivative jump for a delta potential:
The vanishing exterior phase shift introduces no further factor.
Write the incoming and reflected waves as on the left and the transmitted wave as on the right. Matching gives
The denominator is nonzero for real and , so the limiting matching system is well defined. The probability currents on either side have the same factor , making the transmission probability . Thus scattering by a delta potential gives
The limit retains a finite reflection probability despite the barrier's vanishing width because its integrated strength remains finite.
For a normalized trial state in the quadratic-form domain of a self-adjoint Hamiltonian operator bounded below, its expectation of is at least the ground energy . Minimizing this Rayleigh quotient over a tractable trial family therefore supplies an upper bound.
For the two states, set . Their quotient is
The strict assumption excludes , since equality in the overlap Cauchy-Schwarz bound would then force the two states to be proportional and . Thus . The two generalized eigenvalues are , and their difference has the sign of . The symmetric trial is the lower one, so
For the single attractive delta well, the supplied wavefunction has norm one because . Its kinetic energy in quadratic-form notation is , and its delta potential energy is . With this gives
Its derivative jump is , exactly the delta boundary condition, so it is an eigenstate. A decaying bound state on the two half-lines must have a common decay rate fixed by this same jump; there is only this one negative-energy state, confirming it is the ground state.
For , integrate the overlap on the three regions , , . Put . The result is
The second equality is the single-well energy plus the extra delta expectation; the third follows by applying either single-well eigenvalue equation in the cross matrix element. Moreover , so the preceding bound applies. Substitution gives
For widely separated wells, each localized orbital is close to an isolated-well ground state and the low-energy subspace is accurately their span; the symmetric combination captures the leading tunneling energy lowering. In the limit both exact and trial energies approach . This illustrates electronic molecular binding: sharing an orbital between attractive centers can lower its energy below the separated-center value. A full equilibrium bond length would also require the repulsion of the positive centers.
Insert the position and momentum resolutions of the identity to obtain the unitary Fourier transform
The kinetic term is multiplication by , and the potential has momentum kernel . Thus the Schrödinger equation becomes
For the attractive delta potential, , the kernel is the constant , giving the stated positive right-hand side. For with , constancy of the left side requires . Since , its coefficient equation gives
Finally, transforming gives , verifying the same momentum-space form and energy.
The parity operator acts by . For an even potential, differentiating twice and using proves on the parity-invariant domain of the Hamiltonian operator. Therefore each energy eigenspace is invariant under . The projectors and split it into even and odd subspaces, giving a basis of simultaneous eigenstates. An odd state obeys , hence .
The infinite walls impose , with the wavefunction zero outside. It is continuous at the central delta potential. Integrate the stationary Schrödinger equation over ; the energy integral tends to zero and the Dirac delta function contributes . This gives
For the given piecewise positive-energy form, . Replacing by swaps the two pieces and leaves their values unchanged, proving . The value at zero is , and the one-sided derivatives are and . The jump condition gives . The wall at gives , hence
For the undivided wall equation instead gives the usual free even roots. For , set and . In every interval there is exactly one root: writing gives , whose two sides are respectively strictly decreasing and increasing. These are all positive even roots; the quadratic form is nonnegative, so there are no negative eigenvalues. Thus and .
The printed limit of this absolute energy difference is false. The exact quantization equation gives , so but . Consequently
Thus the high-energy shift from a central delta barrier is a finite positive constant, not zero. The corrected vanishing statement is the relative shift , or the wave-number difference . Physically, highly energetic particles are weakly affected in relative terms by a fixed point barrier. The free even wavefunction normalized in has density at the origin, so the barrier's first-order expectation is , matching the exact limit.
Odd states vanish at zero, so the delta term and the derivative jump both vanish. Their normalized wavefunctions and energies are
with zero wavefunction outside. They do not depend on . These even and odd sectors describe the central delta barrier in a symmetric infinite well.
With outgoing-channel index first, the scattering matrix convention is as . An even potential commutes with parity, so its incoming even and odd channels remain separate: .
For the odd channel, on write , while on its prescribed normalization is . A delta potential keeps the wavefunction continuous and imposes , with . Thus the exterior logarithmic derivative is . Matching gives
Multiplying numerator and denominator by and using exponential forms yields
For real , the numerator of the fraction is the complex conjugate of its denominator, and the prefactor has modulus one. Hence , conservation of probability flux in this single elastic parity channel. In particular the minus sign built into makes the free odd-channel matrix element equal to one.
For an odd bound state, continue to , . A pole means the outgoing solution is exponentially decaying on both sides. Its denominator vanishes when
The function rises strictly from to infinity: differentiation reduces positivity to . Thus there is exactly one positive solution if , and none otherwise. The numerator does not vanish at that solution, so the pole is genuine. Consequently
At equality the zero-energy threshold solution is not square integrable.
For repulsive strong barriers put and . The pole equation is . Since , multiplication by gives . Iterating in yields
This lower-half-plane resonance pole is an odd standing wave trapped between the barriers and slowly leaking by tunnelling. Its complex energy is , with . The factor then decays rather than grows; encodes the positive escape rate.
Put and . Integrating the Schrödinger equation across the delta potential gives continuity of and
For incidence from the left, write on and on . Matching gives and , hence
The incident and transmitted probability currents have the same velocity factor, so the reflection coefficient and transmission coefficient are
For , an even bound state has , . The derivative jump now gives , so and
There is exactly one positive value of , hence exactly one even bound state. Its normalized amplitude is . An odd decaying solution would have to vanish at zero and is therefore trivial.
In the shrinking finite square well, its integrated potential is , not : the full width is . With fixed and , . The matching condition then bounds by , so remains bounded. Consequently , and the previous matching condition gives
Using recovers the same bound-state energy. All other even levels disappear as the width shrinks.
A bound state has and , since outside the delta potentials a square-integrable solution must decay exponentially. The wavefunction is continuous at each delta: a jump would produce a delta derivative in with no matching term in the equation. Integrating the time-independent Schrödinger equation across either gives
The symmetry of the symmetric double-delta potential permits even and odd solutions. Up to a nonzero normalization factor, the even candidate is for , and outside. Its derivative jump at gives
The odd candidate is inside and outside. Its jump gives
Parity ensures the jump at as well. Put and . The continuous strictly increasing function ranges from to infinity, so there is exactly one even bound state for every . The continuous strictly increasing function ranges from to infinity, so the odd bound state of a symmetric double-delta potential exists if and only if . Strict increase follows from the allowed monotonicity facts plus the strictly increasing term . Therefore
At equality the odd limiting solution has and is not square integrable, so it is not a bound state. No nonnegative energy can supply another full-line square-integrable state, since its exterior solutions are oscillatory or affine.
In the Parity basis of a one-dimensional S-matrix, parity acts as . Parity invariance requires
and therefore
Independently, unitarity gives
Combining parity and unitarity leaves two phases,
For the symmetric double-delta potential, consider the odd channel on . Choose the asymptotic convention
and write for ; the sine enforces odd parity at the origin. Continuity and the delta potential derivative jump at give
Eliminating yields the Odd-parity S-matrix of a symmetric double-delta potential
Using rewrites this as the requested expression
For real , the two factors in the ratio are complex conjugates, directly confirming .
A bound state occurs at a pole , . Setting the denominator to zero gives
or
For , write . The right-hand side is increasing and concave, starts at zero with slope , and approaches , whereas the left-hand side is the straight line . A positive intersection therefore exists exactly when
This is the odd bound state of a symmetric double-delta potential; equality is its zero-energy threshold.
For , analytically continue the same pole equation into the lower half of the complex -plane. Put
It becomes
When , the lowest odd mode lies near the hard-wall value . Writing and expanding gives
so
In particular, the odd resonance of a strongly repulsive symmetric double-delta potential has
The negative imaginary part identifies a decaying scattering resonance. If
then, to leading nonzero order,
Its time-dependent amplitude is proportional to
so its survival probability decays as and its lifetime is
Physically, this is the lowest odd standing wave trapped for a long time between two strong repulsive barriers and leaking out by quantum tunnelling.
The required dipole matrix element is
where differentiating the elementary Laplace integral for evaluates the last integral. Hence
Put and . The continuum transition probability at long times uses
Since
the total escape probability becomes
For the attractive delta potential,
Substitution gives the requested leading-order result:
Write . An even bound state is proportional to inside the well and to outside, where
Continuity of the logarithmic derivative at gives
For sufficiently small , lies in the first monotone branch, so this equation has exactly one positive root and the even state is unique up to scale.
As , . Hence
The wells converge distributionally to , whose normalized even bound state is
Its derivative jump is precisely the one imposed by the delta potential.
The parity operator acts as , with . For an even potential energy , so energy eigenstates can be chosen even or odd, with parity eigenvalues and . One-dimensional bound states are nondegenerate; the ground state is nodeless and even.
For the finite square well, square integrability excludes the growing exterior exponential, while even parity excludes the interior sine. Thus the interior solution is proportional to , with matching decaying tails and
Continuity at gives and . Their ratio yields
The conditions at follow from the same parity.
For a narrow weak well, , hence and
This is self-consistent because . It is the shallow finite-square-well energy, equivalently the attractive delta potential limit with integrated strength . An arbitrarily weak attractive well of this form still has its even bound ground state.
The first excited state must be odd. Replace the interior cosine by and choose exterior tails with opposite signs. The matching condition becomes . Its first root requires . At the binding threshold , so
is necessary and sufficient for the first odd bound state. Equality gives a zero-energy nonnormalizable threshold solution, not an additional bound state.
Equation and boundary conditions. The delta potential has strength with units of energy times length, because the Dirac delta function has units of inverse length. The time-independent Schrödinger equation is
The wavefunction is continuous at the origin. A finite jump in would produce a derivative of a Dirac delta function in , with no corresponding term in this potential. A jump in does produce a delta contribution and is allowed. Integrate through a symmetric interval:
For continuous finite , the right-hand side vanishes as . Consequently
This derivative jump is the boundary condition associated with the attractive delta potential.
Unique bound state. For , write . Away from the origin, . Normalizability excludes exponentials growing at either infinity, leaving
Continuity gives . The derivative jump then becomes . For a nonzero state,
There is exactly one positive value of , so exactly one negative-energy bound state. An odd solution would have ; the allowed decaying exponentials would then both vanish identically. Normalize the even solution using
Choosing the overall phase real and positive gives the bound state in a delta potential:
It peaks at the origin, has a downward cusp there, and decays exponentially on both sides. Its value, but not its slope, is continuous.
Scattering amplitudes and fluxes. For , put and retain as a strength parameter. A unit-amplitude wave incident from the left has
Continuity gives . The jump condition gives
Solving,
The probability current of a plane wave of amplitude is in its propagation direction. Since the asymptotic wavenumber is the same on both sides, reflection probability and transmission probability are just and :
They satisfy , expressing conserved probability current in this real potential.
Figure 1.
The normalized delta-well bound state and reflection/transmission probabilities as attraction increases
. Details.
The right-hand plot uses , proportional to at fixed incident energy. Then rises monotonically from zero to one, while falls from one to zero; they are equal at . Despite the absence of a classical forbidden barrier, an abrupt attractive interaction reflects a quantum wave. Increasing attraction strengthens the derivative mismatch and can approach complete reflection at fixed positive energy. At fixed , increasing incident energy instead makes the well increasingly transparent. The negative-energy bound state and the positive-energy scattering states are different stationary solutions, not an absorptive loss channel.
For real , taking the modulus squared of the amplitude gives the half-line delta-well amplitude enhancement
The supplied resonance poles have real parts and negative imaginary parts for fixed low at .[ref] Linearizing the denominator at a simple pole gives
Thus the sketch has narrow, tall peaks with a locally Lorentzian line shape, just above positive integer multiples of . Their half-widths are , and the leading peak heights are . In this fixed-low-mode approximation, increasing broadens and lowers the peaks. Exactly at , the value is one. At zero it has the finite limit for .
The peaks represent long-lived standing waves between the wall and the attractive delta potential. Despite the source's transmission-coefficient terminology, is an internal amplitude factor and can exceed one. The interior probability current is zero, and the reflected probability current equals the incident magnitude; the hard wall passes no transmitted flux.