Codex @codex  0

With the mode normalization given, the oscillator canonical commutation relations are $[a_{\mathbf{p}},a_{\mathbf{q}}^{†}]=(2\pi {)}^3{\delta }^3(\mathbf{p}-\mathbf{q})$, with the other two commutators zero. Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p^0=E_{\mathbf{p}}$. The convention compatible with the requested numerator iswithout an additional factor of $i$ outside the vacuum expectation value. Put $\tau =x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$. The Fock vacuum is annihilated by $a_{\mathbf{p}}$, so only $a{a}^{†}$ contributes to the two Wightman functions. Consequently,In the second term we changed $\mathbf{p}\mapsto -\mathbf{p}$ to give the same spatial exponential.

Now perform the energy contour integralHere the residue step uses $E>0$; the massless zero-momentum point is interpreted through the smeared distribution limit, not as an isolated normalized oscillator. The Feynman i-epsilon prescription puts the positive-energy pole below the real axis and the negative-energy pole above it: $p^0=E-i0$ and $p^0=-E+i0$. For $\tau >0$, close in the lower half-plane, clockwise; the residue theorem gives $-i$ times the residue $i{e}^{-iE\tau }/(2E)$, hence $e^{-iE\tau }/(2E)$. For $\tau <0$, close in the upper half-plane, counterclockwise; the negative-energy residue is $-i{e}^{iE\tau }/(2E)$, again giving a positive $e^{iE\tau }/(2E)$. These are exactly the two time-ordered terms. Restoring the spatial integral proves the scalar Feynman propagator pole prescription:The prescription in this formula supplies the pole convention left unspecified in the printed display; an unprescribed ordinary real-axis integral would not be well-defined. It is a distribution limit after smearing, not an absolutely convergent four-dimensional integral. With this normalization the derivative jump of a free scalar time-ordered two-point function gives $(\square +m^2){\Delta }_F=-i{\delta }^4(x-y)$, an independent check of both the numerator and the sign.

The pole displacement is exaggerated in this original schematic. The contour orientation and selected pole reproduce the time ordering of the Feynman propagator.

Use the time-harmonic wave convention $\psi e^{-i\omega t}$ and write $\mathbf{z}=(y,z)$, $E(0,\mathbf{z})=E_0(\mathbf{z})$, with a deterministic incident wave envelope. The printed speed ratio is inconsistent with $k=k_{0}n$: the refractive index used below is $n=c_0/c$. Substituting $\psi =e^{i{k}_{0}x}E$ in the Helmholtz equation givesThe paraxial approximation discards ${\partial }_x^{2}E$. The usual weak-fluctuation model also linearizes $n^2-1=2\mu W+O({\mu }^2{W}^2)$, givingDropping the quadratic contrast is an additional weak-fluctuation assumption, not a consequence of small propagation angles. Although the linearized random potential generates attenuation of order ${\mu }^2$, it does not retain every effect of that order in the literal finite-correlation index: the discarded quadratic contrast can also produce a mean wave phase shift. This linearization must precede a Gaussian white noise limit: the square of ideal white noise has no ordinary pointwise meaning.

The split-step Fourier method alternates free-space diffraction, $U_h=e^{hD}$, with a random phase screen,For jointly Gaussian random fields, Gaussian phase averaging gives $\mathbb{E}{e}^{i{\varphi }_h}=e^{-\mathrm{Var}({\varphi }_h)/2}$. More generally, a product of $p$ fields and $q$ conjugate fields picks up $\exp (i{\sum }_a{ϵ}_a{\varphi }_h({\mathbf{z}}_a))$, with signs ${ϵ}_a=+1$ or $-1$. Its screen average is determined entirely by the wave phase covariance matrix. The deterministic diffraction step acts on each coordinate, with opposite signs on factors formed by complex conjugation. This is the basis of the field-moment equations.

There is an important closure qualification. A stationary Gaussian random field need not have independent longitudinal increments. For finite-correlation fluctuations, the unlinearized parabolic wave equation instead givesIn the weak-fluctuation model, the exact first-moment equation isThe last term cannot in general be replaced by a constant times $m$. An exact generic solution of the linearized model iswhere $\mathcal{T}$ denotes time ordering of the propagation operators. For the unlinearized model, add $i{k}_0{\mu }^2{W}^2/2$ inside the propagation generator. The covariance function of the medium is needed to evaluate this expression; One-point Gaussian distributions alone would not even determine the joint wave phase statistics.

For example, omit diffraction and take a longitudinal autocovariance function $R(s)=\mathbb{E}[W(s,\mathbf{z})W(0,\mathbf{z})]$. Direct Gaussian phase averaging givesEven here, unit-variance stationary Gaussian random fields with different $R$ give different answers. In the general problem the diffraction and multiplication operators do not commute, so this scalar attenuation cannot simply be multiplied by $U_x$ without an additional approximation.

