The Feynman Path Integral Formulation

Wave mechanics evolves a state with a differential operator; matrix mechanics turns operators on a Hilbert space. Feynman's path integral does something that feels almost illegal: it computes the amplitude to go from one place to another by adding up a complex number for every conceivable trajectory between them, classical or not, smooth or jagged. The number attached to each path is a pure phase, eiS[x]/ℏe^{iS[x]/\hbar}, where SS is the same classical action that Lagrangian mechanics extremizes. Nothing is quantized by hand. Nothing collapses. The whole of quantum mechanics is repackaged as interference between histories. Why is this worth a chapter of its own? Because it makes the classical limit obvious rather than mysterious: when S/ℏS/\hbar is enormous the phases of neighbouring paths whirl out of step and cancel, except right around the trajectory where SS is stationary — and that trajectory is Newton's. The principle of least action stops being an axiom and becomes what survives quantum interference as ℏ→0\hbar\to0. It is also the formulation that scales: replace "sum over particle paths" with "sum over field configurations" and you have the engine of quantum field theory, lattice gauge theory, and — in imaginary time — statistical mechanics. You are ready for this. You command the time-evolution operator, the quantum propagator K=⟨xf∣U^∣xi⟩K=\langle x_f|\hat U|x_i\rangle and its composition law, Gaussian integrals, and the classical action. Here we assemble them: we build KK as a genuine integral over paths by slicing time into NN steps, prove the construction is consistent with the propagator you already know, and watch Newton's law fall out of stationary phase. The interpretive knots of the relativistic wave equations — negative energies, indefinite densities — untie cleanly in this language, which is exactly why the chapter ends here.

Key takeaways
  • The amplitude to travel from (xi,ti)(x_i,t_i) to (xf,tf)(x_f,t_f) is K=∫Dx eiS[x]/ℏK=\displaystyle\int\mathcal D x\,e^{iS[x]/\hbar}: a coherent sum over every path, each weighted by the phase of its classical action.
  • The functional integral is defined by time-slicing: cut the evolution into NN steps, insert position completeness, and take an (N−1)(N{-}1)-fold Gaussian integral with per-slice measure (m/2πiℏ ϵ)1/2(m/2\pi i\hbar\,\epsilon)^{1/2}.
  • Slicing is consistent with the propagator you already know: two free short-time kernels convolve into one long-time kernel, with the elapsed times simply adding.
  • The classical path is not put in by hand — it is the path of stationary phase. As ℏ→0\hbar\to0 only δS=0\delta S=0 survives, and δS=0\delta S=0 is mx¨=−V′(x)m\ddot x=-V'(x).
  • Contributing paths are nowhere-differentiable: "velocity" along a path is undefined, which is the uncertainty principle wearing a different coat.
  • KEY IDEA: quantum mechanics is interference between histories — the amplitude is ∑pathseiS/ℏ\sum_{\text{paths}}e^{iS/\hbar}, and the classical world is just the stationary-phase path that does not cancel.
The core insight

The quantum amplitude to go from one spacetime point to another is the sum over all paths of eiS[x]/ℏe^{iS[x]/\hbar}; the classical trajectory is simply the path of stationary phase that survives interference as ℏ→0\hbar\to0.

Where this lives today

The path integral is the working language of modern physics wherever many degrees of freedom interfere. In lattice QCD the sum over gluon and quark field configurations is done by Monte Carlo sampling of the Euclidean weight — Wiener's real, positive measure — on supercomputers, yielding hadron masses from first principles. In condensed matter the same imaginary-time integral computes finite-temperature partition functions. And in quantum computing, the amplitude of a circuit is literally a sum over computational histories, so path-integral intuition — constructive interference on the answer, destructive elsewhere — is how algorithm designers reason about speedups.

The physical picture

