Codex @codex  0

For ordinary bosonic operators at distinct times, time ordering places the later operator to the left:Equivalently, multiply these two products by the corresponding Heaviside step functions. For two odd fermionic operators the second product receives a minus sign; more generally a graded interchange has sign $(-1{)}^{F_1{F}_2}$. This sign is the fermionic sign convention for time-ordered field products, not a minus sign for arbitrary bosonic operators. No equal-time convention is needed for the two distinct times in this clause.

The solution and its first derivative are continuous at $x=1$: a jump would create a Dirac delta distribution or its derivative, absent from the forcing. Put $r=x-1$, $A=a+1$ and $B=b$. The outer expansions, fixed by the respective endpoint boundary conditions, areBalancing diffusion and advection near $r=0$ gives the interior layer at a simple zero of advection with inner variable $z=r/\epsilon$. The leading inner expansion solves $Y_0^{\prime \prime }+2z{Y}_0^{\prime }=0$, so matching to $A$ and $B$ givesFor the next term define the Gaussian drift primitivesWriting $H$ for the Heaviside step function, let $F(z)=zH(-z)+(\sqrt{\pi }/4)\mathrm{erf}|z|$. Its value and first derivative match at zero and $(D^2+2zD)F=2zH(-z)$. Thuswhere the logarithmic overlaps fixHere $C_1,C_2$ are the constants in the supplied large-positive-$z$ limits of $E_0,E_1$. The terms $\epsilon \log \epsilon$ are essential switchback terms; discarding them would not give accuracy through $O(\epsilon )$.

Subtracting the common overlap from the inner expansion and outer expansions gives the additive composite expansionThe standard overlap subtraction leaves $\epsilon rH(-r)$. Replacing $H(-r)$ by $(1-\mathrm{erf}z)/2$ changes the value only by $O({\epsilon }^2)$ uniformly and gives a composite with continuous first derivative. The last term restores the endpoint values through the retained order. When $b=a+1$, $d=0$ and the order-one error-function jump disappears. An $O(\epsilon )$ interior adjustment remains to accommodate the leading outer derivative mismatch, together with logarithmic matching when $A\ne 0$.

Reversing the diffusion sign changes the leading inner equation to $-Y_0^{\prime \prime }+2z{Y}_0^{\prime }=0$. Its nonconstant solution grows like the integral of $e^{z^2}$, so bounded matching forces the same leading value $C$ on both sides. The bulk is $C+r$ on the left and $C$ on the right; both endpoint conditions are instead supplied by decaying endpoint layers for reversed diffusion, of width ${\epsilon }^2$. A weaker width-$\epsilon$ interior adjustment matches derivatives and selects $C$. Indeed, its first-order equation has a non-growing solution only ifThe leading structure is consequentlywith the smaller interior correction understood. Unlike the positive-diffusion case, the endpoint values are not transported into two distinct order-one inner limits.

For $t>2$, the Heaviside step function contributes a constant unit forcing. Solving the resulting first-order linear differential equation from the value at $t=2$ gives

Since $(x^2-1)(-{\delta }_0)={\delta }_0$, write $u^{\prime }=-{\delta }_0+w$, where $(x^2-1)w=0$. The polynomial has two simple zeros. The distribution annihilated by a function with simple zeros lemma givesIntegrating with the Heaviside step function yieldsThis includes every solution: its derivative has exactly the required delta terms, and any two primitives differ only by a constant distribution.

For $f(t)=\sin \omega t$,The Bromwich inversion formula givesWhen $t>x/c$, close the contour to the left. The residues at $p=\pm i\omega$ sum to $\sin (\omega (t-x/c))$. When $t<x/c$, close it to the right, where there are no poles, and the integral vanishes. HenceEquivalently, $y(x,t)=\Theta (t-x/c)\sin (\omega (t-x/c))$ in terms of the Heaviside step function. The boundary oscillation travels to the right at wave speed $c$ without reflection or dispersion; a point at $x$ begins moving after the propagation delay $x/c$.

