Put . Logarithmic differentiation gives
Therefore
so the appropriate potential energy is . The wavefunction never vanishes, so no other time-independent potential gives the same equation for this precise phase choice.
The squared modulus is
It changes width and height with time, so this is not a stationary state. Its integral is the constant ; multiplying the wavefunction by normalizes it. The normalized Gaussian wave packet stays centred at zero but spreads, with position variance . It represents a localized free particle with a spread of momenta, rather than an energy eigenstate. Free Gaussian wave packet evolution broadens the packet even with zero mean momentum.
The GRW model supports robust macroscopic locations, but not universal exact point locations supplied by a list of particle coordinates. It modifies the wave function rather than adding the actual particle configuration of Bohmian mechanics. A claim that objects are located must therefore specify both the collapse dynamics and the proposed ontology. “Localized” is not synonymous with “has compact support” or “is a position eigenstate”.
Using one standard convention, a collapse on constituent with centre acts by the GRW localization operator
Its centre has probability density . Independent constituent collapse rates give a body with relevant constituents an approximately localization rate. Correlations within a macroscopic superposition allow one constituent's localization to suppress an entire alternative configuration: the GRW amplification mechanism.[ref]
For two narrow packets separated by , a hit near one multiplies the other's relative amplitude approximately by . This explains a sharply concentrated pointer reading. Yet the Gaussian function is never zero at finite separation. A nonzero distant tail therefore remains, and a finite-width packet is not an exact position eigenstate. This is the tails problem in collapse theories.[ref] The assertion that each hit produces exact bounded support is false; an initially everywhere nonzero Gaussian wave packet is a concrete counterexample.
A GRW matter-density ontology specifies matter in ordinary SPACE by
A macroscopic body's density is concentrated near one region, making its approximate location physically meaningful. The formula is the proposed matter field itself, not a probability distribution over a hidden set of exact point positions. Residual tails correspond to very low density elsewhere, motivating the further question of why such tails do not constitute ordinary rival objects. Concentration provides a practical account of records; a full ontology must also address that structured-tail objection.
A GRW flash ontology instead identifies the primitive events with collapse centres in spacetime. Individual flashes have definite locations and times, while a persisting macroscopic object is a dense pattern of flashes. There is no continuous particle trajectory between them. An isolated microscopic system may flash rarely, so the model must not be read as assigning it an ordinary point location at every instant. Matter-density and flash accounts express distinct ontologies with related empirical predictions.[ref]
An approximate-location rule, such as requiring almost all matter density inside a suitably sized region, is meaningful if its tolerance and scale are stated. It does not turn the small tail into zero, and arbitrary sharp property thresholds can create artificial counting puzzles. The defensible result is emergent, highly stable macroscopic localization, or exact locations of flashes, depending on the ontology. Exact Bohmian positions are a different theoretical commitment and cannot be attributed to the collapse equation alone.
A localized quantum particle is described by a wave packet: a coherent quantum superposition of momentum components, rather than a single delocalized plane wave. The Born rule assigns the position probability density . A useful normalized initial Gaussian wave packet is
Its mean position and momentum are and , and its standard deviations satisfy , . This unchirped initial state saturates the Heisenberg uncertainty relation.
For a free particle, the momentum-space Schrödinger equation multiplies each amplitude by . The inverse Fourier transform is a Gaussian integral, giving the free Gaussian wave packet evolution
One way to obtain this expression explicitly is to write the initial momentum amplitude as proportional to , multiply by the free evolution phase and complete the square in . Taking the modulus squared removes the position-dependent phase and leaves
The expectation value follows uniform classical motion while the density spreads. A narrower initial packet needs a wider momentum distribution and therefore spreads more rapidly. At long times its width grows linearly with ; its variance grows quadratically, unlike the linear variance growth of ordinary diffusion. The probability continuity equation preserves total probability despite the decreasing peak height. The mean momentum and the mean energy remain constant, since the momentum populations do not change.
The same result has a wave interpretation. For a free nonrelativistic particle , so the group velocity is , equal to at the packet center. The curvature makes different momentum components travel with different group velocities, producing dispersion. The phase velocity is for this energy convention and does not describe the motion of the envelope. The evolving quadratic spatial phase is essential to the spreading; a density snapshot alone does not specify the full quantum state.
A freely spreading Gaussian wave packet remains a pure Gaussian state but develops position-momentum correlations. Its symmetrized covariance is . Consequently
Equality in the Robertson-Schrödinger uncertainty relation survives at every time, whereas the simpler product equals only initially. A chirped Gaussian wave packet can already have . More generally, gives . A sufficiently negative initial covariance permits initial contraction before subsequent spreading.
Figure 1.
Free Gaussian spreading, constant coherent-state width and the breathing width of a squeezed state
. Details.
A potential energy changes this behavior. A quadratic potential preserves the Gaussian form because the evolution of position and momentum is linear. For a quantum harmonic oscillator, the Gaussian covariance in a harmonic potential is
A coherent state has , and zero centered covariance. Its width is constant, and its center follows . In oscillator notation its coherent state label rotates as , apart from the common ground-energy phase of the state vector. A squeezed coherent state generally has a width that breathes at . Thus being Gaussian does not imply either free spreading or constant width; the Hamiltonian operator and initial covariance determine the result.
For a general nonquadratic potential, different parts of the packet experience different nonlinear forces, and an initially Gaussian profile need not remain Gaussian. Ehrenfest theorem gives the exact relations and . Replacing the last mean by is exact for a quadratic potential and approximate for a packet narrow compared with the scale on which the force varies. This connects packet motion to the correspondence principle without treating a quantum state as a definite classical path. Scattering can split a packet into reflected and transmitted components, while confinement can produce oscillations and interference between discrete energy eigenstates. Their relative phases evolve even when their populations remain constant. Localization, spreading and coherent motion all follow from the same linear unitary time evolution.
The natural interpretation is the fundamental transverse normal mode, . Indeed
so only propagates over the main part of the given Gaussian distribution of angular frequencies. Its centre travels at
Differentiate the group velocity with respect to angular frequency:
At the centre this is . The one-standard-deviation group velocity spread is consequently
Over the central travel time the associated spatial shear is
Equivalently, use representative components at , for which and . Their separation after is
At a fixed detector away, the corresponding arrival-time difference is .
The paper asks for an estimate of spreading, without specifying a spatial width convention or the initial spectral phase. A Gaussian distribution has no finite full range. Thus is a one-standard-deviation dispersive contribution, whereas about describes the separation of the representatives. Neither is automatically the increase in the packet's original width. For an initially unchirped Gaussian wave packet in the narrow-band approximation, independent initial and dispersive contributions combine as ; a packet with an initial chirp can behave differently. The modal assumption and width convention must accompany the numerical answer.
Both first moments vanish, so the standard deviations follow directly from the second moments:
The Heisenberg uncertainty relation is satisfied and is not saturated. A Gaussian wave packet can saturate the position-momentum bound; this compactly supported polynomial is not Gaussian.
The normalized initial Gaussian wave packet is real and even, so the integrands for and are odd. Its second position moment is the elementary Gaussian integral . Since and ,
Moreover and . Thus
Quantum harmonic oscillator Created 2026-09-24 Updated 2026-10-13
The quantum harmonic oscillator describes a particle of mass in a quadratic potential energy with angular frequency . In the position representation its Quantum Hamiltonian is
The model is useful far beyond a literal spring: a Taylor series about a stable equilibrium with positive potential curvature has this quadratic term as its leading nonconstant approximation. The pure quartic oscillator has zero curvature at its minimum, so that particular approximation does not apply. [ref].
In units with , creation and annihilation operators give
Here is the number operator. Restoring the reduced Planck constant, the energy eigenvalues are for : neighboring levels are separated by , and the ground-state energy is the zero-point energy .
With length scale and dimensionless coordinate , the ground state is the normalized Gaussian wave packet
For , this is precisely the trial wavefunction in the quantum variational principle. The balance of kinetic energy and potential energy therefore recovers the exact ground-state energy using only Gaussian integrals.
Quantum variational principle Created 2026-10-03 Updated 2026-10-14
The quantum variational principle estimates the lowest energy by choosing a trial wavefunction and computing its expectation value, without solving the Schrödinger equation exactly. Minimizing over a manageable family of trial wavefunctions gives an upper bound; enlarging that family can improve the estimate.
Let be a self-adjoint operator bounded below, and put . For in its operator domain ,
Trial states may also lie in the form domain , where . Then use the quadratic form of a positive quantum Hamiltonian shifted by :
This extends the expectation value beyond without requiring to exist. Equality for a normalized state holds exactly in the energy eigenspace at , if that eigenvalue exists. [ref].
A short proof makes the direction of the bound memorable. If has a complete orthonormal basis of energy eigenstates , write
The expectation value is a weighted average of allowed energies:
Equality requires every occupied energy eigenstate to have energy . [ref]. With continuous spectrum, the spectral theorem replaces the sum by an integral against the scalar spectral measure; the same nonnegative-difference argument works.
For an unnormalized trial wavefunction, divide by its squared norm: is the form version of the Rayleigh quotient. The quantum variational principle supplies an upper bound on ; it supplies neither a lower bound on nor, by itself, a guarantee that a bound state exists.
To see the method at work, choose the real Gaussian wave packet
The Gaussian integral gives wavefunction normalization and the two useful position averages:
Thus its position probability density is a normal distribution with standard deviation : large means narrow confinement. In units , integration by parts and give the kinetic energy
These smooth, rapidly decaying trial wavefunctions belong to the operator domains of both polynomial Quantum Hamiltonians below. Narrowing the Gaussian wave packet increases its kinetic energy; widening it increases the potential energy in either confining potential.
For the quantum harmonic oscillator, take :
The trial values at are respectively . The derivative changes sign at , where kinetic energy and potential energy are both . Here the trial family contains an exact ground state:
The zero-point energy remains positive even at the best width. For comparison, removing the potential leaves as : the free Quantum Hamiltonian on the line has spectral bottom zero but no normalizable zero-energy energy eigenstate. A sequence of ever broader trial wavefunctions need not approach an actual bound state.
Figure 1.
Gaussian trial widths and the kinetic–potential energy tradeoff
. Normalized Gaussian wave packets with have different position probability densities. The bars show the competing kinetic energy and harmonic potential energy ; their sum is smallest at . Details.
For a pure quartic oscillator, keep and set the quartic coefficient to one:
The trial values at are respectively . Setting the derivative to zero gives
The sign of this derivative changes from negative to positive at the optimum, giving the smallest energy within this Gaussian wave packet family:
At this width, and , so , as expected from the power-law quantum virial theorem for a scale-optimized trial wavefunction. The number is an upper estimate of the ground-state energy, not an asserted exact energy eigenvalue. Indeed,
depends on for every , so no member of this trial family is an energy eigenstate. Nevertheless, a single width parameter already turns a difficult Schrödinger equation into a useful, explicit energy estimate.
Figure 2.
Energy landscapes for harmonic and quartic Gaussian trials
. The quantum harmonic oscillator reaches its exact ground-state energy at . The pure quartic oscillator has a Gaussian wave packet optimum at giving only the upper bound . Dashed curves separate kinetic energy and potential energy. Details.
Robertson uncertainty principle Created 2026-09-24 Updated 2026-10-20
The Robertson uncertainty principle bounds the product of the standard deviations of two observables in the same quantum state. For a normalized pure state and self-adjoint operators , write their expectation values and variances as
and similarly for , where is the identity operator. The quantum uncertainties then obey
Here is the commutator, and every expectation value is taken in the same . These standard deviations describe the outcome spreads when and are measured on separate ensembles prepared in that same quantum state. They are not instrument errors or the change caused by a preceding measurement in quantum mechanics. The definition and bound are developed in Barton Zwiebach, Uncertainty Principle and Compatible Observables, 2013, §§1–2.
For unbounded self-adjoint operators, the operator domains matter. Assume , so that both second moments are finite; here means the second moment of the measured probability distribution, equal to , even if is undefined. To use the displayed operator commutator, also assume . Without these product-domain assumptions, the always meaningful version on uses as its right-hand side. These qualifications prevent a formal commutator identity from being applied to an inadmissible quantum state.
The brief proof is Cauchy-Schwarz inequality applied to and :
The final equality uses self-adjointness and the product-domain assumptions. Thus a geometric statement about inner products becomes a constraint on experimental spreads. For nonzero , equality requires to be proportional with a purely imaginary proportionality factor; the equality case of the Heisenberg uncertainty relation turns this condition into a Gaussian wave packet.
The bound is useful because computing a single commutator and its expectation value can constrain two entire outcome distributions. Its right-hand side depends on the quantum state: the spin uncertainty state dependence example gives both a positive saturated bound and a zero bound for the same noncommuting observables. Robertson's paper, titled "The Uncertainty Principle", appeared on 1 July 1929 in Physical Review, volume 34, starting at page 163. [ref].
Schrödinger equation Created 2026-09-24 Updated 2026-10-20
The Schrödinger equation governs the evolution of a quantum state between instances of measurement in quantum mechanics. For a particle of mass moving in a prescribed real potential energy , its position wavefunction satisfies
Here is the Planck constant, is the Laplace operator, and the Hamiltonian operator combines kinetic energy with potential energy. Specify an initial wavefunction, an operator domain and any boundary conditions: the differential expression alone does not completely specify a physical model. The Born rule then turns the evolving probability amplitude into position probabilities through and integration over a region.
This scalar Schrödinger operator describes nonrelativistic quantum mechanics, with relevant momenta much smaller than , where is the speed of light. It omits the constant rest energy , assumes fixed particle number, and suppresses active spin dynamics and coupling to a magnetic vector potential. A free electron with a fixed, uncoupled spin state can still use this scalar model; the Pauli equation supplies the additional spin components and magnetic coupling when needed. A prescribed real describes a closed evolution model or an approximation in which the source of that potential energy is treated as external. Dissipation, particle creation and a general measurement in quantum mechanics need further modelling.
For a time-independent self-adjoint Hamiltonian operator, the spectral theorem defines unitary time evolution by
Initial states in the operator domain satisfy the differential equation in the Hilbert space sense; the unitary operator also evolves all other states in that Hilbert space. Time-dependent Hamiltonian operators need suitable regularity and domain hypotheses; their evolution is not generally the exponential of .
The local reason for conservation is especially tangible. Multiply the Schrödinger equation by the complex conjugate , subtract its conjugate multiplied by , and use the fact that the real potential energy terms cancel. For smooth wavefunctions this gives the probability continuity equation
The divergence theorem identifies the change in probability inside a region with minus the outward probability current. Decay at infinity, or appropriate probability-conserving boundary conditions, preserves wavefunction normalization. Conservation of total probability does not require a stationary probability density.
Its usefulness comes from being a linear partial differential equation: quantum superpositions of solutions are solutions, while different quantum phases can produce interference. For a fixed Hamiltonian operator, the time-independent Schrödinger equation finds the energy eigenstates; each then acquires just its own time-dependent quantum phase. One eigenvalue problem therefore organizes both allowed energies and the evolution of arbitrary quantum superpositions, with a continuous spectrum included when required. This is the attractive bridge between a spatial wave problem and a clock: energy differences determine the rate at which relative quantum phases change. It makes atomic spectral lines, quantum tunnelling and quantum scattering calculable from a specified model.
For an explicit example without confinement, take a free particle on the whole line, , initially in a Gaussian wave packet with position standard deviation and mean momentum :
Unlike a single plane wave, this wavefunction satisfies wavefunction normalization. A Fourier transform writes it as a quantum superposition of plane waves. Their dispersion relation is , so each Fourier component gains the factor . Completing the square in the resulting Gaussian integral gives
with the square root chosen continuously from . Its position probability density and standard deviation are
The centre moves with the group velocity and the wave packet spreads because different wavenumbers have different group velocities. This is unitary wave dispersion, without a random force or loss of probability. Independently, the initial momentum standard deviation is and the initial position-momentum covariance vanishes. Free evolution gives in the Heisenberg picture, hence
This independently checks the spreading law. The Gaussian integral of is one at every time; the decreasing peak is exactly offset by increasing width.
As hypothetical preparation parameters, choose a free electron with and . Using and gives
The rounded electron mass uses the 2022 adjustment.[ref] At , the centre is and ; at , they are and . The mean kinetic energy stays
It includes both mean motion and the momentum spread. The mean speed is about and the momentum spread corresponds to , making the nonrelativistic model accurate for the appreciably occupied part of this ideal Gaussian distribution. Its mathematically unbounded tails are an idealization.
Figure 1. . Calculated electron position probability densities for the hypothetical preparation above, together with the position standard deviation. The area under every curve is one. The dashed width curve illustrates pure translation without spreading for comparison; it is not the free solution. No curve represents measured data. Source. Details.
The complex width also explains where the probability goes. Write , where is the dimensionless quantum phase. Differentiating the exact solution gives
Relative to the moving centre, the probability current points outward on both sides for . This is a flux statement, not an additional assumption about individual particle trajectories. At and , the relative flow speeds are , about . The probability within a moving interval of half-width falls from about initially to at , exactly because of this outward probability current.
Figure 2.
Curved quantum phase drives spreading probability current
. The same free Gaussian wave packet at . With , removing the carrier and a spatially constant quantum phase leaves . The right panel plots : negative to the left, positive to the right. Flux through the marked moving interval endpoints explains its decreasing probability. Source. Details.
Historically, the 1926 Quantisierung als Eigenwertproblem (Erste Mitteilung) by Erwin Schrödinger began with the stationary hydrogen atom problem. It distinguishes discrete negative energies from a positive continuous spectrum; it is not yet the general evolution equation displayed above.[ref] The opening of quantum mechanics discusses its original vibrating-string analogy.
In the 12 December 1933 lecture The fundamental idea of wave mechanics by Erwin Schrödinger, Erwin Schrödinger explains the move from independent rays to waves through diffraction, then applies that analogy to atomic physics. His closing pages explicitly acknowledge the difficulty of reconciling wave language with observed particle tracks. This retrospective physical argument motivates the model; it does not derive the modern Born rule from classical mechanics.[ref]
The dates separate proposal from measurement: the de Broglie wavelength belongs to the 1924 thesis, the Schrödinger equation to 1926, and the Davisson–Germer experiment's electron diffraction beams to 1927.[ref] That linked experiment describes the prepared nickel crystal, evacuated electron gun, angular collector scans, retarding-potential selection and heating controls. Its measured collector electric current tests wave-dependent scattering probabilities; it neither measures a wavefunction directly nor establishes this equation uniquely.
Time-independent Schrödinger equation Created 2026-09-24 Updated 2026-10-20
For a time-independent Hamiltonian operator , an energy eigenstate satisfies
For a one-dimensional particle of mass in a real potential energy ,
Its full time-dependent wavefunction is : substitution in the Schrödinger equation checks both the minus sign and the factor . Thus an energy eigenstate has stationary position probability density, even though its quantum phase evolves. The boundary conditions and operator domain are part of this eigenvalue problem, selecting physically allowed solutions; the infinite square well and Radial Schrodinger equation for the hydrogen atom give complementary confined examples.
For a complete discrete orthonormal basis of energy eigenstates, linearity gives , with . Cross terms in oscillate at angular frequencies . A continuous spectrum requires the corresponding spectral integrals; for a free particle on the whole line, individual plane waves are generalized eigenfunctions, and a normalized Gaussian wave packet is their integral quantum superposition. Adding a spatially constant real to shifts every energy eigenvalue by and multiplies the full state by , a global phase; energy differences, position probabilities and probability current are unchanged.