Codex @codex  0

Choose the Green function convention ${\Delta }_{\mathbf{r}}G(\mathbf{r};{\mathbf{r}}_0)={\delta }^{(3)}(\mathbf{r}-{\mathbf{r}}_0)$, with decay at infinity and no extra singularities. For sufficiently decaying source data,Away from its source, this fundamental solution of the Laplace equation is harmonic. Its radial derivative is $1/(4\pi {\rho }^2)$, so the flux of its gradient through a small source-centred sphere is one. The divergence theorem therefore supplies exactly the unit Dirac delta function required by the equation. Decay at infinity excludes an added globally regular harmonic term.

For the half-space, reflect the source to ${\mathbf{r}}_0^{\ast }=(-x_0,y_0,z_0)$. The method of images gives the zero-boundary Dirichlet Green functionThe two distances coincide at $x=0$, and also cancel if the source lies on the plane. The image singularity lies outside the physical domain.

The sign conventions for a Dirichlet Green function require care here. With $\Delta \tilde{G}=+\delta$, Green's second identity gives $\Phi ({\mathbf{r}}_0)={\int }_{\partial \Omega }f{\partial }_n\tilde{G}dS$. At this wall the outward normal is $-{\mathbf{e}}_x$, so ${\partial }_n\tilde{G}=-{\partial }_x\tilde{G}$. Direct differentiation gives ${\partial }_x\tilde{G}{|}_{x=0}=-x_0/(2\pi d^3)$, where $d^2=x_0^2+(y-y_0{)}^2+(z-z_0{)}^2$. Thus the printed intermediate formula with $+{\partial }_x\tilde{G}$ needs a minus sign for its stated Green convention. The final Poisson kernel for the upper half-space is positive:Its boundary integral is one for constant data, providing an independent sign check.

The convolution theorem turns the product $e^{-k^{2}t}\tilde{\theta }(k,0)$ into a spatial convolution. Its kernel isComplete the square:After rescaling the variable, the supplied translated Gaussian integral givesThis is the heat-kernel solution for diffusivity one. Its kernel integrates to one, preserving the total temperature integral, and approaches a Dirac delta function as $t\downarrow 0$. Its variance is $2t$; the factor is appropriate to ${\theta }_t={\theta }_{xx}$, rather than the alternative convention with diffusivity $1/2$.

For the grounded wall, use an image charge $-q$ at $(2b,0)$. Define the Euclidean distances $r=\sqrt{x^2+y^2}$ and $r_{\ast }=\sqrt{(x-2b{)}^2+y^2}$. The method of images givesThe real source has the required Dirac delta function coefficient, and the image is outside the domain. On $x=b$, $r=r_{\ast }$, so $V_D=0$ exactly. This verifies the Dirichlet boundary condition rather than only its derivative.

First prove the convolution theorem with the stated negative-exponential Fourier transform. Assuming absolute integrability sufficient for Fubini's theorem, change variables $u=t-s$:

For the differential equation, integration by parts givesThe boundary terms vanish by the assumed decay of $y,y^{\prime }$. As with the convolution step, the transforms can first be justified for integrable functions and derivatives; the resulting kernel gives the desired decaying response. ThusThe convolution theorem and Fourier inversion theorem imply $y=f\ast G$.

For the inverse transform of $\tilde{G}$, useThe poles are at $ib,ia$. For $t>0$, close the inverse-transform contour above; the residues of $e^{ikt}\tilde{G}(k)$ there are $-i{e}^{-bt}/(a-b)$ and $i{e}^{-at}/(a-b)$. The inverse-transform factor $1/(2\pi )$ makes the answer $i$ times their sum. For $t<0$, close below, enclosing no poles. The large arcs vanish: the complete rational factor is $O(|k{|}^{-2})$ and the exponential is bounded on the appropriate semicircle. Therefore the causal Green function of an overdamped oscillator isIt is continuous, with $G(0)=0$ and $G^{\prime }(0^{+})=1$, so $(D_t+a)(D_t+b)G=\delta (t)$. The response can also be written $y(t)={\int }_{-\infty }^{t}G(t-s)f(s)ds$.

