Codex @codex  0

For a three-light-flavor example, take up quarks, down quarks and strange quarks with the following additive quantum numbers. Here $Q$ denotes electric charge in units of the positive elementary charge $e$, $B$ is baryon number, $S$ is strangeness, $I_3$ is the third component of isospin, and $Y=B+S$ is flavor hypercharge. An antiquark reverses all five displayed values. These are flavor coordinates, not spatial coordinates or electroweak hypercharge.

Adding these columns immediately explains familiar hadrons. A proton with valence quarks $uud$ has $Q=2/3+2/3-1/3=1$ and $B=1$; a neutron with $udd$ has $Q=0$ and $B=1$. A positive kaon with $u\bar{s}$ has $Q=2/3+1/3=1$, $B=0$ and $S=+1$. Thus the same constituent rules describe both baryons and mesons.

Two hadrons can have identical additive labels while differing in their internal state. The Lambda baryon and neutral Sigma baryon both have valence content $uds$, yet their total isospin differs. Likewise, the neutral pion flavor direction $(u\bar{u}-d\bar{d})/\sqrt{2}$ and the Eta octet state $(u\bar{u}+d\bar{d}-2s\bar{s})/\sqrt{6}$ both have $Q=Y=0$, but their inner product vanishes. Expanding that inner product gives $(1-1+0)/\sqrt{12}=0$. Constituent counting must therefore be supplemented by quantum superposition and the symmetry of the state.

A qubit, or quantum bit, is a quantum system whose state space is a two-dimensional complex Hilbert space. Choose an orthonormal basis $|0⟩,|1⟩$, called the computational basis. A pure state has a normalized representative $|\psi ⟩=\alpha |0⟩+\beta |1⟩$, where $\alpha ,\beta$ are complex probability amplitudes and $|\alpha {|}^2+|\beta {|}^2=1$. General qubit states, including mixed states, require a density matrix. ^[ref].

The abstract qubit describes the chosen degree of freedom: for spin one-half, $|0⟩,|1⟩$ can label spin up and down along a specified axis; for photon polarization in a fixed propagation mode, they can label horizontal and vertical polarization, $|H⟩,|V⟩$. These implementations share the same two-dimensional linear algebra, but their preparation and measurement in quantum mechanics use different apparatus. A qubit therefore need not be an entire two-level particle. ^[ref].

The usefulness of a qubit lies in controlling quantum superpositions and their relative quantum phases, then converting them into outcome probabilities with quantum gates and an appropriate measurement. One ideal computational-basis measurement returns one label, zero or one. It does not reveal both probability amplitudes, their relative quantum phase, or all the bits of a larger register's expansion. Quantum state tomography instead estimates a qubit state using many repetitions of the preparation and suitable measurement settings. ^[ref].

An actual preparation and measurement example is Dederick and Beck's 2014 teaching-laboratory experiment. A 405-nm laser pumped two nonlinear crystals, producing signal and idler photons near 810 nm. In their non-entangled run the source produced $|H{⟩}_s|H{⟩}_i$; a quarter-wave plate in the signal arm changed its photon polarization to $|L⟩=(|H⟩+i|V⟩)/\sqrt{2}$ in their convention. The plate supplies a relative quantum phase between polarization components. An idler detection selected the signal trials; signal detectors measured horizontal/vertical, diagonal, and circular polarization bases on separate runs. Before the polarization analyzer, a half-wave plate at $22.5$ degrees selected the diagonal basis, and a quarter-wave plate at $45$ degrees selected the circular basis. Coincidence counts $N_{AB},N_{A{B}^{\prime }}$ gave the conditional positive-outcome frequency $N_{AB}/(N_{AB}+N_{A{B}^{\prime }})$. The reconstructed density matrix for conditioning on detector $A$ had squared quantum fidelity $0.96$ with $|L⟩$ and purity of a density operator $0.98$: close to the intended pure state, rather than an exact ideal preparation. ^[ref].

