Write , and . The Ising-chain Bogoliubov diagonalization uses the eigenvalues of the two-by-two Nambu matrix. For distinct paired modes , choose , and . The rotation satisfies , and its first transformed component is
The opposite signs of the paired angles preserve the canonical anticommutation relations. Self-paired spinless fermion modes have zero pairing and are treated directly as occupied or empty number levels; they do not require this paired-angle formula.
Both and occur in the sum, so the physical quasiparticle coefficient is , not . Consequently
This is the bulk quadratic result under the stipulated boundary simplification. Restoring the finite-chain parity sectors adjusts allowed modes and global excitation constraints. In the thermodynamic limit, the ground-state energy per site is
For , the bulk quasiparticle gap is . The requested three regimes are shown in the original sketch and discussed separately below.
Figure 1.
Ising-chain quasiparticle dispersion is flat at zero field and gapless at the critical field
. Details.
For the new periodic potential, let be the single-neighbour hopping rate and . It is determined by the barrier between adjacent minima; an unspecified periodic potential does not determine it from the preceding quartic potential's numerical action. The minimum spacing is now .
A path with right hops and left hops has . Summing the ordered-center weights and their direction choices gives
Insert the Fourier representation of a Kronecker delta. The two exponential series sum independently, yielding the dilute-hopping lattice propagator
The two Fourier-sign choices are equivalent by . Equivalently the integral is the modified Bessel function .
The original PDF prints here. Its sign is inconsistent with both the previous double-well result and the requested positive on-site oscillator energy. The correct factor is : a positive oscillator ground-state energy must decay under . The other factors and their derivation remain as displayed. The energy is the positive nearest-neighbour quantum tunnelling matrix-element magnitude, exponentially small relative to the local well scale in the semiclassical regime.
For a normalized trial state , expand it in an orthonormal basis of energy eigenstates, . The Rayleigh quotient is then
This is the quantum variational principle: minimizing the expectation of the Hamiltonian operator over any trial family gives an upper bound on the ground-state energy.
Let be the normalized ground-state wavefunction and form the normalized dilation . Its kinetic and potential expectations are
The true ground state makes the energy stationary at , so
Thus the power-law quantum virial theorem gives
If and , the two sides have opposite signs, because and . If , the scaled energy tends to as , since , so the Hamiltonian is not bounded below. If , the free particle has no normalizable ground state. Hence there is no normalizable ground state for .
For take the normalized Gaussian function
The Gaussian integral gives
and therefore
Its unique minimum satisfies
At this point , consistently with the virial theorem, and the resulting variational estimate is
Since and , the diagonal Hamiltonian is
with the identity on every unlisted qubit. Each computational-basis state is an eigenstate, and its eigenvalue is minus the cost of the corresponding cut. Part i supplies cost four, while the assumed bound rules out a lower energy. Thus the ground-state energy is . For the assignment , one ground state is
Its bitwise complement is another ground state, as are the basis states corresponding to the other maximum cuts.
Use the intended infinite square well interpretation: the kinetic-energy operator on with zero Dirichlet boundary conditions. The variational energy is most safely evaluated through its quadratic form,
The boundary term vanishes because the trial function vanishes at both walls. Its zero extension has derivative jumps there, so the form interpretation avoids assigning an unrestricted full-line second derivative to it.
Set and . Apart from an irrelevant factor , the trial function is . Polynomial integration gives
Hence the polynomial Rayleigh-Ritz trial space for a one-dimensional box has energy
The denominator is strictly positive, since . Differentiation gives
The two stationary parameters are . The smaller is a local maximum; the larger is a local minimum. Since the energy tends to as , the latter is the global minimum:
The upper stationary value is . Exact roots give a slightly tighter bound than substituting the rounded hints.
There are two qualifications to the printed operator and comparison energy. For the box with Dirichlet boundary conditions and width , the actual ground-state energy is
which lies just below the variational result. The coefficient printed as an exact eigenvalue is inconsistent with this box problem. Also, merely declaring the differential expression to be zero outside the interval does not impose hard-wall confinement or specify a self-adjoint domain. Read literally on a space admitting exterior-supported states, it has zero-energy states outside and does not have the claimed positive ground energy. The boxed variational result is the well-defined answer for the intended confined particle; a different full-line operator would need its domain and matching conditions supplied.
The infinite square well has energy levels , with . The longest absorption wavelength corresponds to the smallest available energy difference from the ground state, namely . Its position matrix element is nonzero because the two states have opposite parity about the center of the well, so it is an allowed electric dipole transition.
Equating the photon energy to this gap gives
Here the photon gap is approximately , not the ground-state energy itself.
Write and use the Heisenberg uncertainty relation separately in the three Cartesian directions. Since a position coordinate is confined to an interval of length , Popoviciu's inequality on variances gives . Consequently , and the expectation value of the kinetic energy obeys
Thus an uncertainty-based lower estimate is . This confinement energy bound from position support is not sharp: the exact ground-state energy for a cubical infinite square well is . The Heisenberg uncertainty principle establishes the confinement scale ; identifying a box size with an approximate position standard deviation changes its numerical coefficient.
For either number operator , positivity follows from
Hence an eigenvalue cannot be negative. The commutators
show that , if nonzero, is an eigenstate with eigenvalue . For a normalized eigenstate, .
If were not an integer, successive nonzero lowerings would eventually reach an eigenstate with . Its next lowering would have nonzero norm but a negative eigenvalue, contradicting positivity. The ladder must therefore terminate at , with . Conversely, successive applications of the creation operator never vanish, because
The normalizable bottom state exists explicitly: solving gives . Thus the number operator has exactly the non-negative integer eigenvalues.
Since the two coordinate number operators commute, use the product eigenstates . Their energy eigenvalues are
At fixed , the allowed pairs are . Thus the energy level has quantum degeneracy , and the ground-state energy is .
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.
Variational method Created 2026-09-24 Updated 2026-10-13
The variational method estimates the lowest energy of a Quantum Hamiltonian by choosing a trial wavefunction, evaluating its energy, and optimizing the trial parameters. It uses the quantum variational principle, the Rayleigh-Ritz variational principle expressed for quantum states. Let be a self-adjoint operator bounded below, and define , the bottom of its spectrum. Every normalized admissible trial quantum state satisfies
When is an energy eigenvalue, it is the ground-state energy; equality holds precisely for a quantum state in the ground-state subspace. Otherwise the bottom of the spectrum need not be attained by any normalizable wavefunction.
Admissible means finite form energy. Choose a real and use the shifted quadratic form of a positive quantum Hamiltonian:
The form domain is independent of this choice of and can be larger than the operator domain . For , is the usual expectation value; for nonzero unnormalized trials, divide by to get the form version of the Rayleigh quotient. See the quantum variational principle for the general domain justification.
Thus each trial energy is an upper bound on , and minimizing within a trial family improves that upper bound. It is not a lower bound on . The Gaussian variational energy for the harmonic oscillator reaches because its family contains an exact ground state. The exact variational solution in a two-level system makes the equality condition visible as vanishing excited-state probability amplitude.
The method is useful when solving the Schrödinger equation exactly is difficult but integrals of a trial wavefunction are tractable. Enlarging a trial family cannot increase its optimal energy, although an optimum need not be attained. The beauty is that even an imperfect trial gives a rigorous bound, provided its form energy is computed exactly: numerical integration errors require separate control. The Gaussian variational estimate for a solvable sextic potential provides an existing example beyond a quantum harmonic oscillator.
A short proof explains the direction. With a complete orthonormal eigenbasis of energy eigenstates, write , where and . Then
This is a weighted mean of nonnegative energy differences; equality requires every occupied energy eigenstate to have energy . [ref]. For continuous spectrum, the spectral theorem replaces the sum by an integral against the scalar spectral measure of , with the same nonnegative integrand. This proves the bound without assuming an orthonormal eigenbasis exists.
The historical setting predates quantum mechanics: Rayleigh's The Theory of Sound, volume I, appeared in 1877, and Ritz's paper on solving variational problems of mathematical physics appeared in volume 135 of Crelle's journal, catalogued by its publisher as 1909, pp. 1–61. These are original works associated with the Rayleigh-Ritz method; the dates here identify the publications rather than assert priority. [ref]; [ref].