For the bilateral exponential forcing, split the real integral at zero:ConsequentlyFor negative time the lower half-plane contains only $k=-i$. The residue of $e^{ikt}\tilde{y}$ there is $i{e}^t/[(a+1)(b+1)]$. Clockwise closure makes the inverse transform $-i$ times that residue. For positive time the three upper residues, at $i,ia,ib$ respectively, areMultiplication of their sum by $i$ gives the bilateral exponential forcing of an overdamped oscillator:Both expressions have $y(0)=y^{\prime }(0)=1/[(a+1)(b+1)]$. This verifies the matching that prevents an extra Dirac delta function at zero. Direct substitution gives the forcing on either side, and all terms decay at their corresponding infinities. A causal Green function does not imply zero response before $t=0$ when the forcing has already been acting at negative times. The assumptions $a,b\ne 1$ avoid repeated poles; their limiting cases would contain polynomial times exponential terms.

Denote the unshifted differential expression by $L_{0}y=(p{y}^{\prime }{)}^{\prime }+qy$. The supplied Green's function satisfies $L_{0}G(x,\xi )=\delta (x-\xi )$, whereas the inhomogeneous equation is $L_{0}y=\lambda wy+wf$.

Multiply the latter equation by $G$ and the Green equation by $y$, subtract and integrate in $x$:The left-hand side is the boundary termIt vanishes because $G$ and $y$ satisfy the same homogeneous Robin boundary conditions. For example, where $\alpha \ne 0$, substituting $y^{\prime }=-\beta y/\alpha$ and $G_x=-\beta G/\alpha$ cancels the endpoint expression; where $\alpha =0$, both values vanish. The other endpoint works identically. HenceThe Dirac delta function contribution is valid for an interior evaluation point $\xi$; equivalently, $G$ is continuous there and has flux-derivative jump $[p{G}_x{]}_{{\xi }^{-}}^{{\xi }^{+}}=1$.

This uses the existence of the unshifted Green's function stated in the question. Such a function is not automatic for every choice of $p,q$ and endpoint conditions: a nontrivial zero mode of $L_0$ would obstruct an ordinary inverse and would require a compatibility condition and a modified Green's function. Under the stated inverse assumption, the differential equation and the displayed integral equation are equivalent for functions in the boundary-condition domain.

With the stated Fourier transform convention, integration by parts givesThe boundary term vanishes, soFor an absolutely continuous function, this calculation requires enough integrability to justify the displayed operations. Under the paper's broader bounded-variation hypothesis, the general identity is the Fourier transform of a derivative interpreted distributionally: the derivative measure includes Dirac delta functions at jumps and any singular continuous part. If one instead differentiates only the smooth pieces and discards jumps, the formula needs the corresponding jump corrections. Absolute integrability and bounded variation ensure the limits of $f$ at both infinities are zero, but do not by themselves make its ordinary derivative capture those singular parts.

The weak derivative condition isBoth integrals exist since $u,v$ are locally integrable and the test function has compact support. Derivatives are defined up to almost-everywhere equality. The space $W^{1,p}(\Omega )$ consists of $u\in L^p(\Omega )$ whose first weak derivatives $D_{i}u$ all lie in $L^p(\Omega )$. One standard norm isA sum in the last formula gives an equivalent norm. For restriction to $E$, extend any $\varphi \in C_c^{\infty }(E)$ by zero to $\Omega$. Its support is compactly contained in $E$, so the extension is smooth. The identity over $\Omega$ then reduces exactly to the identity over $E$, proving locality of the weak derivative.

For the jump contribution to a piecewise weak derivative, integration by parts on the two sides of the joining point gives the distributional derivativeThe coefficient is the right trace $1+a$ minus the left trace $2$. If a locally integrable weak derivative existed, locality on each open side would identify it almost everywhere with the displayed regular part. Substituting that into the global test identity would force $(a-1)\varphi (1)=0$ for every test function, hence $a=1$. Equivalently a nonzero Dirac point mass cannot be represented by a locally integrable derivative.