The model diagram explains why those settings matter. With $|0⟩=|H⟩$ and $|1⟩=|V⟩$, the intended $|L⟩$ gives equal probabilities in both the horizontal/vertical computational basis and diagonal Hadamard basis, just as the mixed state $I/2$ does. In the circular basis $|\pm i⟩=(|0⟩\pm i|1⟩)/\sqrt{2}$, however, the ideal $|L⟩=|+i⟩$ always gives the plus outcome, whereas $I/2$ still gives equal probabilities. These are deductions from the Born rule, not the paper's measured counts. Unequal detector efficiencies and finite samples must be accounted for when interpreting experimental frequencies; the ideal comparison assumes equally efficient detection.

A Rabi oscillation is a coherent oscillation of populations in a driven two-level quantum system.

A resonant drive can transfer a quantum state from one energy eigenstate to another and then bring its population back. This is unitary time evolution: the drive exchanges energy with the selected two-level system, while total probability stays one. One complete Rabi cycle has population period $2\pi /\Omega$; halfway through, a Rabi pi pulse transfers the whole population in the ideal model. This makes Rabi oscillations useful for calibrating resonant quantum control, preparing quantum superpositions and supplying the pulses of Ramsey interferometry.

Ramsey interferometry uses two coherent pulses with a free interval between them: the first creates a quantum superposition, the interval accumulates relative quantum phase, and the second converts that quantum phase into measurement probabilities. Its attraction is that an ordinary population readout can reveal the quantum coherence hidden from an energy measurement during the free interval.

Here is an ideal qubit model with an explicit pulse convention. Let $R_y(\theta )=\exp (-i\theta {\sigma }_y/2)$, prepare $|0⟩$, apply $R_y(\pi /2)$, evolve freely with $U_{\varphi }=\exp (-i\varphi {\sigma }_z/2)$, and apply the inverse pulse $R_y(-\pi /2)$ before a computational-basis measurement. The first pulse gives $|+⟩$. Direct matrix multiplication givesThus the Born rule gives $P(0)={\cos }^2(\varphi /2)$ and $P(1)={\sin }^2(\varphi /2)$. Before the second pulse, the same computational-basis measurement would give $1/2$ for both outcomes for every $\varphi$. At $\varphi =2\pi /3$, the second pulse changes these probabilities to $1/4$ and $3/4$. If the two pulses have the same sign instead, the output labels interchange: the pulse convention is part of the prediction.

In a pulse experiment the phase is measured relative to the oscillator supplying the coherent pulses. In a rotating frame, write the free Hamiltonian operator as $H_{\mathrm{rot}}=\hslash \delta {\sigma }_z/2$ with signed angular frequency difference $\delta$ in the chosen labeling convention. A free interval $T$ then accumulates $\varphi =\delta T$; changing the second pulse's phase changes the analysis direction. This comparison makes Ramsey interferometry useful for measuring frequency differences and diagnosing loss of quantum coherence.

For example, suppose the first pulse has prepared a density matrix with equal populations and transverse Bloch vector length $0\le V\le 1$:The same inverse analysis pulse yields $P(0)=(1+V\cos \varphi )/2$. The pure state case $V=1$ gives the full ideal fringe; the mixed state $I/2$ gives the flat result $1/2$. This is a controlled mathematical comparison, not a fit to an experiment. Experimental contrast alone does not identify quantum decoherence: imperfect preparation, pulse errors and readout errors can also reduce it.

The Schrödinger equation governs the evolution of a quantum state between instances of measurement in quantum mechanics. For a particle of mass $m>0$ moving in a prescribed real potential energy $V(\mathbf{r},t)$, its position wavefunction satisfiesHere $h$ is the Planck constant, ${\nabla }^2$ 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 $\rho =|\psi {|}^2$ and integration over a region.

