Measure temperature relative to ambient by ${\theta }_i=T_i-T_{\infty }$. Newton's law of cooling models ordinary cooling by ${\theta }_i^{\prime }=-a{\theta }_i$, since heat loss is proportional to the excess temperature. The additional negative forcing represents the cooling from the milk, with the temperature drop or rate normalized to one in the chosen units.

The Dirac delta function has unit integral and is supported at the instant of addition. It therefore imposes a sudden unit drop for cup 1. The difference $H(t-1)-H(t-2)$ of Heaviside step functions is one during the addition interval and zero outside it, giving a unit cooling rate for cup 2. The second temperature is continuous; only its derivative changes at the switching times. The value assigned to a Heaviside step function exactly at its jump does not change the integrated solution. This is the impulse versus finite-duration cooling comparison.

Taking the conjugate transpose of $U=A+iB$ and using that both parts are Hermitian matrices gives $U^{†}=A-iB$. Addition and subtraction recover the Hermitian parts of a unitary matrix:These depend only on $U$, so any such decomposition has the same $A,B$ and is unique. Conversely the displayed expressions are Hermitian matrices and sum to $U$ in the form $A+iB$, establishing existence as well as uniqueness. Their commutation and quadratic relation are proved in the two labelled parts.

First establish the interior convergence lemma for a power series. If $\sum a_n{w}^n$ converges at some $w\ne 0$, its terms are bounded: $|a_n{w}^n|\le M$ for a finite $M$. For every $|z|<|w|$,Comparison with a convergent geometric series proves absolute convergence at $z$.

Now letThe set contains zero, since at $w=0$ only the constant term remains. If $|z|<R$, the definition of supremum gives a point of convergence $w$ with $|w|>|z|$. The lemma proves absolute convergence at $z$. If $|z|>R$, convergence would place $|z|$ in the set defining $R$, a contradiction. This proves the existence of a radius of convergence, including $R=0$ and $R=\infty$:The proof leaves behavior on $|z|=R$ open; it need not be uniform around that circle. The two examples illustrate opposite possibilities.

The Absolute Frobenius morphism is the identity on the underlying topological space and acts on the structure sheaf by $a\mapsto a^p$. On an affine scheme it is induced by the ring homomorphism $A\to A$, $a\mapsto a^p$. This is additive in characteristic $p$, and the inverse image of a prime ideal $\mathfrak{p}$ is again $\mathfrak{p}$, since $a^p\in \mathfrak{p}$ exactly when $a\in \mathfrak{p}$. These affine maps agree under localization and hence glue.

A perfect ring of characteristic p is one whose Frobenius endomorphism is bijective. Since the Absolute Frobenius morphism fixes points, its map on the stalk at $x$ is precisely the Frobenius endomorphism of ${\mathcal{O}}_{X,x}$. A morphism of sheaves is an isomorphism exactly when it is an isomorphism on all stalks. ThereforeIn the reverse direction the inverse sheaf map, together with the identity on points, gives the inverse morphism of locally ringed spaces. This is perfection in the Frobenius sense, not the unrelated homological use of the word perfect.

For the function field assertions below, interpret the normal schemes as integral, as is implicit in writing a single field $k(X)$ or $k(Y)$. General normal schemes need not be irreducible; without this interpretation those fields are not defined globally as single fields. Set $K=k(Y)\subseteq L=k(X)$. For an affine open $V=\mathrm{Spec}A$ in $Y$, finiteness gives $f^{-1}V=\mathrm{Spec}B$, with $A\subseteq B$ finite and both rings normal domains. In fact $B$ is the integral closure of $A$ in $L$: any element of $L$ integral over $A$ is also integral over $B$ and therefore belongs to $B$. This local description yields the two forms of Frobenius factorization through a finite normal cover.

Let $f$ be probability density per unit solid angle, with ${\int }_{S^2}fd\Omega =1$ and isotropic initial condition $f(\mathbf{p},0)=1/(4\pi )$. In still fluid, bottom-heavy spherical-cell orientation dynamics gives $\dot{\mathbf{p}}=B^{-1}[\mathbf{k}-(\mathbf{k}\cdot \mathbf{p})\mathbf{p}]$. Adding isotropic rotational diffusion gives the gyrotactic orientation Fokker-Planck equationThe derivatives are tangential to the unit sphere. Rotational symmetry about $\mathbf{k}$ is preserved by both the equation and the initial condition, so $f$ depends only on $\theta ,t$. Since $\dot{\theta }=-\sin \theta /B$, the polar probability flux is $j_{\theta }=-\sin \theta f/B-D_r{f}_{\theta }$. Expanding ${\partial }_{t}f+{\sin }^{-1}\theta {\partial }_{\theta }(\sin \theta j_{\theta })=0$ yieldsRegular solutions have no probability flux through the poles. The mean swimming velocity is ${\mathbf{V}}_c=V_s⟨\mathbf{p}⟩$, with only a vertical component by symmetry.