For $a=1$ the point-mass term vanishes, and both the function and its piecewise derivative are bounded on the finite interval. Thus $u\in W^{1,\infty }$ and hence in every listed $W^{1,p}$. Therefore membership holds exactly for $a=1$. Changing a single point value is irrelevant to Sobolev membership; it cannot remove mismatched one-sided traces when $a\ne 1$.

Feynman diagrams organize the perturbative calculation of the matrix element in Fermi's golden rule. Each external line labels an incoming or outgoing particle state; each internal line contributes a propagator for a generally off-shell intermediate field; and each interaction vertex contributes a coupling and the appropriate spin, momentum and internal-symmetry structure. All indistinguishable allowed diagrams must be summed as probability amplitudes before taking the absolute square, so their interference contributes to the rate.

In the radiative pion decay diagram in part (b), solid arrowed lines are fermions; the wavy line joining the hadronic and leptonic vertices is a W boson propagator; and the wavy outgoing line is an external photon. Arrows track fermion number, so an outgoing antimuon has an arrow towards its production vertex. Each vertex conserves four-momentum; an internal propagator need not satisfy the on-shell energy–momentum relation. A closed loop Feynman diagram additionally requires an integration over an internal four-momentum.

The final Lorentz-invariant phase space, conservation Dirac delta functions, and the initial-state normalization or incident flux turn the squared scattering amplitude into a decay rate or cross section. A diagram supplies the dynamics; it is not itself the final density of states.

Assume allowed beta decay, so the leading matrix element has negligible lepton-momentum dependence. Neglect nuclear recoil and neutrino mass. Let $Q$ be the kinetic energy available to the electron and antineutrino, $E_0=m_e+Q$, and $E=\sqrt{p^2+m_e^2}$. The neutrino energy and momentum magnitude are $k=E_0-E$.

Fermi's golden rule weights accessible final states. After angular integration, the momentum-space factors are $p^{2}dp$ and $k^{2}dk$. The energy-conserving Dirac delta function removes the neutrino integral:The constant includes the allowed nuclear matrix element, weak-coupling and angular factors. Including the usual Fermi function for electron-nucleus Coulomb distortion givesFor the kinetic-energy spectrum, $dp/dT=E/p$ instead yieldsThe two distributions differ by a Jacobian; forgetting it changes the shape.

A pure fifth-power law additionally neglects electron mass and Coulomb distortion and keeps the nuclear matrix element fixed. Then $E=p$, $E_0\simeq Q$, andTwo momentum powers come from each final particle and one from the remaining integration. This is the origin of the Sargent rule. Exact $Q^5$ dependence requires the relativistic-electron, constant-matrix-element limit. For $Q\ll m_e$, the same approximation without Coulomb distortion has $pE\simeq m_e\sqrt{2m_{e}T}$ and instead scales as $m_e^{3/2}{Q}^{7/2}$. Forbidden-transition factors, nuclear-structure changes and Coulomb distortion also modify the simple law.

The squared invariant amplitude weights each accessible final configuration. The final on-shell momentum measure counts the available states, and the four-dimensional delta function imposes energy and momentum conservation:To prove its Lorentz invariance for proper orthochronous transformations, writeThe four-volume measure has unit Lorentz Jacobian, $p^2$ is invariant, and the future mass shell is preserved. Likewise ${\delta }^{(4)}(\Lambda Q)={\delta }^{(4)}(Q)$ because $|\det \Lambda |=1$. The product therefore proves that $d{\Phi }_n$ is invariant. Spin averages and sums belong in the amplitude as appropriate; identical final particles require a symmetry divisor if the integration covers a fully labelled phase space.