This scalar Schrödinger operator describes nonrelativistic quantum mechanics, with relevant momenta much smaller than $mc$, where $c$ is the speed of light. It omits the constant rest energy $m{c}^2$, 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 $V$ 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 byInitial 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 $-itH(t)/\hslash$.

The local reason for conservation is especially tangible. Multiply the Schrödinger equation by the complex conjugate ${\psi }^{\ast }$, subtract its conjugate multiplied by $\psi$, and use the fact that the real potential energy terms cancel. For smooth wavefunctions this gives the probability continuity equationThe 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, $V=0$, initially in a Gaussian wave packet with position standard deviation ${\sigma }_0$ and mean momentum $\hslash k_0$: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 $\omega (k)=\hslash k^2/(2m)$, so each Fourier component gains the factor $e^{-i\hslash k^{2}t/(2m)}$. Completing the square in the resulting Gaussian integral giveswith the square root chosen continuously from $t=0$. Its position probability density and standard deviation areThe centre moves with the group velocity ${\omega }^{\prime }(k_0)=v_0$ 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 $\Delta p=\hslash /(2{\sigma }_0)$ and the initial position-momentum covariance vanishes. Free evolution gives $x(t)=x(0)+p(0)t/m$ in the Heisenberg picture, henceThis independently checks the spreading law. The Gaussian integral of $\rho$ is one at every time; the decreasing peak is exactly offset by increasing width.

As hypothetical preparation parameters, choose a free electron with ${\sigma }_0=1\mathrm{nm}$ and $k_0=2{\mathrm{nm}}^{-1}$. Using $\hslash \simeq 1.054571817\times 10^{-34}\mathrm{J}\mathrm{s}$ and $m\simeq 9.1093837\times 10^{-31}\mathrm{kg}$ givesThe rounded electron mass uses the 2022 adjustment.^[ref] At $t=t_s$, the centre is $4\mathrm{nm}$ and ${\sigma }_t=\sqrt{2}\mathrm{nm}$; at $t=2t_s$, they are $8\mathrm{nm}$ and $\sqrt{5}\mathrm{nm}$. The mean kinetic energy staysIt includes both mean motion and the momentum spread. The mean speed is about $7.72\times 10^{-4}c$ and the momentum spread corresponds to $1.93\times 10^{-4}c$, making the nonrelativistic model accurate for the appreciably occupied part of this ideal Gaussian distribution. Its mathematically unbounded tails are an idealization.

The complex width also explains where the probability goes. Write $\psi =\sqrt{\rho }{e}^{i\theta }$, where $\theta$ is the dimensionless quantum phase. Differentiating the exact solution givesRelative to the moving centre, the probability current $j_{\mathrm{rel}}=j-v_0\rho$ points outward on both sides for $t>0$. This is a flux statement, not an additional assumption about individual particle trajectories. At $t=t_s$ and $x-v_{0}t=\pm {\sigma }_0$, the relative flow speeds $j_{\mathrm{rel}}/\rho$ are $\pm {\sigma }_0/(2t_s)$, about $\pm 2.8942\times 10^4\mathrm{m}{\mathrm{s}}^{-1}$. The probability within a moving interval of half-width ${\sigma }_0$ falls from about $0.6827$ initially to $0.5205$ at $t_s$, exactly because of this outward probability current.

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.

Consider an original balanced-readout model, rather than a fit to the experiment. Let $n\in \{0,\dots ,N_{\max }\}$, with a normalized prior probability $p_0(n)\ge 0$ and ${\sum }_n{p}_0(n)=1$. A detected chemical atom has outcome $s\in \{e,g\}$. Define the conditional probabilitiesHere $\chi$ is the accumulated quantum phase per photon, $\varphi$ the controlled analysis quantum phase, and $C$ the fringe contrast. The two probabilities sum to one for every $n$. Actual calibration may also require an offset and asymmetric detection errors; $C$ alone does not specify a quantum instrument.