Picture the particle as an explorer who tries every route between AA and BB at once, each route carrying a little clock whose hand turns at a rate set by the Lagrangian. The amplitude is the vector sum of all the clock hands at the finish. Where many neighbouring routes hand in nearly-aligned arrows — near the path where the action is stationary — the arrows stack up and reinforce. Where the action changes quickly from one route to its neighbour, the arrows fan out in all directions and cancel to almost nothing. The knob that decides how ruthless this cancellation is happens to be ℏ\hbar: it is the unit of action against which SS is measured. For a thrown baseball S/ℏ∼1034S/\hbar\sim10^{34} (the verifier pins this at ≈5.5×1034\approx 5.5\times10^{34} for a 40 m/s40\,\text{m/s} toss), so shifting the path by a hair scrambles the phase completely and exactly one trajectory survives — the classical one. For an electron in an atom S/ℏ∼1S/\hbar\sim1, so a whole tube of paths contributes and the quantum spread is not a correction but the physics.

Figure
Sum over histories: arrows add near the stationary path, cancel elsewhere
A B stationary path δS = 0 near δS=0: arrows aligned → add far from δS=0: arrows scatter → cancel amplitude = Σ over paths of e^(iS/ℏ)

Because the surviving paths are jagged — a typical velocity diverges as the time slice shrinks — the path integral carries the uncertainty principle inside its measure: you cannot pin position and a smooth velocity together, because "velocity" along a contributing history is not even defined. The next figure is the punchline of the whole formulation: the same set of paths, evaluated at large versus small ℏ\hbar, gives a classical world in one limit and a smeared quantum one in the other, with no extra ingredient.

Figure
Same paths, two values of ℏ: the classical limit is stationary phase
weight cos(S/ℏ) vs path deformation q small S/ℏ (large ℏ) broad peak → many paths matter large S/ℏ (small ℏ) sharp spike at q=0 (δS=0) → only the classical path

The core equations

We pin the convention once and keep it for the whole chapter. A path integral means ∫Dx eiS/ℏ\int\mathcal D x\,e^{iS/\hbar} with the nonrelativistic action S=∫(12mx˙2−V) dtS=\int(\tfrac12 m\dot x^2 - V)\,dt, and every free short-time kernel carries the Fresnel prefactor (m/2πiℏt)1/2(m/2\pi i\hbar t)^{1/2}. The first equation is the definition itself; before it we need only the idea that the propagator you met in The Quantum Propagator — the amplitude K=⟨xf∣U^∣xi⟩K=\langle x_f|\hat U|x_i\rangle that already obeys composition and reduces to a delta at coincident times — can be written as a genuine integral over trajectories. The compact continuum statement of Feynman's formulation collapses three semesters of formalism into one line.

K(xf,tf;xi,ti)=∫x(ti)=xix(tf)=xf ⁣ ⁣Dx(t)  e iS[x]/ℏ(1)K(x_f,t_f;x_i,t_i)=\int_{x(t_i)=x_i}^{x(t_f)=x_f}\!\!\mathcal D x(t)\;e^{\,iS[x]/\hbar} \qquad (1)

where: KK — the propagator, the amplitude ⟨xf∣e−iH^T/ℏ∣xi⟩\langle x_f|e^{-i\hat HT/\hbar}|x_i\rangle; Dx(t)\mathcal D x(t) — the functional measure, summing over all paths with the endpoints held fixed; S[x]S[x] — the classical action of each path; ℏ\hbar — the natural unit of action. Read physically: (1) says a quantum amplitude is a democratic vote of histories, and the only thing a path contributes is its phase — every path has weight of unit modulus, so what distinguishes them is purely how their phases interfere. The phase is fixed by the same action that classical mechanics extremizes, which is what makes the classical limit inevitable rather than imposed.

S[x]=∫titf ⁣[12mx˙2−V(x)] dt(2)S[x]=\int_{t_i}^{t_f}\!\Big[\tfrac12 m\dot x^2 - V(x)\Big]\,dt \qquad (2)

where: mm — mass; x˙\dot x — the velocity along a given path (not an operator here, just a slope); V(x)V(x) — potential; L=T−VL=T-V — the Lagrangian. The physical reading: because LL appears in the exponent divided by ℏ\hbar, action differences of order ℏ\hbar produce phase differences of order one — that is the scale at which "which path?" becomes a physical question. Equation (1) is a promissory note until we say what Dx\mathcal D x is. Time-slicing cashes it: chop T=tf−tiT=t_f-t_i into NN pieces of width ϵ=T/N\epsilon=T/N, insert a complete set of position states between every factor of e−iH^ϵ/ℏe^{-i\hat H\epsilon/\hbar}, and you get a concrete multidimensional integral.