The linear damped harmonic oscillator has characteristic exponents $-ϵ\pm i\sqrt{1-{ϵ}^2}$. Imposing both initial conditions gives, for $|ϵ|<1$,For completeness, at $|ϵ|=1$ its continuous limiting form is $t{e}^{-ϵt}$; for $|ϵ|>1$ replace $\sin (\omega t)/\omega$ by $\sinh (\sqrt{{ϵ}^2-1}t)/\sqrt{{ϵ}^2-1}$. The initial derivative is one, as the original PDF specifies; the converted TeX loses this derivative and corrupts the time derivatives in both oscillator equations.

The separate nonlinear oscillator has quadratic drag. Introduce the slow time $T=ϵt$ and use the method of multiple scales,The order-$ϵ$ equation isFor positive amplitude, write $\vartheta =t+\varphi$. Its right side is $-2A_T\cos \vartheta +2A{\varphi }_T\sin \vartheta -A^2|\cos \vartheta |\cos \vartheta$. The cosine coefficient of the last angular function iswhile its sine coefficient is zero. The solvability condition in the method of multiple scales removes the resonant harmonics, givingThe leading initial conditions select $A(0)=1$, $\varphi (0)=0$. Therefore the quadratically damped oscillator hasThe same sign follows from the exact energy balance $d[({\dot{x}}^2+x^2)/2]/dt=-ϵ|\dot{x}{|}^3$. Averaging $|\cos \vartheta {|}^3$ gives $4/(3\pi )$, consistent with the amplitude equation.

For $ϵ>0$, the amplitude decays algebraically and stays at most one. On each fixed interval $0\le ϵt\le T_{\ast }$, the omitted displacement correction is $O(ϵ)$, so this approximation is valid through times of order ${ϵ}^{-1}$. The effective drag strength $ϵA$ decreases rather than increases, and the envelope predicts $A\sim 3\pi /(4ϵt)$ at later times. There is no finite-time pole in the positive-drag approximation. A uniform error assertion for arbitrarily long intervals is stronger than the finite-slow-time expansion just derived; the stated $O(ϵ)$ estimate applies to the specified $ϵt=O(1)$ range.

For $ϵ<0$, write $\delta =1-4|ϵ|t/(3\pi )$. The amplitude grows as $1/\delta$, with a formal pole atOn fixed slow-time intervals strictly before that pole, the displacement error remains $O(|ϵ|)$. This estimate is not uniform as $\delta \to 0$: nonresonant corrections have size comparable to $|ϵ|A^2$, and the required weak-drag condition is $|ϵ|A=|ϵ|/\delta \ll 1$. Thus the approximation necessarily fails when $\delta =O(|ϵ|)$, corresponding to a distance of order one in fast time from the formal pole, and may require additional error control when approaching it. The plot therefore stops short of the pole. Its position is a prediction of the leading amplitude equation, not an exact blow-up time for the original nonlinear equation; continuing the formula past it is invalid.

For a slender viscous bubble, approximate the interface normal by the radial direction and its curvature by $1/a$. Axial curvature and axial derivatives of the leading external motion are smaller by the slenderness factors. The exterior leading Stokes flow is locally radial. Incompressibility gives ${\partial }_r(r{u}_r)=0$, while the kinematic boundary condition gives $u_r(a)=a_t$ to leading order, henceThe radial viscous Laplacian of this $1/r$ velocity vanishes, so exterior pressure is constant to leading order; set its reference value to zero. Its radial normal stress at the interface is $2\mu {\partial }_r{u}_r{|}_a=-2\mu a_t/a$. Inner viscous normal stress is smaller by the viscosity ratio $\lambda$, and the bubble's leading normal stress is $-P$. The interfacial stress balance with variable surface tension, here with constant $\gamma$, consequently givesThese are local leading-order relations: the outer axial flow and tangential traction enter at higher slenderness order. The sign also has a useful check: a uniform cylinder with no imposed excess pressure contracts, $a_t=-\gamma /(2\mu )$, under surface tension.