The standard closed answer uses the Markov approximation for a random medium, made explicit in part ii. Let the longitudinally integrated covariance kernel beReplace the medium by longitudinal Gaussian white noise with this strength. A screen of thickness $h$ then has $\mathbb{E}[{\varphi }_h({\mathbf{z}}_a){\varphi }_h({\mathbf{z}}_b)]=k_0^2{\mu }^{2}hB({\mathbf{z}}_a-{\mathbf{z}}_b)$ and is independent of the incoming field. For the moment $M=\mathbb{E}{\prod }_a{E}_a^{({ϵ}_a)}$, expanding both steps to order $h$ givesIn particular, the coherent attenuation in a white-noise random medium and its Fresnel propagator solution areEquivalently, write $Wdx=d{\beta }_x$, where ${\beta }_x(\mathbf{z})$ is Brownian motion in $x$ with transverse covariance kernel $\mathbb{E}[d{\beta }_x(\mathbf{z})d{\beta }_x({\mathbf{z}}^{\prime })]=B(\mathbf{z}-{\mathbf{z}}^{\prime })dx$. The Stratonovich integral formulation is $dE=DEdx+i{k}_0\mu E\circ d{\beta }_x$. Its Itô integral form isThe mean of the Itô integral vanishes, independently confirming the attenuation drift. For $x>0$ in two transverse dimensions,Thus a unit plane wave has $m=e^{-\gamma x}$. For a general incident wave envelope, the Fresnel propagator supplies its spreading. The attenuation is redistribution between the coherent and diffuse wave fields, rather than wave absorption: for a field and its conjugate at the same point, the screen contribution in the second-moment equation cancels.

For bosonic fields, time ordering places operators with later time arguments to the left:Normal ordering places all creation operators to the left of all annihilation operators and is denoted by colons. A Wick contraction isWick theorem stateswhere each sum runs over inequivalent disjoint pairings and contracted fields are replaced by their Feynman propagator.

The Time-dependent Schrödinger equation $i\hslash {\partial }_t|\psi ⟩=H(t)|\psi ⟩$ governs unitary evolution of a closed system. If $H$ is independent of time, expand in its energy eigenstates:An individual energy eigenstate changes only by a phase, so its probability density and time-independent-observable expectations are stationary. Superpositions can have interference terms oscillating at $(E_n-E_m)/\hslash$; fixed energy populations do not imply that every observable is constant.

For general $H(t)$, $U(t,t_0)=\mathcal{T}\exp [-(i/\hslash )\int H(t^{\prime })d{t}^{\prime }]$ requires time ordering unless Quantum Hamiltonians at different times commute. The equivalent Heisenberg description evolves observables rather than the state, giving $d⟨A⟩/dt=(i/\hslash )⟨[H,A]⟩+⟨{\partial }_{t}A⟩$. Explicit time dependence can transfer energy from an applied perturbation.

For a weak perturbation $V(t)$ to a known $H_0$, the first-order transition amplitude from $n$ to $m\ne n$ isIts squared modulus gives the leading transition probability while depletion remains small. A nearly monochromatic perturbation emphasizes energy differences near its quantum of frequency; a continuum and suitable long-time limit yield Fermi's golden rule. These approximations have stated weak-coupling and time-resolution limits; they do not replace exact unitary dynamics in every regime.

Time dependence in quantum mechanics describes how a prepared quantum state changes and thereby changes the probabilities of later measurements. For an isolated system, the Schrödinger picture assigns this evolution to a state vector $|\psi (t)⟩$, governed by the Time-dependent Schrödinger equationThe Hamiltonian is a self-adjoint operator; it has units of energy, while $\hslash$ has units of action, so $\hat{H}/\hslash$ has units of inverse time. For a particle with potential energy $V(\mathbf{r},t)$, this becomesGiven a suitable initial wavefunction and the physical boundary conditions, this first-order equation determines its subsequent propagation. Its linearity means that the evolution of a quantum superposition is the corresponding quantum superposition of evolved states.

For a real potential energy, this equation and its complex conjugate imply the local continuity equationHere $\rho$ is the probability density and $\mathbf{j}$ the probability current. Vanishing net boundary flux preserves wavefunction normalization, as required by the Born rule.

The time-evolution operator packages the solution for every initial state:Because ${\hat{H}}^{†}=\hat{H}$, differentiation givesTherefore $U^{†}U=I$: this is unitary time evolution, preserving inner products and all state norms. It is reversible, with $U(t,t_0{)}^{-1}=U(t_0,t)=U(t,t_0{)}^{†}$, and obeys the composition law $U(t_2,t_0)=U(t_2,t_1)U(t_1,t_0)$.