K=lim⁡N→∞(m2πiℏ ϵ) ⁣N/2 ⁣∫ ⁣∏j=1N−1 ⁣dxj  exp⁡ ⁣[iℏ∑j=1N(m(xj−xj−1)22ϵ−ϵV(xˉj))](3)K=\lim_{N\to\infty}\Big(\frac{m}{2\pi i\hbar\,\epsilon}\Big)^{\!N/2}\!\int\!\prod_{j=1}^{N-1}\!dx_j\;\exp\!\Big[\frac{i}{\hbar}\sum_{j=1}^{N}\Big(\frac{m(x_j-x_{j-1})^2}{2\epsilon}-\epsilon V(\bar x_j)\Big)\Big] \qquad (3)

where: NN — number of slices; ϵ=T/N\epsilon=T/N; xjx_j — the intermediate position at time ti+jϵt_i+j\epsilon, with x0=xi, xN=xfx_0=x_i,\,x_N=x_f; xˉj\bar x_j — the midpoint where VV is sampled. Read it: the exponent is a discrete Riemann sum for ∫L dt\int L\,dt, and the crucial detail is the per-slice normalization (m/2πiℏϵ)1/2(m/2\pi i\hbar\epsilon)^{1/2} — get that wrong and the continuum limit diverges or vanishes.

Derivation AThe functional measure from time-slicing

We want to turn the formal continuum object ∫Dx eiS/ℏ\int\mathcal D x\,e^{iS/\hbar} into an honest product of ordinary Gaussian integrals, and in doing so derive the per-slice measure rather than postulating it. The strategy is standard but the payoff is precise: factor the evolution operator, insert position completeness at each junction, and evaluate one short-time matrix element exactly by passing through momentum space. Watch the free-particle prefactor (m/2πiℏϵ)1/2(m/2\pi i\hbar\epsilon)^{1/2} appear as the Jacobian of a single Gaussian pp-integral — this is the piece the measure must carry.

Step 1. Split the evolution into NN identical short steps, e−iH^T/ℏ=(e−iH^ϵ/ℏ)Ne^{-i\hat HT/\hbar}=\big(e^{-i\hat H\epsilon/\hbar}\big)^N with ϵ=T/N\epsilon=T/N. Why: each factor acts over a time short enough that the noncommutativity of T^\hat T and V^\hat V costs only O(ϵ2)O(\epsilon^2) (the Trotter splitting), which vanishes in the N→∞N\to\infty limit.

Step 2. Insert ∫∣x⟩⟨x∣ dx=1\int|x\rangle\langle x|\,dx=\mathbb 1 between every pair of factors, giving K=∫∏j=1N−1dxj  ∏j=1N⟨xj∣e−iH^ϵ/ℏ∣xj−1⟩K=\int\prod_{j=1}^{N-1}dx_j\;\prod_{j=1}^{N}\langle x_j|e^{-i\hat H\epsilon/\hbar}|x_{j-1}\rangle. Why: this trades an operator on Hilbert space for a chain of cc-number matrix elements linked by integrations over the intermediate positions xjx_j — exactly the "stopovers" the propagator's composition law anticipates.

Step 3. Evaluate one short-time element by going to momentum space: ⟨xj∣e−iH^ϵ/ℏ∣xj−1⟩=∫dp2πℏ eiℏ[p Δx−ϵH(p,xˉj)]\langle x_j|e^{-i\hat H\epsilon/\hbar}|x_{j-1}\rangle=\int\frac{dp}{2\pi\hbar}\,e^{\frac{i}{\hbar}[p\,\Delta x-\epsilon H(p,\bar x_j)]} with Δx=xj−xj−1\Delta x=x_j-x_{j-1}. Why: with H=p2/2m+VH=p^2/2m+V, the pp-integral is a pure Gaussian; completing the square in pp and doing the Fresnel integral produces the prefactor and the discretized kinetic term.