For a fixed fluid domain with prescribed velocity on its whole boundary, the Minimum-dissipation theorem for Stokes flow compares the actual incompressible Stokes flow $u$ with every admissible incompressible velocity field $v$ having the same boundary values. The trial fields need not themselves satisfy the momentum equation. With $e(u)=(\nabla u+\nabla u^T)/2$ and constant dynamic viscosity $\mu$, define the viscous dissipationAssume no nonconservative body force, or absorb a conservative force into the pressure. Set $w=v-u$, so $\nabla \cdot w=0$ and $w=0$ on the boundary. The Stokes stress tensor $\sigma =-pI+2\mu e(u)$ obeys $\nabla \cdot \sigma =0$. Because pressure has zero contraction with $e(w)$,Expanding the square provesEquality requires $e(w)=0$, so $w$ is a rigid motion; the homogeneous boundary values force it to vanish. The theorem concerns fixed boundary velocities, not fixed applied forces or arbitrary comparisons between different fluid domains.

For the concentric-sphere motion, seek $u=F(r)\mathbf{\Omega }\times x$. Incompressibility is automatic. The choice $F=A+B/r^3$ gives a harmonic velocity: the linear term is harmonic, and each component of $x/r^3$ is a derivative of the harmonic function $1/r$. The no-slip conditions $F(a)=1$, $F(b)=0$ giveThus the Rotational Stokes flow between concentric spheres isIt satisfies the Stokes equation with constant pressure, and the uniqueness implied by the theorem identifies it as the required flow.

In the Unscaled Papkovich–Neuber representation, take $\Phi =-u/2$ and $\chi =0$. Both potentials are harmonic and $x\cdot \Phi =0$, so the representation gives $u=-2\Phi$. Alsoand therefore $p=0$, after choosing the arbitrary pressure constant.

This pressure conclusion also follows from symmetry and Linearity of Stokes flow. Pressure is a scalar linear in the axial vector $\mathbf{\Omega }$, and isotropy leaves only a possible term proportional to $\mathbf{\Omega }\cdot x$. That is a pseudoscalar, forbidden by reflection symmetry of the concentric spherical domain. Hence the pressure perturbation vanishes. Equivalently, axisymmetry and reflection with flow reversal make the motion purely azimuthal, while its meridional momentum balance has no pressure gradient.

Let $v=\mathbf{\Omega }\times x$ in the following stress calculation. The symmetric velocity gradient isThe full stress tensor, in the zero-pressure gauge, is thereforeFor $n=x/r$, the traction on a radial surface is $\sigma n=-3\mu B(\mathbf{\Omega }\times n)/r^3$. This completely specifies the stress field; with the polar axis along $\mathbf{\Omega }$, its only nonzero spherical shear component is ${\sigma }_{r\varphi }=-3\mu B\Omega \sin \theta /r^3$.

The fluid couple on the inner sphere is the integral of $x\times \sigma n$ with the normal pointing from the solid into the fluid. Using $\int n_i{n}_{j}d\omega =(4\pi /3){\delta }_{ij}$,Hence the sustaining applied Torque in rotational Stokes flow between concentric spheres isFor $a\ll b$, it approaches the isolated rotating sphere in Stokes flow result $8\pi \mu a^3\mathbf{\Omega }$. For a narrow gap $d=b-a\ll a$, it becomes $8\pi \mu a^4\mathbf{\Omega }/(3d)$: the divergent resistance is that of a thin Couette shear layer, with speed scale $a\Omega$ and gradient scale $a\Omega /d$.

Now add the force-free, couple-free rigid particles and denote the resulting velocity by $u_p$. Extend $u_p$ through each particle as its rigid translational and rotational velocity. No slip makes this extension continuous across every particle surface; it is incompressible and has zero rate-of-strain tensor inside each particle. It therefore gives an admissible trial velocity in the original particle-free annulus, with exactly the same sphere boundary velocities. Its full-annulus dissipation is the actual fluid dissipation ${\mathcal{D}}_p$, since the added interior integrals are zero. The Minimum-dissipation theorem for Stokes flow consequently gives ${\mathcal{D}}_p\ge {\mathcal{D}}_0$.

The boundary-work identity follows by integrating $\nabla \cdot (\sigma u)$: dissipation equals the mechanical power supplied at all rigid boundaries. The outer sphere is stationary; the inner sphere has no translation, so any force needed to hold its centre does no work; and each added particle contributes zero power because its net force and couple are zero. Thus ${\mathcal{D}}_p={\mathbf{G}}_p\cdot \mathbf{\Omega }$, while ${\mathcal{D}}_0={\mathbf{G}}_0\cdot \mathbf{\Omega }$. ThereforeFor $\mathbf{\Omega }\ne 0$ and at least one particle of positive volume in the open annulus, the inequality is strict. Equality would make the extended velocity identical to the particle-free solution, yet that solution has nonzero strain on every open region away from the rotation axis and cannot be rigid throughout any particle interior. Hence the applied couple component along the imposed angular velocity increases. This is the fixed-velocity extra dissipation due to a rigid inclusion argument; it does not assert an increase in every transverse torque component.