Assume no photon number change during the record and fresh independently prepared atomic meters. For observed outcomes $s_1,\dots ,s_k$ at recorded settings ${\varphi }_1,\dots ,{\varphi }_k$, the Born rule and Bayes' theorem giveThe posterior probability is defined only for $Z_k>0$. Summing its numerator over $n$ proves normalization. This is also the update of diagonal elements of a field density operator under a measurement in quantum mechanics; it does not assume that an initial quantum superposition was an unknown pre-existing integer. The observed settings are part of the record, especially when selected adaptively.

For an exact finite example take $p_0=(1,1,1,1)/4$, $C=4/5$, and $\chi =\pi /2$. An $e$ outcome at $\varphi =0$ has likelihood $(9,5,1,5)/10$, giving $p_1=(9,5,1,5)/20$. Another $e$ at $\varphi =-\pi /2$ gives $p_2=(9,9,1,1)/20$. A subsequent $g$ at $\varphi =0$ givesThus three atomic records raise the posterior probability of $n=1$ from $1/4$ to $45/68$, while leaving finite uncertainty. No one of these outcomes identifies $n$ with certainty.

The periodic likelihood function imposes an identifiability limit. If $\chi =2\pi /q$ for an integer $q>0$, then $L_s(n+q;\varphi )=L_s(n;\varphi )$ for every setting: the posterior odds between $n$ and $n+q$ never change. Varying $\varphi$ can resolve a cosine reflection ambiguity within a period, but cannot remove this exact period. Changing $\chi$ can supply additional information. A finite cutoff is therefore a physical assumption to justify, not evidence that excluded high photon numbers were absent. At $C=0$ no outcome changes the prior probability; at $0<C<1$ neither binary outcome has an exactly zero likelihood.

The no-loss assumption also has a time limit. In an original zero-temperature loss model, let $\kappa =1/T_c$. A Fock state $|n⟩$ loses its first photon at rate $n\kappa$, so its survival probability is $\exp (-n\kappa t)$ and its mean lifetime is $T_c/n$ for $n>0$. Several records must be acquired within that lifetime to identify a plateau. When losses matter, first predict $p_k^{-}(n)={\sum }_m{T}_{nm}(\Delta t)p_{k-1}(m)$ using a normalized transition matrix, then apply the atomic Bayesian inference update to $p_k^{-}$. Thermal excitation requires upward transitions too. A delayed or broadened inference after a jump need not indicate that an atom caused the jump.

A single-photon state lies in the one-excitation sector of the quantized electromagnetic field. It can be a quantum superposition over spatial ports or temporal modes while its total photon number remains one. Squared coefficients give exclusive event probabilities only when the corresponding wavepacket modes are orthonormal; overlapping temporal histories must be added coherently.

A Slater determinant extends two-particle antisymmetrization to $N$ fermions. Here $x_i$ contains position and spin or other internal variables. Distinct orthonormal one-particle states give the normalized wavefunction displayed in the title. Repeating a one-particle state repeats a column of the determinant, so the wavefunction vanishes. Exchanging particle slots exchanges rows and reverses its sign. A general fermionic quantum state can be a quantum superposition of Slater determinants; it need not be one Slater determinant.

Measure the commuting stabilizer generators $g_1=Z_1{Z}_2$ and $g_2=Z_2{Z}_3$ in that order. Use error syndrome bit $0$ for eigenvalue $+1$ and bit $1$ for eigenvalue $-1$. For $X_j$ the Pauli X gate on physical qubit $j$, every row holds for arbitrary normalized $\alpha ,\beta$. This ordered parity-check convention follows ^[ref]:

Each pair of computational basis words within a row has the same two parities. Consequently the projective measurement reveals the row without distinguishing its two amplitudes. For example, $|{\psi }_L⟩=(\sqrt{3}/2)|000⟩+(i/2)|111⟩$ becomes $(\sqrt{3}/2)|010⟩+(i/2)|101⟩$ under $X_2$. Both checks give $-1$ with probability one, and applying $X_2$ restores both coefficients exactly, including the relative phase $i$.