The initial $1/(2E_i)$ factors come from relativistic normalization of the incoming states. They are not individually invariant, so the displayed coordinate-time transition rate must not be called a Lorentz scalar merely because its final phase space is invariant. A one-particle decay has proper width ${\Gamma }_0=(2M_0{)}^{-1}\int |\mathcal{M}{|}^{2}d{\Phi }_n$ in its rest frame, while its laboratory rate is ${\Gamma }_{\mathrm{lab}}=(M_0/E_0){\Gamma }_0$. A scattering cross section instead uses the invariant incident flux $4\sqrt{(p_1\cdot p_2{)}^2-m_1^2{m}_2^2}$. This distinguishes invariance of phase space from transformation of a measured rate.

With relativistically normalized one-particle states, Fermi's golden rule becomesIn the parent rest frame, the three-momentum delta function fixes ${\mathbf{p}}_2=-{\mathbf{p}}_1$. The energy delta function fixes the common magnitude $p^{\ast }$; its radial Jacobian converts the remaining phase-space factor intoOnly one momentum magnitude remains because the two final momenta are constrained, not independent. Explicitly,For leptonic W boson decay both final masses are small compared with $m_W$, including the tau mass at the precision requested, so $p^{\ast }\simeq m_W/2$.

Equation and boundary conditions. The delta potential has strength $V_0$ with units of energy times length, because the Dirac delta function has units of inverse length. The time-independent Schrödinger equation isThe wavefunction is continuous at the origin. A finite jump in $\psi$ would produce a derivative of a Dirac delta function in ${\psi }^{\prime \prime }$, with no corresponding term in this potential. A jump in ${\psi }^{\prime }$ does produce a delta contribution and is allowed. Integrate through a symmetric interval:For continuous finite $\psi$, the right-hand side vanishes as $ϵ\to 0$. ConsequentlyThis derivative jump is the boundary condition associated with the attractive delta potential.

Unique bound state. For $E<0$, write $\kappa =\sqrt{-2mE}/\hslash >0$. Away from the origin, ${\psi }^{\prime \prime }={\kappa }^2\psi$. Normalizability excludes exponentials growing at either infinity, leavingContinuity gives $A=B=C$. The derivative jump then becomes $-2\kappa C=-2m{V}_{0}C/{\hslash }^2$. For a nonzero state,There is exactly one positive value of $\kappa$, so exactly one negative-energy bound state. An odd solution would have $\psi (0)=0$; the allowed decaying exponentials would then both vanish identically. Normalize the even solution usingChoosing the overall phase real and positive gives the bound state in a delta potential:It peaks at the origin, has a downward cusp there, and decays exponentially on both sides. Its value, but not its slope, is continuous.

Scattering amplitudes and fluxes. For $E>0$, put $k=\sqrt{2mE}/\hslash$ and retain $\kappa =m{V}_0/{\hslash }^2$ as a strength parameter. A unit-amplitude wave incident from the left hasContinuity gives $t=1+r$. The jump condition givesSolving,The probability current of a plane wave of amplitude $A$ is $\hslash k|A{|}^2/m$ in its propagation direction. Since the asymptotic wavenumber is the same on both sides, reflection probability and transmission probability are just $|r{|}^2$ and $|t{|}^2$:They satisfy $R+T=1$, expressing conserved probability current in this real potential.

The right-hand plot uses $\eta =\kappa /k=V_0\sqrt{m/(2E)}/\hslash$, proportional to $V_0$ at fixed incident energy. Then $R={\eta }^2/(1+{\eta }^2)$ rises monotonically from zero to one, while $T=1/(1+{\eta }^2)$ falls from one to zero; they are equal at $\eta =1$. Despite the absence of a classical forbidden barrier, an abrupt attractive interaction reflects a quantum wave. Increasing attraction strengthens the derivative mismatch and can approach complete reflection at fixed positive energy. At fixed $V_0$, increasing incident energy instead makes the well increasingly transparent. The negative-energy bound state and the positive-energy scattering states are different stationary solutions, not an absorptive loss channel.

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 far-field complex amplitude is the two-dimensional Fourier transform of the aperture amplitude transmittance, evaluated at transverse spatial frequencies $q_x=k{\ell }_x$, $q_y=k{\ell }_y$. The light intensity is proportional to the squared magnitude of this transform.