Step 4. The Gaussian pp-integral gives ⟨xj∣e−iH^ϵ/ℏ∣xj−1⟩=(m2πiℏϵ)1/2exp⁡ ⁣[iℏ(mΔx22ϵ−ϵV(xˉj))]\langle x_j|e^{-i\hat H\epsilon/\hbar}|x_{j-1}\rangle=\big(\tfrac{m}{2\pi i\hbar\epsilon}\big)^{1/2}\exp\!\big[\tfrac{i}{\hbar}\big(\tfrac{m\Delta x^2}{2\epsilon}-\epsilon V(\bar x_j)\big)\big]. Why: the coefficient (m/2πiℏϵ)1/2(m/2\pi i\hbar\epsilon)^{1/2} is precisely the Jacobian of one Fresnel integral — the origin of the functional measure's N/2N/2-power normalization in (3). The verifier reproduces this exactly by completing the square symbolically.

Step 5. Multiply the NN short-time elements. The exponents add to iℏ∑jϵ(12m(Δx/ϵ)2−V)→iℏ∫L dt=iℏS[x]\frac{i}{\hbar}\sum_j\epsilon\big(\tfrac12 m(\Delta x/\epsilon)^2-V\big)\to\frac{i}{\hbar}\int L\,dt=\frac{i}{\hbar}S[x] as ϵ→0\epsilon\to0, and the product of prefactors defines Dx\mathcal D x. Why: the Riemann sum in the exponent is the action; everything not in the exponent is bookkeeping that we bundle into the symbol Dx\mathcal D x — which is why (1) is only meaningful together with the slicing rule (3). *Verified: the short-time free kernel =m/2πiℏϵ exp⁡(imΔx2/2ℏϵ)=\sqrt{m/2\pi i\hbar\epsilon}\,\exp(im\Delta x^2/2\hbar\epsilon) is confirmed symbolically by the Fresnel pp-integral, and ∣prefactor∣=m/2πℏϵ|{\rm prefactor}|=\sqrt{m/2\pi\hbar\epsilon}.*

The definition would be worthless if it disagreed with the propagator you already trust. It does not, and the cleanest proof is a consistency check: two free short-time kernels must convolve into one longer-time kernel, with the elapsed times simply adding. This is the sum-over-paths saying nothing more than "a path from ii to ff is a path from ii to a midpoint followed by a path from the midpoint to ff."

Go deeper · Derivation B · 1 minConsistency with the propagator: times add under convolution
Derivation BConsistency with the propagator: times add under convolution

The claim to establish is that the slicing in (3) is not an arbitrary discretization but is forced to agree with the known propagator. We test the smallest nontrivial case: glue two free kernels of duration ϵ\epsilon and show they produce the single free kernel of duration 2ϵ2\epsilon. The route is a Gaussian convolution over the shared midpoint xmx_m; the tell that everything is consistent is that the two durations add in the final kernel and the prefactor tracks t=2ϵt=2\epsilon automatically. This is the microscopic version of the composition law that The Quantum Propagator states as an axiom.

Step 1. Write the free kernel K0(xf,xi;t)=m/2πiℏt  eim(xf−xi)2/2ℏtK_0(x_f,x_i;t)=\sqrt{m/2\pi i\hbar t}\;e^{im(x_f-x_i)^2/2\hbar t} and form the convolution ∫dxm K0(xf,xm;ϵ)K0(xm,xi;ϵ)\int dx_m\,K_0(x_f,x_m;\epsilon)K_0(x_m,x_i;\epsilon). Why: concatenating a step i→mi\to m with a step m→fm\to f, then summing over the free intermediate xmx_m, is exactly the two-slice case of (3) with V=0V=0.

Step 2. Collect the exponent: (xf−xm)2+(xm−xi)2=2 xm2−2(xf+xi)xm+xf2+xi2(x_f-x_m)^2+(x_m-x_i)^2=2\,x_m^2-2(x_f+x_i)x_m+x_f^2+x_i^2, a quadratic in xmx_m with minimum at xm⋆=(xf+xi)/2x_m^\star=(x_f+x_i)/2. Why: the midpoint that dominates is the average of the endpoints — the straight-line classical path — and the leftover after completing the square is the residual action.

Step 3. The residual at xm⋆x_m^\star is (xf−xi)2/2(x_f-x_i)^2/2, so the surviving exponent is im2ℏϵ⋅(xf−xi)22=im(xf−xi)22ℏ(2ϵ)\dfrac{im}{2\hbar\epsilon}\cdot\dfrac{(x_f-x_i)^2}{2}=\dfrac{im(x_f-x_i)^2}{2\hbar(2\epsilon)}. Why: the denominator is now 2ϵ2\epsilon — the two durations have added. The verifier confirms the residual is exactly (xf−xi)2/2(x_f-x_i)^2/2 and the exponent matches the single kernel at t=2ϵt=2\epsilon.