In the Boussinesq approximation, density departures from a constant reference density ${\rho }_0$ are small enough to neglect in inertia and mass continuity, but their gravitational effect is retained. Write the total buoyancy acceleration as $b=-g(\rho -{\rho }_0)/{\rho }_0$, choose a hydrostatic background $\bar{b}(z)$, and set $\sigma =b-\bar{b}$. Thus $\sigma$ is the buoyancy perturbation, positive for a parcel lighter than the background at its present height. The squared buoyancy frequency isStable stratification has $N^2>0$. A small vertical displacement $\xi$ of a parcel conserving its buoyancy gives $\sigma =-N^2\xi$, explaining the restoring acceleration and the name of $N$.

Let $P=p^{\prime }/{\rho }_0$ be pressure with the background hydrostatic pressure removed, $\mathbf{u}=(u,v,w)$, and $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ the material derivative. The ideal, nonrotating Boussinesq equations areThe last equation is equivalent to $Db/Dt=0$. There is neither viscosity nor buoyancy diffusion. The reference density occurs in the pressure normalization, while the small density anomaly occurs only through buoyancy.

For the two-dimensional disturbance, linearization about rest givesA plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$ therefore satisfiesThese internal-wave polarization relations and the dispersion relation determine the individual responses in parts i and ii. A radiating response has no incoming internal gravity wave from $x=+\infty$; an evanescent response is bounded there. The harmonic solutions describe the response after excluding arbitrary initial transients.

For the localized oscillatory boundary, define $\alpha =\omega /\sqrt{N^2-{\omega }^2}$ and $I(z)=ϵ{\int }_0^{\infty }{\hat{u}}_0(m)e^{imz}dm=f(z)+ig(z)$. Both signs of vertical wavenumber must radiate into $x>0$. Their horizontal wavenumber is $k=\alpha |m|$, so the positive-$m$ contribution involves $Z_1=z+\alpha x$ and the negative-$m$ contribution involves $Z_2=z-\alpha x$. Consequently the complex horizontal velocity isTaking its real part provesThe positive-$m$ internal-wave polarization has $w=-\alpha u$ and the negative-$m$ one has $w=+\alpha u$. HenceIn particular, $2f(z)=ϵu_0(z)$, so the boundary condition is recovered exactly. This is two-beam radiation from a localized oscillating boundary.

Take a smooth localized profile, for example $u_0(z)=e^{-z^2/H^2}$. Its effective support is a neighbourhood of width a few $H$; increasing that neighbourhood makes the tails arbitrarily small. The boundary displacement is $ϵu_0(z)\sin \omega t/\omega$ to first order, so at $t=0$ it is flat and moving at maximum speed. At that instant,There are two narrow oblique beams centered on $z=-\alpha x$ and $z=+\alpha x$. Away from their overlap, the velocity is parallel to the relevant beam, in the directions $(1,-\alpha )$ and $(1,+\alpha )$. The figure uses $\alpha =3/4$ and measures velocities in units of $ϵ$.

At other phases the $g$ terms contribute. For the positive-frequency convention above, $g$ is the Hilbert transform of $f$; it generally has longer tails than $f$. Thus even when the boundary is nearly stationary, the interior need not be stationary or confined to the bright cores shown here. Half a period later all velocities reverse. The beam directions remain the same, but the instantaneous velocity profile and its quadrature component change.

A short pulse has a broad frequency spectrum. Frequencies below $N$ produce outgoing internal gravity wave packets with frequency-dependent ray slopes $\pm \omega /\sqrt{N^2-{\omega }^2}$; frequencies above $N$ contribute an evanescent near-boundary response. Moreover the group velocity depends on wavelength as well as direction. The radiated packet therefore spreads into a dispersive fan or wake, with different frequencies traveling at different angles and different scales separating along those directions. After the forcing ends the propagating packet travels away, rather than maintaining the two permanently driven monochromatic beams.