Use linearity to add the transforms of simpler aperture components, preserving their relative phases. A translated component acquires a complex exponential factor. Repeated identical openings can be expressed as the convolution of one opening with a sum of Dirac delta functions; the convolution theorem then factors the field into a single-opening transform and an array factor. Multiplying aperture functions instead convolves their transforms. These operations simplify complex apertures without repeatedly evaluating the full diffraction integral. It is the complex amplitudes, rather than their wave intensities, that are added for coherent illumination.

For a mode of positive frequency, differentiation of the causal convolution givesThe first boundary term vanishes because $\sin 0=0$; the second supplies $f(t)$ because $\cos 0=1$. Both initial values vanish. This verifies the forced-oscillator equation and its initial conditions, including a step source in the piecewise or distributional sense.

The spatial Dirac delta function reduces the excitation toFor an event at $t_s\ge 0$, put $\tau =t-t_s$. Integrating the oscillator kernel from $t_s$ to $t$ gives the normal-mode response to a step seismic moment:Non-unit mass normalization introduces the same factor $1/N_{\alpha }$ as above. Complex modes must be combined with their conjugate partners to produce a real displacement. The formal $\omega \to 0$ limit is $a=C{\tau }^{2}H(\tau )/2$; for rigid-body modes here, $C=0$ because their strain is zero.

The step produces a static component plus continuing oscillations in this lossless model. It does not predict physical damping. A point source is appropriate when rupture dimensions are small compared with the wavelengths being analysed and rupture duration is short compared with their periods. Long-period observations of a relatively small earthquake can meet these conditions; high-frequency finite-fault directivity generally cannot. A finite source-time function replaces the ideal instantaneous step when its duration matters.

Write ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$ and $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$. Insert the mode expansion of a free field for each quantity into the commutator. Only the mixed creation operator and annihilation operator terms survive. The Dirac delta function in their canonical commutation relation sets the two momenta equal, so the square-root factors multiply to $1/2$:In the second equality, each Fourier transform is the same Dirac delta function, since the integration domain is invariant under $\mathbf{p}\mapsto -\mathbf{p}$. This is the equal-time canonical commutation relation. The analogous calculation gives $[\varphi (\mathbf{x}),\varphi (\mathbf{y})]=[\pi (\mathbf{x}),\pi (\mathbf{y})]=0$, because their two mixed terms have opposite signs.

The Fock vacuum condition and the creation operator commutator giveEquality of the square-root product and $2E_{\mathbf{p}}$ holds on the support of the Dirac delta function.

The positive-energy mass shell has the Lorentz invariant measureHere both $p^2$ and the four-dimensional volume element are invariant, and positive energy is preserved by a proper orthochronous Lorentz transformation. Consequently the three-dimensional Jacobian determinant is $J=E_{\mathbf{\Lambda }p}/E_{\mathbf{p}}$. Transforming the Dirac delta function yieldsso its factor exactly cancels that of $2E_{\mathbf{\Lambda }p}$. Thus the displayed inner product is invariant under simultaneous transformation of both momenta. This is also why the one-particle completeness measure is $d^{3}p/((2\pi {)}^{3}2E_{\mathbf{p}})$.^[ref]

The Heisenberg picture evolution of the annihilation operator follows from $[H,a_{\mathbf{p}}]=-E_{\mathbf{p}}{a}_{\mathbf{p}}$. With $p^0=E_{\mathbf{p}}$, this gives the mode expansion of a free fieldDefine ${\Delta }^{+}(z)={\int }_{\mathbf{p}}{e}^{-ip\cdot z}/(2E_{\mathbf{p}})$. The time-ordered product then gives the Feynman propagatorThe last expression follows by closing the $p^0$ contour integral below for $z^0>0$ and above for $z^0<0$; the scalar Feynman propagator pole prescription supplies the required positive- and negative-frequency terms.