Step 4. The Gaussian xmx_m-integral supplies exactly the prefactor conversion (m/2πiℏϵ)2⋅2πiℏϵ/2m=m/2πiℏ(2ϵ)\big(\sqrt{m/2\pi i\hbar\epsilon}\big)^2\cdot\sqrt{2\pi i\hbar\epsilon/2m}=\sqrt{m/2\pi i\hbar(2\epsilon)}, so the whole convolution equals K0(xf,xi;2ϵ)K_0(x_f,x_i;2\epsilon). Why: prefactor and phase both land on the t=2ϵt=2\epsilon kernel with no leftover constant — the discretization is self-consistent, which is what licenses taking N→∞N\to\infty. *Verified: the symbolic exponent match and a numerical Fresnel convolution both reproduce K0(xf,xi;2ϵ)K_0(x_f,x_i;2\epsilon) to <5×10−3<5\times10^{-3}.*

Finally, the payoff that makes the whole reformulation feel inevitable: classical mechanics is not assumed, it is extracted by stationary phase. When S/ℏS/\hbar is large the integrand oscillates violently except where SS is stationary, and there the vote is unanimous.

δSδx(t)=0  ⟹  mx¨=−∂V∂x(4)\frac{\delta S}{\delta x(t)}=0\;\Longrightarrow\;m\ddot x=-\frac{\partial V}{\partial x} \qquad (4)

where: δS/δx\delta S/\delta x — the functional derivative of the action; the arrow is the Euler–Lagrange equation for L=12mx˙2−VL=\tfrac12 m\dot x^2-V. Read physically: (4) is Newton's second law, recovered as the only path that does not phase-cancel as ℏ→0\hbar\to0. The dedicated treatment of how the leading correction organises itself — the stationary-phase expansion that turns this into semiclassics — is taken up by the next article, Stationary Phase and the Classical Limit; here we only need that δS=0\delta S=0 selects the classical trajectory.

Limiting cases

caseresultmeaning
ℏ→0\hbar\to0only δS=0\delta S=0 survivesclassical mechanics = stationary phase; least action is emergent
free particle, V=0V=0K0=m/2πiℏt eim(xf−xi)2/2ℏtK_0=\sqrt{m/2\pi i\hbar t}\,e^{im(x_f-x_i)^2/2\hbar t}single Gaussian; convolutions add durations
N=1N=1 (one slice)K≈K0(ϵ)K\approx K_0(\epsilon)exact for the free particle; O(ϵ2)O(\epsilon^2) error with VV
macroscopic S/ℏ≫1S/\hbar\gg1classical path onlybaseball: S/ℏ∼5×1034S/\hbar\sim5\times10^{34}
S/ℏ∼1S/\hbar\sim1many paths contributeelectron in an atom; quantum spread is physical
imaginary time t→−iτt\to-i\taueiS/ℏ→e−SE/ℏe^{iS/\hbar}\to e^{-S_E/\hbar}oscillatory → convergent; the Euclidean chapter ahead

Worked examples

Example 1 — Two slices really do add their times. Take a free particle (V=0V=0, m=1m=1, ℏ=1\hbar=1) and check the two-slice construction against the known kernel, choosing ϵ=0.37\epsilon=0.37, xi=−0.3x_i=-0.3, xf=0.8x_f=0.8.

Step 1. The single-step kernel is K0(xf,xi;ϵ)=1/(2πi⋅0.37) ei(1.1)2/(2⋅0.37)K_0(x_f,x_i;\epsilon)=\sqrt{1/(2\pi i\cdot0.37)}\,e^{i(1.1)^2/(2\cdot0.37)}, with (xf−xi)2=1.21(x_f-x_i)^2=1.21.

Step 2. Convolve two such kernels over the midpoint: the exponent collapses to i (1.1)22(2⋅0.37)=1.21 i1.48\dfrac{i\,(1.1)^2}{2(2\cdot0.37)}=\dfrac{1.21\,i}{1.48}, i.e. the same (xf−xi)2(x_f-x_i)^2 but with duration 2ϵ=0.742\epsilon=0.74.

