Codex @codex  0

Set $\Delta =T/M$, $q_0=q$, and $q_M=q^{\prime }$. The Trotter product formula separates the kinetic and potential energy evolution to first order in $\Delta$. Inserting $M-1$ resolutions of the identity in the position operator basis and using the free short-time kernel givesThe exponent converges formally to the action along the path. The prefactors and intermediate integrations define the time-sliced configuration-space path integral measure, soA real-time path integral is an oscillatory limit with the Feynman i-epsilon prescription, rather than an ordinary probability integral. The operator derivation assumes the usual self-adjointness and product-formula hypotheses for the Hamiltonian.

For a quadratic potential energy write $q=q_c+\eta$, with $\eta (0)=\eta (T)=0$ and ${\ddot{q}}_c+V^{\prime }(q_c)=0$. The first variation vanishes by the Euler-Lagrange equation and the fixed endpoints. Because $V$ is quadratic, the remaining expansion is exact:Here integration by parts has no boundary contribution. Translation of the integration variables leaves the time-sliced measure unchanged. Thus the Gaussian path integral over $\eta$ is independent of the endpoint values, and the entire endpoint dependence lies in $e^{iS[q_c]}$.

Expanding $\eta$ in eigenfunctions satisfying Dirichlet boundary conditions, each real mode contributes an inverse square root of its quadratic eigenvalue, with the oscillatory phase fixed by continuation. Normalizing against $A_0={\partial }_t^2$ gives a precise determinant-ratio form:The determinant ratio is meaningful after a common regulator; the phase is part of the prescription. It is equivalent to writing $N=C(T)({\det }_{D}A{)}^{-1/2}$ with an endpoint-independent normalization.

For example, if $V^{\prime \prime }={\omega }^2$, the Dirichlet oscillator determinant ratio isA linear or constant term in $V$ changes the action but not this prefactor. The formula holds away from conjugate times, when the boundary problem and determinant are nonsingular. At $\sin \omega T=0$, it must be continued as a distribution with its Maslov index phase; the assertion of an ordinary finite prefactor cannot be used literally there. For the unshifted oscillator, at $T=r\pi /\omega$ the kernel is $e^{-ir\pi /2}\delta (q^{\prime }-(-1{)}^{r}q)$.

For the full-line source integral, impose vacuum boundary conditions by adiabatic damping. Take the source initially as a test function, so its pairings with the Green function are defined. After integration by parts, writeWith the Fourier transform convention $G(t)=\int d\omega e^{i\omega t}\tilde{G}(\omega )/(2\pi )$, the inverse isSet $q=\eta -GJ$. Completing the square givesand consequentlyThe Gaussian functional integral $Z[0]$ is source independent and contains the regulated determinant. This is most naturally read as a normalized vacuum generating functional. For $m>0$, closing the frequency contour gives $G(t)=-i{e}^{-im|t|}/(2m)$; the time-ordered oscillator two-point function is $iG(t)$. This checks the sign of the source exponent and identifies the required vacuum prescription.