The causal Green function for the heat operator iswhere $H$ is the Heaviside step function. It vanishes for $t<\tau$, satisfies the homogeneous heat equation away from the source, and approaches $\delta (x-\xi )$ as $t\downarrow \tau$. Superposing the responses to all infinitesimal sources gives Duhamel principle:The lower limit and causality give the required homogeneous initial data.

After separating the conserved horizontal factor as $\psi =e^{ipx}f(z)$, the Helmholtz equation becomeswhere $H$ is the Heaviside step function. The outgoing Green function for $d^2/d{z}^2+q_0^2$ isThe Born approximation at first order replaces $f$ on the right-hand side by the incident profile $f_i(z^{\prime })=e^{-i{q}_0{z}^{\prime }}$. For $z>0$ this givesThe integral is understood with the usual outgoing-wave convergence factor. Sincewe obtainTherefore the Born reflected field is

Letwhere $\mathcal{H}$ is the Heaviside step function. Thenso the Eliassen equation for residual circulation isTake no normal residual flow at $y=0,L$ and $z=0$, decay as $z\to \infty$, and choose the streamfunction constant on the connected rigid boundary to be zero. Thus

The required Fourier series in sine modes isDefineFor each mode, the vertical equation isThe Dirac delta function requiresThe solution satisfying the boundary conditions is thereforewhere

The momentum equation givesThe jump of ${\overline{\chi }}_a^{\ast }$ supplies an equal positive delta function in $f_0{\overline{\chi }}_{a,z}^{\ast }$, so the singular terms cancel. The regular acceleration iswhereFinally, the transformed density equation givesThese exponentially decaying modes are the balanced mean response to wave-activity deposition at $z=H$.

Laplace transform. Use $Y(p)={\int }_0^{\infty }{e}^{-pt}y(t)dt$. The initial value is zero and the Dirac delta contributes $e^{-p}$, so the transformed equation is $pY+Y=e^{-p}$. Therefore

Bromwich inversion. The Bromwich inversion formula iswith the vertical Bromwich contour traversed upwards and its symmetric improper limit understood. For $t>1$, close to the left, where the exponential decays; the enclosed simple pole at $p=-1$ has residue $e^{-(t-1)}$. For $t<1$, close to the right, which contains no pole. Thuswhere $\Theta$ is the Heaviside step function. This agrees with the Laplace transform time-shift rule. At $t=1$, inversion returns the average $1/2$ of the two one-sided values. The value assigned at that single point does not affect the differential equation in the distributional sense.

Fourier transform. State the convention explicitly:The Fourier transform of the distributional derivative is $ik\tilde{y}(k)$, including the jump at $t=1$. The Dirac delta has transform $e^{-ik}$, giving

Inverse Fourier integral. Substitute into the inverse Fourier transform:The only pole is at $k=i$. For $t>1$, close in the upper half-plane. Jordan's lemma removes the large arc and the residue theorem gives $i{\mathrm{Res}}_{k=i}=i[e^{-(t-1)}/i]=e^{-(t-1)}$. For $t<1$, close in the lower half-plane, where there are no poles, giving zero. This reproducesIntegrating the differential equation across the impulse checks the unit jump, $y(1^{+})-y(1^{-})=1$. Away from the impulse the solution obeys $y^{\prime }=-y$, and it satisfies the stipulated zero past.

Now the inverse integrand can be writtenFor $0<\lambda <1/2$, the exponent has a positive real time parameter. Close the Bromwich contour to the left; the residue at $s=-1$ is $\frac{1}{2}(2\lambda )=\lambda$. For $1/2<\lambda <1$, close to the right, where there are no poles. The reversed orientation there still gives zero. ConsequentlyThe same result follows from the Laplace transform time-shift rule: $H(s)=\frac{1}{2}{e}^{-s\log 2}/(s+1)$ corresponds to $f(t)=e^{-t}\Theta (t-\log 2)$, where $\Theta$ is the Heaviside step function.