One implementation uses two quantum ancillas initialized to $|00⟩$. Two controlled-NOT gates from data qubits 1 and 2 to the first quantum ancilla write their parity; two more from data qubits 2 and 3 to the second quantum ancilla write the other parity. Measuring only the quantum ancillas in the computational basis gives $(s_1,s_2)$. Measuring the individual data qubits instead would reveal the encoded computational basis value and destroy the unknown quantum superposition. The existing coherent syndrome extraction for the three-qubit repetition code uses the alternative ordered checks $(Z_2{Z}_3,Z_1{Z}_3)$, so its bit labels differ from this table.

The recovery also handles coherent errors. Put $E_0=I,E_1=X_1,E_2=X_2,E_3=X_3$. Their images of the code are mutually orthogonal, giving $P{E}_a^{†}{E}_{b}P={\delta }_{ab}P$, the Knill--Laflamme condition for this family. For any linear combination $E={\sum }_{a=0}^3{c}_a{E}_a$, its component in sector $s(E_a)$ is $c_a{E}_a|{\psi }_L⟩$. Measuring that error syndrome has probability $|c_a{|}^2/{\sum }_b|c_b{|}^2$ on the normalized erroneous state, and recovery $E_a$ returns $|{\psi }_L⟩$. More generally one recovery quantum channel corrects every noise quantum channel whose Kraus operators belong to this span. The assumption concerns the whole jointly correctable family; it does not extend to the span of arbitrary single-site Pauli errors. ^[ref].

A fully computed coherent example is the unitary operator $U=\exp (-i\theta X_2)=\cos \theta I-i\sin \theta X_2$. At $\theta =\pi /6$ it givesThe error syndromes $00$ and $11$ occur with probabilities $3/4$ and $1/4$, independently of $\alpha ,\beta$. Doing nothing for $00$ and applying $X_2$ for $11$ restores the same logical qubit in both branches; the scalar $-i$ in the second branch is an irrelevant global phase. This is an example of linearity of coherent quantum error correction, rather than an assumption that the original error was a classical random bit flip.

For a time-independent Hamiltonian operator $H$, an energy eigenstate satisfiesFor a one-dimensional particle of mass $m$ in a real potential energy $U(x)$,Its full time-dependent wavefunction is $\Psi (x,t)=\psi (x)e^{-iEt/\hslash }$: substitution in the Schrödinger equation checks both the minus sign and the factor $\hslash$. 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 $\Psi (x,t)={\sum }_n{c}_n{\psi }_n(x)e^{-i{E}_{n}t/\hslash }$, with ${\sum }_n|c_n{|}^2=1$. Cross terms in $|\Psi {|}^2$ oscillate at angular frequencies $(E_n-E_m)/\hslash$. 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 $C$ to $U$ shifts every energy eigenvalue by $C$ and multiplies the full state by $e^{-iCt/\hslash }$, a global phase; energy differences, position probabilities and probability current are unchanged.

Trapped-ion quantum control prepares an ion's internal quantum state, couples it to quantized motion with laser pulses, and measures the result through state-dependent fluorescence. Its appeal is access to both a controllable two-level system and a quantum harmonic oscillator in one apparatus.

Confinement, preparation, control and readout solve different problems. The Paul trap's RF electric field confines the ion through its RF-averaged trapping potential; it still produces Paul trap micromotion. Doppler cooling and resolved sideband cooling reduce motional excitation. Optical pumping selects an internal quantum state. A coherent pulse then changes amplitudes, while electron shelving converts an internal population into a detectable stream of photons. A bright/dark count by itself does not reveal a quantum superposition's relative phase; an additional controlled rotation before detection can make that phase affect the measured population.