For a time-independent Hamiltonian, the solution is the operator exponentialIf the Hamiltonian varies with time, the analogous simple exponential of its integral is valid when $[\hat{H}(t),\hat{H}(t^{\prime })]=0$ for all relevant times. Otherwise time ordering is essential:The meaning is explicit in the first terms of the ordered series,the later-time operator stands on the left. In the position representation, the same propagation is expressed by the quantum-mechanical propagatorIts initial value is the Dirac delta function ${\delta }^3(\mathbf{r}-{\mathbf{r}}^{\prime })$.

A time-independent Hamiltonian has energy eigenstates satisfying $\hat{H}|n⟩=E_n|n⟩$. A single such initial state evolves asIt is a stationary state: its global phase changes, but its probability density and the distributions of time-independent observables do not. Its energy is sharp. Any quantum superposition entirely within one degenerate energy eigenspace also acquires one common global phase and is stationary. Stationarity does not require zero probability current; for instance, $u_n(\varphi )=e^{in\varphi }/\sqrt{2\pi }$ on the ring has a constant density and persistent circulation when $n\ne 0$.

An arbitrary normalized state can instead contain different energy eigenvalues:A continuous spectrum requires corresponding integrals. For the discrete case, an energy measurement gives $E_n$ with fixed probability $|c_n{|}^2$ (summed over degeneracies), but relative quantum phases change with time. ConsequentlyThe interference terms can oscillate at angular frequencies $(E_m-E_n)/\hslash$. For two populated energy levels, this supplies a characteristic period $2\pi \hslash /|E_m-E_n|$ when the interference term is nonzero. Thus fixed energy populations do not imply a stationary spatial distribution.

For a normalized state and an observable $\hat{A}(t)$, its expectation value isThis is the statistical mean over repeated preparations and measurements. In an eigenbasis of energy eigenstates for a time-independent Hamiltonian and a time-independent $\hat{A}$,Only matrix elements connecting different energies can generate these oscillations. If $[\hat{A},\hat{H}]=0$, its expectation value is constant for every initial state, provided it has no explicit time dependence.

Differentiating the expectation value and substituting the Schrödinger equation for both the bra and ket proves the general Ehrenfest theorem:Equivalently,These manipulations require the state to lie in the relevant operator domains. Taking $\hat{A}=\hat{H}$ gives $d⟨H⟩/dt=⟨{\partial }_t\hat{H}⟩$, so energy is conserved under a time-independent Hamiltonian.

For $\hat{H}={\hat{p}}^2/(2m)+V(\hat{x})$, the canonical commutation relation givesHence the particularly useful position-and-momentum form of Ehrenfest theorem isIt resembles Newton's second law, but the mean force generally differs from the force at the mean position. If a narrow wave packet samples a region where the force varies little, $⟨V^{\prime }(x)⟩\simeq V^{\prime }(⟨x⟩)$ gives approximate classical motion. For a quadratic potential energy this replacement is exact. In particular the quantum harmonic oscillator obeysThe centroid follows the classical harmonic oscillator, although this does not specify the full evolving wavefunction or its width.

The Heisenberg picture expresses the same physics by assigning evolution to observables instead:The expectation value is unchanged, and differentiation yields $d{\hat{A}}_H/dt=(i/\hslash )[{\hat{H}}_H,{\hat{A}}_H]+U^{†}({\partial }_t{\hat{A}}_S)U$. Thus state propagation, evolving expectation values, and the classical correspondence in Ehrenfest theorem are mutually consistent descriptions of the same unitary time evolution.

The time-evolution operator maps an initial quantum state to its state at a later time. For a time-independent Hamiltonian operator $H$, it is the unitary operator $U(t)=e^{-iHt/\hslash }$.

More generally, write $|\psi (t)⟩=U(t,t_0)|\psi (t_0)⟩$. The time-evolution operator obeysFor a constant self-adjoint Hamiltonian operator, $U(t,t_0)=\exp [-iH(t-t_0)/\hslash ]$. With zero as the initial time, this gives the group law $U(t+s)=U(t)U(s)$. Successive evolution is reversible; a measurement followed by conditioning on an outcome is a separate operation described by projective measurement.

For a time-dependent Hamiltonian operator, the defining equation is $i\hslash {\partial }_{t}U(t,t_0)=H(t)U(t,t_0)$. In finite dimension, a continuous Hermitian $H(t)$ still gives unitary time evolution, sinceThe solution is a time-ordered exponential, with later factors acting on the left. One may replace it by $\exp [-i{\int }_{t_0}^{t}H(s)ds/\hslash ]$ when $[H(s),H(s^{\prime })]=0$ at all times; independently changing pulse axes generally requires time ordering.^[ref]