Step 3. The prefactor lands on 1/(2πi⋅0.74)\sqrt{1/(2\pi i\cdot0.74)} automatically. Evaluating the convolution numerically gives the exact K0(xf,xi;0.74)K_0(x_f,x_i;0.74) to within 5×10−35\times10^{-3}.

Answer: the two durations add to 0.740.74 and the result is the single free propagator over 2ϵ2\epsilon — the slicing is consistent with the propagator you already trust.

Example 2 — Why the baseball has no quantum spread. Estimate S/ℏS/\hbar for a thrown baseball and show it kills all but the classical path.

Step 1. Take m=0.145 kgm=0.145\,\text{kg}, v=40 m/sv=40\,\text{m/s}, flight time t=0.5 st=0.5\,\text{s}. The action scale is S∼12mv2 t=12(0.145)(402)(0.5)=58 J ⋅ sS\sim\tfrac12 m v^2\,t=\tfrac12(0.145)(40^2)(0.5)=58\,\text{J}\,\cdot\,\text{s}.

Step 2. Divide by ℏ=1.0546×10−34 J ⋅ s\hbar=1.0546\times10^{-34}\,\text{J}\,\cdot\,\text{s}: S/ℏ=58/1.0546×10−34≈5.5×1034S/\hbar=58/1.0546\times10^{-34}\approx5.5\times10^{34}.

Step 3. A path deformation changing SS by even one part in 103410^{34} shifts the phase by ∼1\sim1 radian and cancels. Only deformations with δS≲ℏ\delta S\lesssim\hbar survive, i.e. essentially none.

Answer: S/ℏ≈5.5×1034S/\hbar\approx\mathbf{5.5\times10^{34}}, so the interference window around the classical path is astronomically narrow — the baseball has a single, sharp trajectory, exactly as Newton says.

Origins · 1 minNorbert Wiener and how this was discovered
NW
Norbert Wiener
1894–1964 · MIT, Cambridge MA

Long before Feynman, Wiener gave a rigorous meaning to "integrating over paths" — but in imaginary time. His 1923 Wiener measure assigns a genuine probability to sets of continuous, nowhere-differentiable Brownian trajectories, and it is a mathematically honest measure precisely because the weight e−SE/ℏe^{-S_E/\hbar} is real and positive. Feynman's oscillatory eiS/ℏe^{iS/\hbar} is not a measure in Wiener's sense — its "cancellation, not smallness" is what makes the real-time path integral notoriously hard to define — but Wick-rotating back to imaginary time lands you squarely in Wiener's world. That is why the Euclidean path integral, not the Minkowski one, is where mathematicians and lattice physicists actually compute.

History
1923
— Wiener defines the Wiener measure on Brownian paths (rigorous imaginary-time path integration)
1933
— Dirac notes that eiS/ℏe^{iS/\hbar} "corresponds to" the short-time amplitude
1942
— Feynman's Princeton PhD thesis builds least action into quantum mechanics
1948
— Feynman publishes the space-time (path-integral) formulation of non-relativistic QM
1951
— Morette gives the semiclassical (Van Vleck–Pauli–Morette) prefactor a firm footing

Going further

  • Prerequisites: the time-dependent Schrödinger equation; The Quantum Propagator; Gaussian integrals; Lagrangian and Hamiltonian mechanics.
  • Leads to: The Propagator from Path Integrals; Free-Particle and Harmonic-Oscillator Path Integrals; Imaginary-Time and Euclidean Path Integrals; and Stationary Phase and the Classical Limit.
  • Provenance: upgraded from the AllQFT `qm` base (QUANT.RQM.PATHI.1092), Griffiths §-independent (Feynman & Hibbs, Quantum Mechanics and Path Integrals). The historical Feynman band is owned by Time-Dependent Schrödinger Equation; the propagator's composition law, spectral form, and free/oscillator kernels are owned by The Quantum Propagator — recalled here by name, not re-derived.
Is the path integral really equivalent to the Schrödinger equation, or just an approximation?