The original transform cannot fix a value at the single point $\lambda =1/2$. The Bromwich inversion formula chooses the average of the limiting values $1/2$ and zero, giving $\boxed{h(1/2)=1/4}$ under that convention. Indeed,regardless of the assigned value at the jump.

For causal functions the convolution is $(f\ast g)(t)={\int }_0^{t}f(\tau )g(t-\tau )d\tau$. In a half-plane of absolute convergence, changing variables $u=t-\tau$ in its Laplace transform and using Fubini's theorem givesThis proves the convolution theorem for Laplace transforms. Transforming the differential equation gives $(s^2+1)Y(s)=sy(0)+y^{\prime }(0)+F(s)$. The direct transform pairs $\mathcal{L}(\cos t)=s/(s^2+1)$ and $\mathcal{L}(\sin t)=1/(s^2+1)$ then giveFor a continuous forcing, call the last term $w(t)$. The Leibniz integral rule gives $w^{\prime }={\int }_0^{t}f(\tau )\cos (t-\tau )d\tau$ and $w^{\prime \prime }=f(t)-w(t)$. Also $w(0)=w^{\prime }(0)=0$. Thus both the original equation and the initial values hold.

For an impulse at positive time $T$, the Dirac delta evaluates the convolution, yieldingHere $H$ is the Heaviside step function. The displacement remains continuous at $T$, while its derivative jumps upward by one, exactly as required by integrating the differential equation across the unit impulse.

Take $x\ge 0$ and a source position $\xi >0$. Before the source, the zero initial data require $G=0$. After it, set $G=e^{-2(x-\xi )}h(x-\xi )$; the homogeneous equation becomes $h^{\prime \prime }+a^{2}h=0$. Continuity of $G$ and a unit jump in $G_x$ at $x=\xi$ give $h(0)=0$, $h^{\prime }(0)=1$. Thus the causal Green function of a damped oscillator isHere $H$ is the Heaviside step function. For $a=0$, take the continuous limit $\sin (as)/a\to s$, giving $G=H(x-\xi )(x-\xi )e^{-2(x-\xi )}$. The formula is also unchanged when a nonzero $a$ changes sign. Its continuity removes a delta derivative, and its unit first-derivative jump supplies the required Dirac delta function.

For the forced equation use $a=2$. The zero-initial-data Green function representation isThe sine-product identity reduces this to elementary exponential-trigonometric integrals. Evaluating givesFor an independent check, the steady particular part is $\sin 2x/20-\cos 2x/10$, while the decaying homogeneous part has coefficients $1/10$ and $1/20$ that impose $y(0)=y^{\prime }(0)=0$.

Use $f(x)=(2\pi {)}^{-1}{\int }_{\mathbb{R}}\tilde{f}(k)e^{ikx}dk$. Jordan's lemma says that for $x>0$, if a function $h$ has maximum modulus $M_R\to 0$ on upper semicircles, then ${\int }_{C_R}h(k)e^{ikx}dk\to 0$; a bound is $\pi M_R/x$. For $x<0$ use lower semicircles instead. The rational factors here meet this decay condition, so the residue theorem evaluates their inverse transforms.

For $1/(a+ik)$, only the upper half-plane contains a pole, at $k=ia$. Its residue after multiplication by $e^{ikx}$ is $-i{e}^{-ax}$. The lower contour encloses no pole. ThereforeHere $H$ is the Heaviside step function, with midpoint value $H(0)=1/2$ for symmetric Fourier inversion.

For $1/(a^2+k^2)$, the upper-pole residue is $e^{-ax}/(2ia)$ when $x>0$, while the clockwise lower contour gives the matching expression for $x<0$. ThusThe final rational factor has double poles at $\pm ia$. At the upper pole,The upper inverse integral is $i$ times that residue. The lower contour for negative $x$ supplies $-x{e}^{ax}/2$, and the value at zero is zero. Consequently