Consider one chosen model at a fixed total energy, with two open scattering channels, no absorption and no omitted outgoing modes. Channel 1 leaves the target unchanged; channel 2 excites it. Let $v_1=v>0$ and $v_2=v/2$. For the same projectile mass, conservation of energy then requires a target excitation $\Delta =3m{v}^2/8$. These are hypothetical parameters, not a fit to a particular atom. Each incoming or outgoing mode is normalized to unit probability current; the outgoing modes will have the angular profiles specified below.

Choose the S-matrixThis phenomenological unitary matrix specifies the channel conversion; it is not derived here from an interaction potential energy. For incidence in channel 1, its first column gives the outgoing flux amplitudesTo recover ordinary wave coefficients relative to an incident coefficient one, undo the velocity normalization:Thus $|c_2{|}^2=32/25>1$ is allowed: channel 2 travels more slowly. Its outgoing probability current fraction is $16/25$, not $32/25$. Renormalizing the unweighted squares would give $9/41$ and $32/41$, the wrong probabilities. More generally the ordinary coefficient map is $A=V^{-1/2}S{V}^{1/2}$, with $V=\mathrm{diag}(v,v/2)$; it obeys $A^{†}VA=V$. The conserved quantity is the squared coefficient norm weighted by velocity, which represents probability current.

Now specify an original angular realization of these two outgoing modes, with scattering angle $\theta$ and no azimuth dependence:The first normalization follows from $2\pi {\int }_{-1}^1(1+u^2)du=16\pi /3$, where $u=\cos \theta$. These angular profiles define the modes of this model; they are not asserted to follow from Born's approximation or an atomic potential energy. The joint directional probability densities are

Suppose a detector accepts the cone $D:0\le \theta \le \pi /3$, with full azimuth and unit efficiency. Its angular fractions within the channels areMultiplying the angular fractions by the outgoing flux weights givesFor $1600$ independent incident events, the expected values are therefore $171$ and $256$ accepted detections, compared with $576$ and $1024$ over all directions. The remaining expected $1173$ events fall outside the cone. The chosen incident number makes the expectations integers; actual counts fluctuate. Conditional on a detection inside the cone, the channel probabilities are $171/427$ and $256/427$, rather than $9/25$ and $16/25$. Angular acceptance changes the observed channel proportions even with equal detector efficiencies.

The final atomic states are orthogonal, so ignoring the atomic readout adds the two channel intensities; it does not produce $|b_1{g}_1+b_2{g}_2{|}^2$. Coherent paths that reach the same final state can instead interfere, as the existing probability amplitude and quantum superposition examples demonstrate. The distinction is physical distinguishability of final states, not whether the detector operator chooses to display their labels. Flux normalization and angular integration make the scattering interpretation useful: they turn complex outgoing waves into conserved event weights and then into the statistics of a specified detector.

For a deliberately simplified neutrino oscillation calculation, take two neutrino flavors related to two mass eigenstates by a real rotation matrix with angle $\theta$. A quantum superposition acquires different phases during propagation; projecting back onto the initial neutrino flavor gives the probabilityHere $L$ is path length and $E$ is neutrino energy; the numerical coefficient incorporates unit conversion.

The paper states the corresponding conversion formula.^[ref]

The path dependence has a geometrical origin. In an original spherical-Earth model, let $R$ be the radius and $h$ a fixed production height, neglecting detector depth. Let $\Theta$ be the zenith distance of the direction in the sky from which the neutrino arrives, so $c=\cos \Theta =+1$ means downward propagation. The backward ray from detector to production obeys $R^2+L^2+2RLc=(R+h{)}^2$. Its positive solution isThus $L(+1)=h$ and $L(-1)=2R+h$. The calculation explains why opposite directions test very different baselines, rather than attributing every observed direction to the same travel distance. Actual production heights and the angle between the outgoing charged lepton and incoming neutrino smear that relationship.