Exactly equivalent. The article Path Integral and the Schrödinger Equation proves the reverse direction: evolving a wavefunction by one short-time kernel and Taylor-expanding reproduces iℏ∂tψ=H^ψi\hbar\partial_t\psi=\hat H\psi term by term at O(ϵ)O(\epsilon). The two formulations are the same physics written with different primitives — a differential operator versus a sum over histories. Neither is more fundamental; you pick whichever makes a given problem easy.

If every path has the same modulus, in what sense is the classical path "preferred"?

Not in modulus — in phase coherence. Each path contributes a unit-modulus arrow eiS/ℏe^{iS/\hbar}; the classical path is preferred because its immediate neighbours have almost the same action, so their arrows point the same way and add. Away from stationarity the action changes quickly, arrows fan out, and the sum cancels. So "preferred" means "surrounded by constructively-interfering neighbours," which is a statement about the second variation of SS, not about any single path being bigger.

Why the factor of ii in the measure (m/2πiℏϵ)1/2(m/2\pi i\hbar\epsilon)^{1/2}?

It comes from the Gaussian momentum integral in each short-time step: ∫dp e−iϵp2/2mℏ\int dp\,e^{-i\epsilon p^2/2m\hbar} is a Fresnel integral, and completing the square produces 2πmℏ/iϵ\sqrt{2\pi m\hbar/i\epsilon}-type factors. The ii is the fingerprint of real-time (oscillatory) evolution. Wick-rotate to imaginary time and the ii disappears, leaving the real diffusion prefactor m/2πℏτ\sqrt{m/2\pi\hbar\tau} — the sign that you have entered Wiener's convergent world.

The dominant paths have infinite velocity — doesn't that break the kinetic term?

It looks alarming but it is consistent. On a slice, ⟨Δx2⟩∼ℏϵ/m\langle\Delta x^2\rangle\sim\hbar\epsilon/m, so Δx∼ϵ\Delta x\sim\sqrt{\epsilon} and the "velocity" Δx/ϵ∼ϵ−1/2\Delta x/\epsilon\sim\epsilon^{-1/2} diverges as ϵ→0\epsilon\to0. The kinetic term mΔx2/2ϵ∼ℏ/2m\Delta x^2/2\epsilon\sim\hbar/2 stays finite, though, because it is Δx2/ϵ\Delta x^2/\epsilon (not (Δx/ϵ)2ϵ(\Delta x/\epsilon)^2\epsilon) that appears. The measure is built precisely so this works: it is the same Δx∼ϵ\Delta x\sim\sqrt{\epsilon} scaling that encodes the uncertainty principle.

Does the midpoint prescription xˉj=(xj+xj−1)/2\bar x_j=(x_j+x_{j-1})/2 for VV matter?

For a smooth potential the choice of where to sample VV (left endpoint, right endpoint, or midpoint) differs only at O(ϵ2)O(\epsilon^2) and washes out in the continuum limit. It starts to matter for velocity-dependent or singular potentials, and for the phase-space (Hamiltonian) form where operator ordering ambiguities live. For the standard V(x)V(x) problems in this chapter, the midpoint rule is a convenient symmetric choice and nothing depends on it.

How does this "sidestep" the negative-energy problems of relativistic wave equations?

By changing what is fundamental. In the relativistic single-particle picture, negative-energy plane waves and indefinite densities are unavoidable. The path integral over field configurations promotes the wavefunction to an operator field, and the negative-energy solutions reappear reinterpreted as antiparticles — with a positive-definite, causal amplitude. The construction is the same ∫DeiS/ℏ\int\mathcal De^{iS/\hbar} with SS a spacetime integral of a Lagrangian density; the later chapter on field quantization builds it explicitly. The path integral does not "fix" single-particle relativistic QM so much as show you were asking a field-theory question all along.

Was Dirac's 1933 remark really the origin, or is that a myth?

It is genuine, and Feynman was scrupulous about it. Dirac's 1933 paper observed that the short-time amplitude "corresponds to" eiS/ℏe^{iS/\hbar}, using that cautious phrase without committing to equality. Feynman, prompted by the remark at a Princeton beer party (as the story goes), asked what "corresponds to" precisely meant, worked out that it was a strict proportionality with a definite measure, and built the whole formulation from there. So the credit is layered: Dirac saw the exponential-of-action; Feynman turned "corresponds to" into "equals, with this normalization."

Before this

Where this leads