The Pusey-Barrett-Rudolph theorem excludes overlapping preparation distributions for distinct pure states in an ontological model of a quantum system satisfying preparation independence and the specified quantum predictions. It excludes neither all hidden variables nor every informational use of a quantum state. The precise overlap criterion is essential: a psi-epistemic model lets at least two distinct pure-state preparations assign positive common probability mass to the same underlying states; a psi-ontic model has mutually singular measures for each distinct pair.[ref]
Let preparation sample from . A fixed quantum measurement has nonnegative response functions with . Averaging over the preparation must reproduce the Born rule. The response is fixed by and the complete underlying state, not by additional knowledge of the preparation label. No outcome determinism or Bell local causality is assumed; the joint measurement response need not factor into local responses. For independently operated preparations, preparation independence further identifies the joint state space with the appropriate product and assigns
Operationally independent devices do not logically entail this ontic factorization; it is an additional physical premise.[ref]
The two-qubit example can be checked directly. Set . In computational basis order , take the following normalized columns as an orthonormal basis:
Their pairwise inner products vanish. Direct multiplication shows that the respective preparations have zero probability amplitude for respective outcomes . This is the PBR exclusion measurement for zero and plus: the measurement does not perfectly distinguish the preparations, but each outcome rules one out.
Suppose a common submeasure of and has mass ; with densities one can take . Preparation independence makes a common submeasure of all four product preparations, of mass . The zero probability for outcome under its excluded preparation and nonnegativity force
Summing and using normalization yields , a contradiction. This measure formulation avoids treating a set with merely nonempty topological support as a positive-probability overlap.
For arbitrary distinct nonorthogonal pure states, write , with overlap and . The general proof constructs an exclusion measurement on independently prepared copies when . A finite always exists. The same common-submeasure contradiction then excludes positive overlap. Orthogonal preparations are already disjoint if their perfectly distinguishing measurement is reproduced. The original paper also bounds residual overlap for imperfect exclusion data; exact disjointness concerns the ideal theorem, not a finite experiment's proof of exact zero.[ref]
The scope is consequently sharp. Bohmian mechanics is not excluded: its underlying state includes the wave function, making distinct pure preparations disjoint even though extra configuration variables remain. A psi-ontic model need not be complete with the wave function alone, and the theorem does not identify the wave function with a three-dimensional material field. Models abandoning preparation independence, or accounts not representing preparations as uncertainty over observer-independent states of this specified kind, are outside this conclusion. Epistemic uncertainty between mixed preparations is not banned. Nor does the theorem establish a collapse mechanism, choose between Bohmian and Everettian theories, or solve the measurement problem. Its achievement is a conditional impossibility theorem about a precise form of pure-state epistemic overlap.
Quantum decoherence explains the practical stability of records and suppression of observable interference; it does not alone select one outcome. The measurement problem arises when universal unitary time evolution, a complete wave function description and a unique definite outcome are all demanded. A theory must explain which commitment is changed or supplemented, not merely replace a pure composite description by a reduced one.
Model an ideal quantum measurement and its environment by
Taking the partial trace over the environment gives the system–apparatus reduced density matrix
Environmental records make very small for macroscopically distinct alternatives. If independent fragments have overlaps , their combined overlap is , explaining the rapid suppression. The apparatus interaction also favours robust records, giving environment-induced pointer-basis selection rather than an arbitrary choice of basis.[ref]
What follows is approximate agreement, for local apparatus observables, with a classical probabilistic mixture. What does not follow is that either term has ceased to exist. Even with orthogonal environmental records, the global quantum state is a pure state with quantum entanglement. For orthogonal branches , the global observable has expectation in the superposition and zero in the corresponding incoherent mixture. Its practical inaccessibility is not its mathematical disappearance. The reduced state is an improper mixed state, not proof that an unknown single branch was selected.
Different interpretations use this result differently. In Bohmian mechanics, disjoint environmental configurations make one occupied branch an effective description of the actual apparatus. In the many-worlds interpretation, robust branches supply approximately autonomous relative records, without a unique global outcome. In the GRW model, a separate stochastic change of the wave function selects a dominant component; environmental decoherence helps explain ordinary records but is not the collapse law. These are different completions of the argument.
Nor does diagonality derive the Born rule. Interpreting the diagonal entries as experienced outcome probabilities requires a probability rule or an interpretive argument. Extremely small overlaps also leave approximations and possible quantum recoherence; a finite controlled environment can in principle be manipulated to reveal the global coherence. The appropriate verdict is that decoherence resolves much of the preferred-record-basis and effective-classicality problem while leaving the unique-outcome and probability problems dependent on further theory. It is an essential physical mechanism in many accounts, not a self-contained solution by relabelling a partial trace as collapse.
An empty wave is a nonzero component of the wave function in Bohmian mechanics containing no actual particle configuration. Its persistence is a question about ontology and explanation, not an inconsistent prediction of duplicate detected particles. In de Broglie-Bohm theory, the wave function evolves by the Schrödinger equation while the actual configuration evolves by the guidance equation. For spinless particles without a magnetic vector potential,
The physical state includes both elements, so a component of need not contain another set of particles.[ref]
After a position-recording interaction, write the universal wave function as , with apparatus configurations supported in disjoint regions. If the actual lies in the first region, and its relevant local derivatives vanish in the ideal disjoint-support model. Evaluating the guidance equation at therefore gives the velocity field of the occupied component alone. For realistic nearly disjoint packets this is an approximation justified by decoherence. The other component is empty in configuration space, even when the system packets overlap in ordinary SPACE. This is the effective wave function in Bohmian mechanics.
The empty component cannot simply be deleted from the fundamental theory. In a coherent interferometer, separated packets can later overlap and both contribute to the actual probability current. Their relative quantum phase can change the outgoing particle's path. An empty arm is therefore capable of future dynamical significance. When an environment records the alternatives in very many degrees of freedom, restoring their configuration-space overlap is overwhelmingly difficult; that explains effective collapse without a literal change to the universal wave function.
There are two serious objections. First, if every branch has the detailed functional organization of a detector and observer, what additional work do the particles do in making only the occupied branch a world with records? A many-worlds critic can question whether the particle ontology is doing explanatory work rather than merely labelling one branch. The Bohmian response is that ordinary objects and records are constituted by actual configurations, not by every formally similar pattern in . Assessing that response requires an account of material constitution; counting amplitudes alone settles neither side.
Second, the universal wave function affects particles without the usual reciprocal influence of particles on it. This can look unlike a material field. A nomological interpretation of the wave function treats its role more like law than additional matter, although time-dependent effective subsystem waves make that position philosophically demanding.[ref]
Neither objection licenses assigning ordinary particle energy separately to every empty component as though it contained a second classical apparatus. Such an assignment would require further assumptions. The definite mathematical theory reproduces standard equilibrium statistics while retaining empty components. The problem is the cost and intelligibility of that ontology, especially relative to alternatives, not an unexplained failure of its specified guidance equation.
For the oppositely travelling components, write and after removing the common time phase. Then . In the interference terms contribute no imaginary part, so
A standing-wave interference pattern can therefore change density without adding a cross-term to the net current.
For barrier-top scattering, the stationary equation inside the barrier is . Thus its spatial part is linear: . Take incident amplitude one, and . Continuity of the spatial wave function and its derivative at the two finite steps gives
Solving these four equations yields
The equal exterior velocities mean flux ratios are squared amplitude ratios:
This is the threshold limit of finite-barrier scattering, not exponential decay with a nonzero interior decay constant.
The interior constants are
Since ,
It is positive and independent of position, matching the transmitted flux. Although the classical kinetic energy inside is zero, the complex linear quantum wave function carries conserved probability current; its phase varies with position.
Figure 1.
At the barrier-top energy the complex amplitude is linear while its probability current remains constant
. Source. Details.
identical particles have the same intrinsic properties, including mass, charge and spin. Their labels are bookkeeping devices rather than measurable identities. When wave functions overlap, quantum mechanics supplies no observable trajectory that permanently identifies “particle 1” independently of “particle 2”. A many-particle description must therefore treat exchange consistently. The Hamiltonian and physical observables are invariant under relabelling identical particles.
For two particles, the exchange operator obeys . The usual three-dimensional symmetrization postulate selects states symmetric under exchange for bosons and antisymmetric for fermions. In relativistic quantum theory, the Spin-statistics theorem associates integer spin with bosons and half-integer spin with fermions under its locality and positivity assumptions; the connection is not derived merely from in elementary nonrelativistic mechanics. For orthonormal one-particle states and , the two-particle spatial combinations are
The sign concerns the complete state, including spin and other internal variables, not necessarily the spatial factor alone.
The antisymmetric state vanishes if both fermions occupy the same one-particle state: this is the Pauli exclusion principle. For fermions in orthonormal spin orbitals , the normalized Slater determinant supplies the required exchange sign. Bosonic occupation has no corresponding upper limit, and symmetrized products describe many bosons. Fermionic occupation is zero or one per complete one-particle state; two electrons can still occupy the same spatial atomic orbital with opposite spin.
Spin and spatial symmetry are linked. A two-electron spin singlet state is antisymmetric, so its spatial factor must be symmetric. A symmetric triplet state instead requires an antisymmetric spatial factor. Exchange then changes pair correlations and the energies available when interactions depend on particle separation, even without introducing a new classical “exchange force”. It contributes to atomic shell structure, bonding and magnetic behaviour.
Figure 1.
Symmetric and antisymmetric two-particle spatial states have different joint-position correlations
. Source. Details.
In an occupation-number description, bosonic creation and annihilation operators satisfy , whereas fermionic operators satisfy and anticommute among themselves. The latter imply . This formalism implements the same exchange constraints without assigning permanent particle labels.
For noninteracting particles in thermal equilibrium, the mean occupation of a state of energy is
Fermi-Dirac statistics produces a filled Fermi sea at zero temperature and degeneracy pressure; Bose-Einstein statistics permits macroscopic Ground state occupation when the density of available states and thermodynamic conditions allow condensation. The different statistics therefore have measurable collective consequences, not just different notation.[ref]
When phase-space occupation is dilute, both distributions reduce to the Maxwell–Boltzmann statistics limit. Likewise, well-separated localized wave packets have negligible overlap and small exchange corrections, so effectively labelled classical trajectories can become a useful approximation. In two spatial dimensions more general exchange statistics are possible; the standard boson/fermion account here concerns ordinary three-dimensional particles and their total states.
In the time-independent Schrödinger equation, is the kinetic-energy operator acting on the state, its potential-energy contribution and the energy-eigenvalue term. The complex wave function is a probability amplitude: the Born rule assigns density for a normalized single-particle state. A plane-wave beam uses an incident-amplitude or flux normalization instead of a unit-integral wavefunction.
For a free plane wave, , so and its time factor obeys . Above the upward potential step, the bounded subthreshold solution has decay constant
For the PDF's unit incident amplitude, write on the left and on the right, omitting their common time phase. Continuity of both and at this finite step gives
Consequently
The evanescent amplitude versus transmitted current distinction matters: the second quantity is a squared amplitude, not a transmitted-flux probability. The decaying tail has zero probability current, despite nonzero density and the possible value .
Using and the electron mass and elementary charge specified above, the energy difference1eV gives . Thus at the stated position
This is relative to the unit incident density in the source figure, not an absolute normalized probability per metre. [ref]
Here is potential energy, is the stationary state's total energy, and is its spatial probability amplitude. The wave function has time factor ; is a position probability density for a normalized state. Stationary plane-wave scattering solutions are idealizations described by probability-current ratios, or components of a normalizable wave packet.
Above an upward step, a classical particle is fully transmitted with reduced speed. Quantum mechanically the abrupt wavelength change generally produces both reflected and transmitted amplitudes even though both regions are classically allowed. With unit incident amplitude and no wave arriving from the right,
The reflected component has signed wavevector and the same wavenumber magnitude as the incident component. A finite potential jump with constant mass and no delta interaction requires continuity of and :
The probability current gives probability reflection and transmission and , with . The amplitude may exceed one without violating conservation of total probability, because the transmitted group speed is smaller.
For , let . The bounded right-hand solution is , an evanescent wave tail with amplitude penetration scale . It carries no outgoing probability current in an infinitely long forbidden half-space, so reflection is total. This differs from transmission through a finite potential barrier, where a second boundary can support quantum tunnelling.
For a downward step, replace by . The same boundary conditions give a nonzero negative , corresponding to a reflected amplitude with phase reversal. If reflection were forbidden, would require both and , impossible for a nonzero step. Reflection therefore does not require a classical potential barrier; the discontinuity of wavenumber mismatch is enough.
For the specified electron energies, , so the requested amplitude coefficients are
The corresponding probability coefficients are approximately and . electron mass cancels from these ratios; there is no need to invent an additional material effective mass.
The ground-state energies are and , respectively. For a particle in a box of width , the wave function must vanish at both infinite walls. The lowest nonzero standing wave therefore has wavelength , giving the kinetic energy eigenvalue
Using the electron mass and converting joules to electronvolts gives
Doubling the confinement length divides this kinetic energy by four. A zero potential energy inside the box does not imply zero total energy: a constant, zero-momentum wave function cannot satisfy both wall conditions unless it vanishes everywhere. This is a concrete manifestation of quantum confinement and the Heisenberg uncertainty principle.
The one-dimensional time-independent Schrödinger equation is
For the infinite square well on , inside and the wave function is zero outside. The boundary conditions are . For positive integer , satisfies them, and . Substitution gives the energy levels
Normalization requires . The requested real positive constant is therefore
For the second isotope, consecutive gap differences are approximately MHz, but the smallest gap is about times this value. The higher- values for the five gaps are therefore . The six levels run from to . The electronic state is unchanged, so as in d; the endpoints give
For example, fitting to with the same equal-weight convention gives MHz. The nearly constant gap differences are compatible with the same leading model of the magnetic-dipole hyperfine interaction and the printed rounding.
Within an unchanged electronic state, and neglecting isotope-dependent changes of the electronic wave function and hyperfine anomaly, . Thus
or the reciprocal magnitude is about . This is a ratio of nuclear g-factors, not directly a ratio of nuclear magnetic moments: , so the moment-magnitude ratio would additionally contain and is about .
The positive transition frequencies do not fix either coupling sign. If a signed ratio is required, equal signs must be supplied as an additional assumption or determined by an ordering-sensitive measurement. No absolute nuclear g-factor is inferred without a calibrated electronic field coefficient. The second ladder in b illustrates the inferred half-integer sequence without turning that sign ambiguity into a falsely measured ordering.
Take parity symmetry in quantum field theory to preserve free time evolution, so the conjugate initial data reflect along with the field. Substitute the mode expansion of a free field into the defining parity operator relation and change the integration variable . Since , coefficient comparison, or comparison of both the field and its conjugate initial data, gives the intrinsic-even convention
Here parity symmetry in quantum field theory reverses spatial momentum and preserves time. The creation operator formulas follow by taking Hermitian adjoints. The given field relation has fixed the intrinsic phase to .
For one real scalar field, let act on one-particle wave functions by . It is a Hermitian operator with , so is the orthogonal projector onto odd wave functions. Its second quantization is the odd-mode number operator
The full momentum space is integrated in both expressions; the factors prevent double counting the pairs . A suitable parity operator is
For verification, . Thus the even combination is unchanged, and the odd combination changes sign under conjugation by , exactly reproducing . The same holds for creation operators. The odd-mode number operator has integer eigenvalues, so , and . Multiplying by an overall constant phase would not change its conjugation action. This is parity as odd-mode number parity.
In the scalar two-particle model, the fermionic exchange symmetry and bosonic exchange symmetry conditions are
For particles with spin, exchange acts on all their labels; the stated coordinate-only model fixes the spatial symmetry sector. Define . Equal masses make the kinetic operator invariant under exchange, and the assumed potential is invariant too, so . Therefore implies . The components
are respectively symmetric and antisymmetric eigenfunctions whenever nonzero. Their sum is the original nonzero quantum state, so at least one is nonzero. An eigenspace may contain both sectors if the energy is degenerate.
The normal modes of two harmonically coupled identical particles provide the separation. The coordinate transformation to is orthogonal, with absolute Jacobian one. Exchanging the particles leaves fixed and changes to . Hence bosonic quantum states are even and fermionic quantum states odd in the relative coordinate. The chain rule gives
Also and . Combining the two trapping potentials with the interaction therefore gives the separated Hamiltonian operator
These are independent oscillators of frequencies and . Both frequencies are positive in the stated range. Their product eigenstates have energies
The parity of an oscillator level is . Consequently bosons allow even and fermions odd , while is unrestricted.
The bosonic ground state uses , and the fermionic ground state uses . The oscillator length in the relative mode is
Using the normalized Gaussian for each mode gives
Indeed by differentiating the Gaussian, so it is the normalized odd first excited state. The ground energies are and .
Choose normalized factors for a normalized product quantum state; this can always be done by reciprocal rescaling of its factors. Since the coordinate Jacobian has absolute value one,
If the supplied factors are not individually normalized, the general expression has an additional factor ; choosing normalized factors is the convention in the displayed formula.
For the ground states, put . Then and . Thus, for differentiable with suitable Gaussian integrability and vanishing boundary terms,
This is integration by parts, and every polynomial satisfies the required endpoint conditions. The relative-coordinate moments in even and odd oscillator ground states follow by substituting :
For odd , both expectations vanish by parity, so this is an equality of zeros rather than an increase in signed mean separation. The even moments show the physical difference. In particular and , so the fermionic root-mean-square particle separation is times the bosonic value. The fermionic probability density vanishes at coincidence , expressing Pauli exclusion principle in this spatial symmetry sector. This is an exchange constraint on the wave function, not an additional classical force. As approaches , the relative oscillator softens and its width diverges; the excluded endpoint would no longer have this normalizable oscillator ground state.
For identical fermions the full wave function is antisymmetric:
In particular it vanishes when the two complete one-particle arguments coincide. Applying this to the position-only model assumes that no independent spin symmetry changes the required spatial parity.
For identical bosons the full wave function is symmetric under the particle exchange operator:
In the given position-only model this is the spatial parity under the particle exchange operator. If additional spin variables are present, the symmetry applies to the complete state including them.
The orthogonal coordinate change preserves the kinetic energy and takes the particle exchange operator to . Thus the spatial wave function is even in for bosons and odd for fermions.
The total potential energy and Hamiltonian separate as
Assume so the relative oscillator is confining. Its spectrum is , and its parity in the two-particle relative coordinate is .
The lowest boson state uses , whereas the lowest fermion state uses , . Therefore
For an oscillator of frequency , . Since ,
The factor three is the extra relative-coordinate excitation required by antisymmetry. The supplied oscillator wave functions correspond to kinetic energy ; the extra factor printed before that term in the hint is inconsistent with those functions and the stated energies. The calculation uses the conventional kinetic energy dictated by particle mass .
For a normalized state the position probability density is . With a real potential, the Time-dependent Schrödinger equation and its conjugate give
This is the continuity equation . Integrating over the line yields under the stated decay. A real-valued wave function gives , so its instantaneous current vanishes; the same is true up to a spatially constant phase.
Put . Complete the square in the overlap:
The Gaussian integral is , also for complex displacement by analytic continuation. Hence
The product of square roots fixes the branch convention implicit in the printed square root of . At the norm squared is , which becomes one for the stated . In particular
The centre is and its position variance is .
When differentiating in , keep fixed on the chosen constant-radius orbit. Then
Substitution gives . Since , the chain rule gives . Continue continuously along the time orbit; its sign after a full orbit supplies the usual oscillator zero-point phase. The case of zero displacement is obtained as a limit, since the separate printed and expressions are singular at .
Finally . For real nonzero this is
The printed final current omits the Gaussian normalization factor . Retaining that factor is required both by the normalized state and by the continuity equation. This current is the Gaussian probability density times the velocity of its moving centre, .
For bound hydrogen atom states, the radial wave function is regular and finite at the origin, and decays at infinity so that . Equivalently the reduced radial function has and is square integrable on the half-line. A singular branch is not an admissible regular Coulomb eigenfunction; its derivative energy diverges and it introduces an origin singularity absent from the Hamiltonian domain.
Write , so . Multiplication by turns the radial equation into
Since and , substitution gives
Choose to remove the constant coefficient and factor out the decaying exponential. For , the lowest power has coefficient , so . The behavior would give , violating regularity; choose .
Matching successive powers then yields
A nonterminating tail eventually has one sign, and its coefficient ratio is asymptotic to . More precisely, for any , the ratio is eventually at least . Its same-sign tail consequently dominates a positive constant times the exponential series of rate , up to finitely many terms. The factor cannot then make decay. Thus square integrability forces termination.
If the last nonzero term is , termination requires . Therefore
For , , is proportional to , and
The radial volume element is , so the mean radius is the ratio of the radius-weighted norm integral to the norm integral. With , the common angular and normalization factors cancel:
In the ground state .
A permutation of identical particles multiplies a fermionic state by ; a bosonic state is invariant. Particle exchange keeps fixed and sends to . Consequently the spatial wave function is even in for the stated bosonic model and odd for the stated fermionic model. This assumes no additional spin variables in the supplied model; exchange of full particle states is the underlying rule.
The coordinate change is orthogonal, with unit Jacobian determinant, so the kinetic energy becomes
on the square with zero boundary values. Product states have energy
The one-dimensional eigenfunction is even for odd and odd for even . Bosons therefore allow odd , whereas fermions require even . The lowest bosonic pair is ; its first excited pair is . The lowest fermionic pair is . Thus
The printed bosonic ground-energy claim omits the factor ; the value above follows directly from the supplied one-dimensional spectrum and kinetic energy.
The separation is . The normalized factor integrates out. In the fermionic ground state, the relative-coordinate density is , giving
For the bosonic ground state it is , giving
Both integrals follow by replacing squared trigonometric functions by their constant and cosine parts and integrating once by parts. Hence
and the difference is . The odd relative state vanishes at coincidence, which increases its expected separation.
A permutation of identical particles multiplies a fermionic state by ; a bosonic state is invariant. Particle exchange keeps fixed and sends to . Consequently the spatial wave function is even in for the stated bosonic model and odd for the stated fermionic model. This assumes no additional spin variables in the supplied model; exchange of full particle states is the underlying rule.
The coordinate change is orthogonal, with unit Jacobian determinant, so the kinetic energy becomes
on the square with zero boundary values. Product states have energy
The one-dimensional eigenfunction is even for odd and odd for even . Bosons therefore allow odd , whereas fermions require even . The lowest bosonic pair is ; its first excited pair is . The lowest fermionic pair is . Thus
The printed bosonic ground-energy claim omits the factor ; the value above follows directly from the supplied one-dimensional spectrum and kinetic energy.
The separation is . The normalized factor integrates out. In the fermionic ground state, the relative-coordinate density is , giving
For the bosonic ground state it is , giving
Both integrals follow by replacing squared trigonometric functions by their constant and cosine parts and integrating once by parts. Hence
and the difference is . The odd relative state vanishes at coincidence, which increases its expected separation.
In the scalar two-particle model, the fermionic exchange symmetry and bosonic exchange symmetry conditions are
For particles with spin, exchange acts on all their labels; the stated coordinate-only model fixes the spatial symmetry sector. Define . Equal masses make the kinetic operator invariant under exchange, and the assumed potential is invariant too, so . Therefore implies . The components
are respectively symmetric and antisymmetric eigenfunctions whenever nonzero. Their sum is the original nonzero quantum state, so at least one is nonzero. An eigenspace may contain both sectors if the energy is degenerate.
The normal modes of two harmonically coupled identical particles provide the separation. The coordinate transformation to is orthogonal, with absolute Jacobian one. Exchanging the particles leaves fixed and changes to . Hence bosonic quantum states are even and fermionic quantum states odd in the relative coordinate. The chain rule gives
Also and . Combining the two trapping potentials with the interaction therefore gives the separated Hamiltonian operator
These are independent oscillators of frequencies and . Both frequencies are positive in the stated range. Their product eigenstates have energies
The parity of an oscillator level is . Consequently bosons allow even and fermions odd , while is unrestricted.
The bosonic ground state uses , and the fermionic ground state uses . The oscillator length in the relative mode is
Using the normalized Gaussian for each mode gives
Indeed by differentiating the Gaussian, so it is the normalized odd first excited state. The ground energies are and .
Choose normalized factors for a normalized product quantum state; this can always be done by reciprocal rescaling of its factors. Since the coordinate Jacobian has absolute value one,
If the supplied factors are not individually normalized, the general expression has an additional factor ; choosing normalized factors is the convention in the displayed formula.
For the ground states, put . Then and . Thus, for differentiable with suitable Gaussian integrability and vanishing boundary terms,
This is integration by parts, and every polynomial satisfies the required endpoint conditions. The relative-coordinate moments in even and odd oscillator ground states follow by substituting :
For odd , both expectations vanish by parity, so this is an equality of zeros rather than an increase in signed mean separation. The even moments show the physical difference. In particular and , so the fermionic root-mean-square particle separation is times the bosonic value. The fermionic probability density vanishes at coincidence , expressing Pauli exclusion principle in this spatial symmetry sector. This is an exchange constraint on the wave function, not an additional classical force. As approaches , the relative oscillator softens and its width diverges; the excluded endpoint would no longer have this normalizable oscillator ground state.