Assume global four-dimensional N=1 supersymmetry with canonical kinetic terms and no explicit gaugino mass. Write $W_i={\partial }_{i}W$ and $W_{ij}={\partial }_i{\partial }_{j}W$. To match the displayed stationarity equation, use $F^i=\overline{W_i}$; this differs by a minus sign from the eliminated auxiliary field in the component convention of Solution 2. The real D-terms are $D^a=-g^a({\bar{\varphi }}_i(T^a{)}^i{}_j{\varphi }^j+{\xi }^a)$, with Hermitian gauge generators. A constant Fayet–Iliopoulos term is permitted only for an appropriate Abelian factor. The scalar potential isDifferentiate with respect to ${\varphi }^i$, holding $\bar{\varphi }$ independent:Thus stationary point conditions gives the first requested identity. Gauge invariance of the superpotential gives $W_i(T^a{)}^i{}_j{\varphi }^j=0$. Conjugating it and using Hermiticity gives ${\bar{\varphi }}_j(T^a{)}^j{}_i{F}^i=0$, the second condition in the form needed below.

At the vacuum, define $H_{ij}=⟨W_{ij}⟩$ and $G_{ia}=g^a⟨{\bar{\varphi }}_j⟩(T^a{)}^j{}_i$. The two conditions combine asHere the upper equation is stationarity and the lower equation is conjugated gauge invariance. The chiral-superfield fermion mass matrix contributes $-\frac{1}{2}{H}_{ij}{\psi }^i{\psi }^j$, and the gauge Yukawa interaction contributes the mixed chiral-fermion/gaugino bilinear. With canonical Weyl fermion and gaugino kinetic terms, and choosing the gaugino phase to put the mixed bilinear in real-sign convention, these terms areThe factor $\sqrt{2}$ is important: the matrix directly acting on the printed $(F,D)$ is the same mass bilinear with rescaled gauginos, rather than literally the canonical matrix. In detail, for $S=\mathrm{diag}(I,\sqrt{2}I)$,Since supersymmetry breaking makes this vector nonzero, there is at least one zero mass eigenvalue. A normalized massless fermion isWith the consistent gaugino phase, the inhomogeneous vacuum supersymmetry transformation of $(\psi ,\lambda )$ is proportional to the same vector $(F,D/\sqrt{2})$. Therefore $\eta$ shifts by a nonzero constant times the transformation parameter: it is the Goldstino, the Goldstone fermion of broken global supersymmetry. The argument guarantees at least one Goldstino zero mode, not the absence of other massless fermions.

For the second calculation, first take only canonical chiral superfields and F-term breaking. At a stationary vacuum $v$, let $h=\varphi -v$, $m_{ij}=W_{ij}(v)$ andThese follow by differentiating $V={\sum }_k|W_k{|}^2$ twice. The quadratic scalar potential isThe doubled scalar mass matrix has blocksTransforming to the $2N$ canonically normalized real components of $h$ preserves its trace, so the sum of real scalar mass squares is $\mathrm{tr}{M}_s^2=2\mathrm{tr}A$. The Weyl fermion mass squares are the eigenvalues of $m^{†}m=A$ and their sum is $\mathrm{tr}A$. Each Weyl mass eigenstate has two physical spin states. The printed supertrace has scalar sign $-1$ and fermion weight $+2$, henceThe frequently used boson-minus-fermion definition has the opposite overall sign and gives the same zero. The mixing block $B$ can split scalar masses but cannot change this trace. This proves the tree-level supertrace mass sum rule; it uses canonical kinetic terms and a renormalizable global theory, rather than a general supergravity or noncanonical effective theory.

In the MSSM, trying to make all squarks and sleptons heavy through direct canonical tree-level F-term breaking encounters this constraint. In an unbroken electric-charge and color sector, the scalar mass-squared trace is twice the trace for its ordinary fermionic partners. More explicitly, take a unit vector $u$ minimizing $u^{†}Au$. For $h=e^{i\alpha }u$, choose $\alpha$ so that $\mathrm{Re}(e^{2i\alpha }{u}^{T}Bu)\le 0$. The curvature along this real scalar direction is then at most $u^{†}Au$, so at least one scalar mass square is no larger than the smallest fermion mass square in that sector. For one isolated pair the scalar masses are $m_f^2\pm |B|$. Thus light quarks can force an unacceptably light squark, or a tachyon, under these assumptions. A hidden supersymmetry-breaking sector with mediation produces effective soft supersymmetry breaking; radiative effects, noncanonical kinetic terms, or supergravity can change the assumptions and allow a viable spectrum. The trace over an arbitrary enlarged theory alone would not locate the light state in the visible sector; the MSSM conclusion uses the conserved-charge block structure.