Use the zero initial displacement and velocity given in the PDF. The characteristic roots of the homogeneous equation are $-1\pm i\mu$. Its causal Green function of a damped oscillator iswhere $H$ is the Heaviside step function. The positive-time kernel has $G(0+)=0$ and $G^{\prime }(0+)=1$, so its derivative jump generates a unit impulse. The zero-data forced solution is consequently the causal convolutionDifferentiating once gives $x^{\prime }(0)=0$; differentiating twice supplies the endpoint term $g(t)$ and verifies the differential equation. No homogeneous transient is required with these initial data.

The result of part (i) vanishes identically on the open interval $t<0$, so its first and second derivatives also vanish there. Substituting into the equation of motion immediately givesThus the force has no negative-time prehistory, consistent with the causal Green function of a damped oscillator. In fact, direct inversion or the positive-time substitution above identifies $f(t)=-e^{-2t}H(t)$; the given numerator $+i$ corresponds to a negative exponential force with this Fourier convention. Here $H$ is the Heaviside step function; its value exactly at the jump does not affect the transform.

The inverse Fourier transform isFor $t>0$, close in the upper half-plane, where $e^{i\omega t}$ decays. The single enclosed pole has residue $-i{e}^{-t}$, so counterclockwise integration gives $f=e^{-t}$. For $t<0$, close clockwise in the lower half-plane; there is no enclosed pole and the result is zero. ThusHere $H$ is the Heaviside step function. At the jump, symmetric Fourier inversion gives $f(0)=1/2$, the average of the one-sided limits. The transform itself cannot fix a function’s value at a single point, so this is the conventional representative rather than extra frequency information.

A direct check is ${\int }_0^{\infty }{e}^{-(1+i\omega )t}dt=(1+i\omega {)}^{-1}=-i/(\omega -i)$. Distributionally, $f^{\prime }=\delta (t)-e^{-t}H(t)$.

A causal Green function vanishes before the impulse. For $t>\tau$ it solves the homogeneous linear ordinary differential equation, whose characteristic roots are $-q\pm 2i$. Thus it has the formContinuity at $t=\tau$ is necessary to avoid a delta derivative in $G_{tt}$, so $A=0$. Integrating its forced equation across the impulse gives the jump condition for an impulse, $G_t({\tau }^{+},\tau )-G_t({\tau }^{-},\tau )=1$. The derivative from below is zero and the derivative from above is $2B$, givingHere $\Theta$ is the Heaviside step function. The zero initial data exclude any additional homogeneous response. By the superposition principle, the output for a general locally integrable input isThe upper endpoint causes no unwanted term in the first derivative because $G(t,t)=0$; in the second derivative its unit derivative jump reproduces the input.

The Heaviside step function restricts the integral defining the Fourier transform to $t\ge 0$, soThe convergence condition is essential. For real $\alpha$ it means $\alpha >0$, which is the decay assumption used in the subsequent parts. At $\alpha =0$ a distributional Fourier transform is needed; for negative real $\alpha$ the defining improper integral diverges.

For $T>0$ the difference of Heaviside step functions is an interval indicator, and $v=1/2$ on $(-T,T)$ and zero elsewhere. Its Fourier transform isThe apparent singularity at $\omega =0$ is removable. The formula remains valid for real negative $T$, because the original step-function difference and the oriented integral both reverse sign; for $T=0$ both vanish. This is the Fourier transform of an interval indicator and produces a scaled sinc function.

For the zero-data initial value problem, choose a causal Green function satisfying $L_{x}G(x,\xi )=\delta (x-\xi )$ and $G=0$ for $x<\xi$. Continuity at $x=\xi$ prevents a derivative of the Dirac delta distribution, while integration across the source givesFor $x>\xi$, part (i) gives the shifted homogeneous solutions. The zero value at the join removes the cosine term, and the unit derivative jump fixes the sine coefficient. Therefore the causal Green function of a damped oscillator isEquivalently, multiply the first expression by the Heaviside step function $H(x-\xi )$. The required solution isIts initial value and first derivative vanish: the interval has zero length at $x=0$ and the kernel vanishes at its upper limit. Differentiating twice shows that the upper-limit contribution to $y^{\prime \prime }$ is precisely $g(x)$, while the remaining homogeneous kernel terms cancel in $Ly$.