For an original chosen example, set $\theta =\pi /4$ and $\Delta m^2=2.5\times 10^{-3}{\mathrm{eV}}^2$. At $E=1\mathrm{GeV}$, $L=20\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.9960$, while $L=500\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.000279$. The first minimum is at $L/E\simeq 494.7km/GeV$. These are stipulated inputs, not a fit to the 1998 data. A broad distribution of true $L/E$, together with energy and directional uncertainty, smooths the sharp minima.

A direct amplitude check makes the mechanism explicit. DefineThen $S_{11}={\cos }^2\theta +{\sin }^2\theta e^{-2i\phi }$. Taking its squared absolute value yields the displayed survival formula because $1-\cos (2\phi )=2{\sin }^2\phi$. Unitarity gives $|S_{11}{|}^2+|S_{21}{|}^2=1$: disappearance from one channel becomes conversion to the other within this closed model.

This illustration does not describe all solar neutrino propagation. Coherent forward scattering in matter changes the effective propagation operator, and three active neutrino flavors require a larger mixing description. Even in the simplified model, neutrino oscillation measures a squared-mass difference; it does not determine the absolute neutrino mass.

A wavefunction is the position representation of a pure state. In the one-dimensional model used here, $\psi (x,t)$ takes values in the complex numbers, and the Born rule givesThe probability density $\rho$ has units of inverse length; the wavefunction has units of inverse square root of length. A value of $\rho$ at a point is not itself a probability: probabilities come from an integral over a region. For the infinite square well ground state, $\rho (L/2)=2/L$, whereas the probability of the exact single point $L/2$ is zero. The interesting question is how much probability a finite interval contains.

Two wavefunctions differing only by a global phase represent the same pure state. Relative quantum phases between components of a quantum superposition affect matter-wave interference, so retaining the full complex number amplitude is essential. This scalar description covers a particle without active spin degrees of freedom; a general electron wavefunction can have two spin components, with the position probability density obtained by adding the squared modulus of each component.

Wavefunction normalization makes the total Born rule probability equal to one. In an orthonormal basis of energy eigenstates, $\psi ={\sum }_n{c}_n{\varphi }_n$ is normalized precisely when ${\sum }_n|c_n{|}^2=1$.

For a proposed infinite square well mode $A\sin (n\pi x/L)$, the integral of its squared modulus is $|A{|}^{2}L/2$. Thus wavefunction normalization requires $|A|=\sqrt{2/L}$; choosing $A$ real and positive fixes a convenient global phase. For the equal quantum superposition $({\varphi }_1+{\varphi }_2)/\sqrt{2}$, orthogonality makes the integrated cross term zero, leaving $1/2+1/2=1$. The cross term can nevertheless be nonzero locally: normalization constrains total probability, not the shape of its probability density.

For $M$ complete energy eigenstates of a one-particle Quantum Hamiltonian for identical noninteracting fermions, order their energies ${ϵ}_1\le \cdots \le {ϵ}_M$. At fixed particle number $0\le N\le M$, the Quantum Hamiltonian and the particle-number constraint areThe ground state occupies the $N$ lowest complete one-particle states. A short minimization argument suffices: moving an occupant from a higher energy to an empty lower energy decreases the total energy. This filling describes the ground state at temperature $T=0$; exclusion by itself also applies to excited quantum states and at positive temperature. The underlying thermal description is the Fermi-Dirac distribution.^[ref]

If the last occupied energy belongs to a group of $g$ one-particle states at the same energy and $q$ of them must be occupied, with all lower states filled, the ground state subspace has dimension $\left(\genfrac{}{}{0ex}{}{g}{q}\right)$. Exclusion does not select a particular basis or a particular quantum superposition in that subspace. The factor of two for an electron spatial orbital assumes both spin projections are included; a model restricted to one spin projection has only one state per spatial orbital. A magnetic field can split the degenerate energy levels, and interactions can change the lowest-energy electron configuration. Thus this finite noninteracting filling rule is not a complete prescription for interacting chemical atoms.