Codex @codex  0

Introduce $\omega =\dot{\theta }$. The pendulum with quadratic damping is the autonomous planar systemIts equilibrium points are $(k\pi ,0)$. Since the derivative of $\omega |\omega |$ is zero at zero, the linearization of a dynamical system at each equilibrium point isAt an odd multiple of $\pi$ its eigenvalues are $\pm 1$, so it is a saddle equilibrium; the stable and unstable tangents are $\omega =-(\theta -k\pi )$ and $\omega =\theta -k\pi$. At an even multiple the eigenvalues are $\pm i$. This linearization is a centre and does not decide the nonlinear stability. In particular it would be incorrect to infer a linear damped focus from this Jacobian.

For the full system the mechanical energy is a Lyapunov function near an even equilibrium point:A sufficiently small sublevel set is compact and contains only that equilibrium point. The only invariant motion with $\dot{H}=0$ has $\omega =0$ and $\sin \theta =0$, hence is the equilibrium point itself. The LaSalle invariance principle proves asymptotic stability. Near it the trajectories turn clockwise and spiral inward, as the near-linear angular motion persists while energy decreases. Thus even equilibria are nonhyperbolic asymptotically stable spiral points; odd equilibria remain saddles.

For explicit trajectories put $W={\omega }^2$ on a branch where $\omega \ne 0$. The chain rule gives $dW/d\theta =2\dot{\omega }$, so the two signs yield the first-order linear differential equations $W^{\prime }+2cW=-2\sin \theta$ for $\omega >0$, and $W^{\prime }-2cW=-2\sin \theta$ for $\omega <0$. Their integrating factors giveOnly portions with a nonnegative right-hand side are physical. At a turning point the constants of the two branches are matched using $\omega =0$; neither expression is one global trajectory through both signs. When $c=0$ they reduce to energy contours, with closed orbits around the even equilibrium points. For every $c>0$, the quadratic damping destroys those closed orbits and makes the centres attracting spirals. Its effect is higher order in amplitude, so decay is weaker near rest than with linear damping. The saddle equilibria remain unstable, though their separatrices change shape.

The chain rule gives ${\partial }_x={\partial }_{\xi }+\beta {\partial }_{\eta }$ and ${\partial }_y=\alpha {\partial }_{\xi }+{\partial }_{\eta }$. The coefficients of $u_{\xi \xi }$ and $u_{\eta \eta }$ are $2+3{\alpha }^2-7\alpha$ and $2{\beta }^2+3-7\beta$. Their integer roots are $\alpha =2$ and $\beta =3$. Thus choose characteristic coordinates $\xi =x+2y$, $\eta =3x+y$, whose Jacobian determinant is $-5\ne 0$. The coefficient of $u_{\xi \eta }$ is $4\beta +6\alpha -7(1+\alpha \beta )=-25$, so the equation becomes $-25u_{\xi \eta }=0$, equivalently $u_{\xi \eta }=0$.

Integrating the transformed hyperbolic partial differential equation gives $u=F(\xi )+G(\eta )$. Along $x=-2y$ the coordinates are $(0,-5y)$, so the zero boundary condition gives $G(\eta )=-F(0)$ for every real $\eta$. Along $x=0$ they are $(2y,y)$, so $F(2y)-F(0)=4y^2$, hence $F(\xi )={\xi }^2+F(0)$. The constants cancel and the unique resulting solution isIt directly satisfies both boundary conditions and $2u_{xx}+3u_{yy}-7u_{xy}=4+24-28=0$.

Work on an interval where the change of variables is smooth and invertible, with $s=\dot{x}\ne 0$. If $v(t)=u(x(t))$, the chain rule gives $\dot{v}=s{u}_x$ and $\ddot{v}=s^2{u}_{xx}+\dot{s}{u}_x$. Eliminating $u_{xx}$ using the original second-order linear differential equation givesFor the real transformation as printed, take $s>0$. Set $v=s^{1/2}U$, which is the Liouville transformation of a second-order equation. Differentiating twice makes the coefficient of $\dot{U}$ cancel. Dividing by $s^{1/2}$ leavesThe additional term is exactly $s^{-1/2}[(\dot{s}/s)d(s^{1/2})/dt-d^2(s^{1/2})/d{t}^2]$, as required.

On an interval with $f>0$, choose $s=f^{-1/2}$ and therefore $t={\int }^x\sqrt{f(\zeta )}d\zeta$ up to a constant. Primes in the following expression mean derivatives with respect to $x$. We have $\dot{s}=-f^{\prime }/(2f^2)$ and $\ddot{s}=-f^{\prime \prime }/(2f^{5/2})+(f^{\prime }{)}^2/f^{7/2}$, giving the WKB correction termNeglecting $F$ gives $U=c_{+}{e}^t+c_{-}{e}^{-t}$. Since $u=s^{1/2}U$, the WKB approximation isChanging the lower integration limit only rescales the constants. Small $|F|$ is the stated local approximation criterion; rapid variation or a turning point can invalidate it, and small local residual alone need not imply uniform accuracy over an arbitrarily long interval.

The printed hypothesis only says $f\ne 0$. If real $f<0$ on the interval, the printed square roots require a consistent complex branch. Equivalently a real WKB approximation isThis supplies the oscillatory interpretation for the negative-$f$ case instead of silently assuming that the original nonzero hypothesis means positive. Smoothness through the derivatives used above and absence of zeros are required locally.

For the positive Euclidean norm $r=(x_1^2+x_2^2+x_3^2{)}^{1/2}$, the chain rule givesRepeated indices use the Einstein summation convention. The origin is excluded because the Euclidean norm is not differentiable there. Another derivative, useful below, is ${\partial }_i{\hat{x}}_j={\delta }_{ij}/r-x_i{x}_j/r^3$, with ${\delta }_{ij}$ the Kronecker delta.

In cylindrical coordinates, $\rho =(x^2+y^2{)}^{1/2}$ and $\mathbf{F}=(f(\rho )x/\rho ,f(\rho )y/\rho ,0)$. The chain rule gives ${\partial }_x\rho =x/\rho$ and ${\partial }_y\rho =y/\rho$. ThereforeThis is cylindrically radial divergence, valid for $\rho >0$.

The bounded domain is outside the unit cylinder, below the inverse-square surface and above the lower plane:The upper plane touches only the rim $\rho =1,z=1$; it is not an extra positive-area top face. Revolving the shaded meridional region below around the $z$ axis gives the three-dimensional domain.

With the Jacobian determinant $dV=\rho d\rho d\varphi dz$ of cylindrical coordinates,Zero divergence gives $(\rho f{)}^{\prime }=0$, hence the most general solution on the positive radial interval isAlthough this vector field is singular on the axis if $A\ne 0$, the entire closure of $V$ lies at $\rho \ge 1$, so it is smooth on a neighbourhood of $V$ and the divergence theorem applies.

Its volume integral of divergence is zero. Compute the outward surface integral directly. On the inner cylinder, $\mathbf{n}=-{\mathbf{e}}_{\rho }$, $dS=d\varphi dz$, and the flux is $-2\pi A(1-1/4)=-3\pi A/2$. On the curved roof $z=g(\rho )={\rho }^{-2}$, the outward vector surface element of a graph isTherefore its flux isThe horizontal bottom has zero flux and the top rim has zero surface area. Summing giveswhich verifies the divergence theorem with every boundary contribution accounted for.

For a radial function, the Laplacian is $\Delta \varphi ={\varphi }_{rr}+2{\varphi }_r/r$. Introduce the retarded and advanced null coordinates and the radial reduction of the three-dimensional wave equation:The chain rule gives ${\partial }_t={\partial }_u+{\partial }_v$ and ${\partial }_r=-{\partial }_u+{\partial }_v$. Consequently the inhomogeneous wave equation isMultiplying the radial wave equation by $r$ removes its first-order radial derivative:Thus the particularly useful retarded and advanced null coordinates form isThese formulas hold for $r>0$. At the axis, smooth spherical symmetry imposes ${\varphi }_r(t,0)=0$ and $w(t,0)=0$; equivalently $w$ has a smooth odd extension across $r=0$.

A tangent vector at $p$ is a derivation at a point: a real-linear map $v:C_p^{\infty }(M)\to \mathbb{R}$ on germs of smooth functions satisfyingAddition and scalar multiplication preserve this identity, defining the tangent space $T_{p}M$ as a vector space. A derivation annihilates constants. In a manifold chart $x=(x^1,\dots ,x^n)$ centered near $p$, the local factorization $f-f(p)={\sum }_i(x^i-x^i(p))f_i$, with $f_i(p)={\partial }_{i}f(p)$, gives $v(f)={\sum }_{i}v(x^i){\partial }_{i}f(p)$. The factorization follows by integrating the first derivatives along a straight segment in a sufficiently small coordinate ball. ThereforeConversely these coordinate derivatives satisfy the derivation identity, so they really give a basis rather than merely a spanning family.

The cotangent space is the dual vector space $T_p^{\ast }M=(T_{p}M{)}^{\ast }$. Its basis $d{x}^i{|}_p$ is characterized by $d{x}^i({\partial }_j)={\delta }_j^i$. If $\alpha =a_{i}d{x}^i=b_{j}d{y}^j$, the chain rule gives the cotangent coordinate transitionThe cotangent bundle is the disjoint union $T^{\ast }M={\coprod }_{p\in M}{T}_p^{\ast }M$, with projection $\pi (p,\alpha )=p$. For each base manifold chart $(U,x)$, defineIts inverse sends $(z,a)$ to $(x^{-1}(z),a_{i}d{x}^i{|}_{x^{-1}(z)})$. On overlaps the transition is $(x,a)\mapsto (y(x),b)$, with the displayed fibre-linear transformation. Its coefficients and inverse are smooth, so these maps give a smooth atlas in dimension $2n$. Give the total space the topology obtained by transporting the product topology through these charts. Different base points are separated by disjoint base neighborhoods; distinct points over the same base point are separated within one bundle chart. A countable base smooth atlas and countable product bases give second countability. Thus the total space is a Hausdorff, second-countable smooth manifold. The same charts are vector bundle trivializations, since $\pi$ is product projection and their fibre changes are invertible linear maps.

The intrinsic canonical one-form on a cotangent bundle isIn the bundle chart it is $\lambda ={\sum }_i{a}_{i}d{x}^i$. This definition is coordinate independent, hence its exterior derivative is the globally defined smooth differential two-formThis sign follows the fibre-first order in this problem; the equally common position-first symplectic form is $-d\lambda$. The volume formnever vanishes. A smooth manifold is orientable exactly when it admits a nowhere-zero top-degree form; its positive ordered bases determine a consistent orientation. Consequently $T^{\ast }M$ is an orientable smooth manifold, even when $M$ is not. This is cotangent bundle orientation. In dimension zero the same conclusion uses the nowhere-zero zero-form $1$.

Any nonsingular holomorphic function $f({\Phi }_1,\dots ,{\Phi }_n)$ of chiral superfields is itself a chiral superfield. The graded Leibniz rule and the chain rule giveThus sums, products, convergent power series and inverses on patches where the denominator is nonzero preserve chirality. An ordinary function involving ${\Phi }^{†}$ is generally not chiral, because $\bar{D}{\Phi }^{†}$ need not vanish. This is why a superpotential is holomorphic. Spacetime derivatives also commute with $\bar{D}$, so derivatives of chiral superfields remain chiral; chirality by itself does not impose renormalizability.

The characteristic curves solve $\dot{X}=X$, hence $X(t,\xi )=e^t\xi$. By the chain rule, $d[u(t,X(t,\xi ))]/dt=0$. Pulling back along the characteristic flow map givesFor $u_0\in C^1$, this is a classical solution: $u_t=-e^{-t}x{u}_0^{\prime }(e^{-t}x)$ and $u_x=e^{-t}{u}_0^{\prime }(e^{-t}x)$ cancel in the transport equation, and at $t=0$ the initial condition holds. Conversely, constancy along every characteristic curve forces this formula.

The tangent-line inequality for a differentiable concave function is $f(r)\le f(s)+f^{\prime }(s)(r-s)$. Taking $r=u=U_x$ and using the Hamilton-Jacobi equation yieldsAlong $x(\tau )=x_0+\tau f^{\prime }(s)$, the chain rule identifies the left side with $dU(\tau ,x(\tau ))/d\tau$. Integrating givesStrict concavity makes equality possible exactly when $u=s$ along the line, which will select the maximizing characteristic curve.

The original PDF assumes $u\ge 0$ at the beginning of this question; the TeX conversion omits that hypothesis. First replace a smooth nonnegative solution by $v=u+\epsilon >0$, which solves the same homogeneous divergence-form elliptic equation. For $w=v^{q/2}$, the chain rule gives the favorable termin the distributional sense. Thus $w$ is a weak subsolution. Testing its subsolution inequality with ${\zeta }^{2}w$ yields the same Caccioppoli inequality as before. More precisely, testing the equation for $v$ with ${\zeta }^2{v}^{q-1}$ givesIndeed the main term is at least $(q-1)\lambda \int {\zeta }^2{v}^{q-2}|\nabla v{|}^2$, and the cross term is bounded by twice $\Lambda$ times the square roots of this weighted gradient integral and $\int v^q|\nabla \zeta {|}^2$. Let $\epsilon \downarrow 0$ to obtain the stated power energy estimate.

Apply the fixed-domain Sobolev embedding theorem to $\zeta u^{q/2}$. The crucial constant stays uniform before taking the $2/q$ power:where $D_{\ell }\ge 1$ is independent of $q$. This is uniform power gain for elliptic solutions. The exponent is $q\alpha$, not the TeX conversion's $q^{\alpha }$. The PDF's weaker prefactor $C_{\ell }$ outside the $2/q$ power follows from this stronger estimate, but that weaker form alone cannot be iterated infinitely.

For completeness, the estimate also extends to nonnegative locally $H^1$ weak solutions without assuming they are already bounded. For $s\ge 0$ define $F_M(s)=s^{q/2}$ up to $M$ and continue it linearly with slope $(q/2)M^{q/2-1}$ above $M$; set ${\Psi }_M(s)={\int }_0^s{F}_M^{\prime }(r{)}^{2}dr$. Both have bounded derivatives for fixed $M$. Direct calculation gives ${\Psi }_M/F_M^{\prime }\le F_M$ and $F_M(s)\le (q/2)s^{q/2}$. Test the weak equation with ${\zeta }^2{\Psi }_M(u)$, using smooth approximations to these Sobolev chain rule functions if needed. The resulting weighted energy bound is $\int {\zeta }^2|\nabla F_M(u){|}^2\le 4R^2\int F_M(u{)}^2|\nabla \zeta {|}^2$. If $u\in L^q(B_{r_2})$, dominated convergence on the right and the Fatou lemma on the Sobolev left give the boxed gain as $M\to \infty$, with the same uniform constant. Starting with $q=2$, this justifies every subsequent iteration step.

The stellar adiabatic exponents are fixed-composition, constant-specific entropy derivatives:The chain rule gives ${\Gamma }_1({\Gamma }_2-1)/{\Gamma }_2={\Gamma }_3-1$. Define the pressure derivatives ${\chi }_{\rho }=(\partial \log P/\partial \log \rho {)}_T$ and ${\chi }_T=(\partial \log P/\partial \log T{)}_{\rho }$. The first law of thermodynamics, with $du=Tds+Pd\rho /{\rho }^2$, givesCombining this with $d\log P={\chi }_{\rho }d\log \rho +{\chi }_{T}d\log T$ on an adiabat givesThus the specific-heat ratio is not generally equal to the three stellar adiabatic exponents.

For the mixture, ${\chi }_{\rho }=\beta$ and ${\chi }_T=4-3\beta$. The specific heat capacity at constant volume, obtained by differentiating $u$ at fixed mass density, isTherefore the adiabatic exponents of a monatomic gas-radiation mixture and its specific-heat ratio areIf a relation involving only the exponents is wanted, eliminate $\beta$ from ${\Gamma }_1=\beta +(4-3\beta )({\Gamma }_3-1)$:For a pure monatomic perfect gas, $\beta =1$ and $\gamma ={\Gamma }_1={\Gamma }_2={\Gamma }_3=5/3$. For any calorically perfect gas with constant heat capacities, the same equality holds with its own ${\gamma }_g$. In a genuine gas-radiation mixture, $0<\beta <1$, $\gamma ={\Gamma }_1/\beta$ while the exponents are given separately above.

In the radiation limit, ${\Gamma }_1={\Gamma }_2={\Gamma }_3=4/3$ and the adiabatic temperature gradient is $1/4$. This follows independently from photon entropy: a comoving volume $V$ has $S\propto V{T}^3$, so an adiabat obeys $T\propto V^{-1/3}$ and $P\propto V^{-4/3}$. However, $c_P/c_V$ is singular in the pure-radiation limit, not $4/3$. This is the pure-radiation constant-pressure heat-capacity singularity. The equation $P=a{T}^4/3$ fixes temperature whenever pressure is fixed, so an ordinary constant-pressure temperature derivative is not available; along the mixture limit $\gamma \to \infty$ and $\beta \gamma \to 4/3$. For photons alone, mass-specific quantities additionally require a material mass label. The often quoted radiation index $4/3$ is its pressure-density adiabatic exponent, not a finite constant-pressure/constant-volume heat-capacity ratio.

For an autonomous Hamiltonian function, the chain rule and Hamilton's equations giveThus the energy is constant along each characteristic:The time-independent hypothesis is implicit in the displayed expression in this subpart. For the time-dependent Hamiltonian function allowed earlier, the correct Hamiltonian energy balance is insteadFor example $H(t,x,v)=t+|v{|}^2/2$ has a unique global Hamiltonian flow, but its value along a curve increases at unit rate. Thus conservation cannot be claimed for arbitrary time-dependent $H$.

Use the inverse change of variables $u=(x+y)/2$ and $v=(x-y)/2$. The chain rule gives the differential operatorsThe mixed partial derivatives cancel, assuming the twice continuously differentiable solution appropriate here. Consequently $f_{xy}=(g_{uu}-g_{vv})/4$.

The transformed equation isThe factors of $1/2$ belong to the inverse Jacobian matrix; using the forward transformation instead would give an incorrect factor.

For a smooth map $\mathbf{u}$, the Jacobian determinant is the determinant of its Jacobian matrix $A_{ai}={\partial }_i{u}_a$. In Einstein summation convention,Differentiate the supplied expression for $V_i$ by the product rule. The term differentiating ${\partial }_j{u}_a$ vanishes because ${\partial }_i{\partial }_j{u}_a$ is symmetric in $i,j$, while ${ϵ}_{ijk}$ is antisymmetric; the term differentiating ${\partial }_k{u}_b$ vanishes similarly. Only differentiation of $u_c$ remains:The last equality follows by the cyclic relabelling $(i,j,k)\mapsto (j,k,i)$, which has positive sign. This proves the Jacobian determinant as a divergence identity.

For a composition $\mathbf{w}=\mathbf{u}\circ \mathbf{v}$, the chain rule givesTaking determinants yieldsFor the first-order perturbation, use spatial derivative control in perturbations of the identity map: the usual smooth-family interpretation gives $D{\mathbf{u}}_t=I+tD\mathbf{F}+o(t)$ locally. Multilinearity of the determinant shows that its linear term comes from choosing one column of $tD\mathbf{F}$ and all other columns from $I$. This gives the trace:The regularity qualification matters: a pointwise $o(t)$ remainder for maps does not by itself give an $o(t)$ remainder for spatial derivatives. For instance $(x+t^2\sin (x/t^2),y,z)$ has a uniform $o(t)$ displacement and $\mathbf{F}=0$, but its Jacobian determinant at $x=0$ is $2$ for $t\ne 0$. Thus the first expansion needs the remainder in $C^1$ locally, or an equivalent smoothness assumption on the family.

Under the subsequent group law, ${\mathbf{u}}_0$ is the identity and ${\mathbf{u}}_{-t}$ is the inverse of ${\mathbf{u}}_t$. Write $j(t,\mathbf{x})=J[{\mathbf{u}}_t](\mathbf{x})$. Applying the composition formula to ${\mathbf{u}}_{t+h}={\mathbf{u}}_h\circ {\mathbf{u}}_t$ and the short-time expansion gives the Jacobian evolution of a smooth flow:HenceIf only the pointwise generator expansion is assumed initially, the group law still gives ${\partial }_t{\mathbf{u}}_t(\mathbf{x})=\mathbf{F}({\mathbf{u}}_t(\mathbf{x}))$. Uniqueness and smooth dependence for this smooth ordinary differential equation make it a smooth flow map, thereby justifying the differentiated expansion used above. Thus a divergence-free generator hasBy the change of variables formula, these maps preserve ordinary Euclidean volume and orientation. This is the Euclidean instance of a volume-preserving vector field.

Away from the origin, the given vector field is the gradient of the globally single-valued radius:The sphere has centre $(5,3,2)$, at distance $\sqrt{38}$ from the origin, and radius one. Its intersection circle is therefore wholly away from the origin, since every point on the sphere has distance at least $\sqrt{38}-1>0$. Along the closed curve the chain rule gives $\mathbf{F}\cdot d\mathbf{x}=dr$. Consequently, independently of orientation,The intersection is nondegenerate: the centre's distance from the plane is $4/\sqrt{35}<1$. No parametrization of the resulting circle is needed because $r$ is a potential of a conservative vector field throughout the punctured space.

The chain rule along a smooth curve isDifferentiating the circular parametrization gives the tangent vectorThus the chain rule turns the angular derivative of $u$ into $-y{u}_x+x{u}_y$. The given partial differential equation makes this equal to $u$, so along every portion of the circle inside the upper half-plane, $U(t)=u(\mathbf{x}(t))$ satisfies the ordinary differential equation $U^{\prime }=U$. Its solution is $U(t)=U(t_0)e^{t-t_0}$.

For the two specified points take $a=\sqrt{2}$. The initial point corresponds to $t_0=\pi /4$ and the terminal point to $t_1=3\pi /4$. The entire intervening arc lies in $y>0$, so no continuation outside the domain is used. The requested value isThese circular characteristic curves also show locally that the general solution has the polar-coordinate form $u(r,\theta )=C(r)e^{\theta }$, with $0<\theta <\pi$; there is no requirement of angular periodicity on the upper half-plane.

A conservative vector field on ${\mathbb{R}}^n$ is a vector field $\mathbf{V}=\nabla F$ for a globally defined smooth potential of a conservative vector field $F$. The chain rule givesso its line integral is path independent.

Green theorem states that, for a bounded planar region $D$ with piecewise smooth boundary and continuously differentiable $P,Q$ on a neighbourhood of its closure,The boundary is positively oriented: the region lies on the left while it is traversed, so outer components run counterclockwise and hole boundaries clockwise.

For the proposed potential of a conservative vector field, the fundamental theorem of calculus and differentiation under the integral sign giveSince $Q_x=P_y$, the last integral equals $P(x,y)-P(x,0)$, including negative $y$ with the usual signed-integral convention. Hence the proposed potential is global andThe fact that the vector field is defined on all of ${\mathbb{R}}^2$ ensures both straight integration segments stay inside its domain; a hole in the domain could obstruct a global potential.

Set $s=\log x$ and $Y(s)=y(e^s)$, so $0\le s\le 1$. By the chain rule, the equation becomes $Y^{\prime \prime }+2Y^{\prime }+(1+\lambda )Y=0$. The substitution $Y=e^{-s}v$ reduces it further toMultiplication by $\overline{v}$ and integration by parts gives $\lambda {\int }_0^1|v{|}^2={\int }_0^1|v^{\prime }{|}^2>0$ for a nonzero solution, ruling out zero and negative eigenvalues. For $\lambda =k^2>0$, the first boundary condition gives $v=A\sin (ks)$, and the second gives $k=n\pi$, $n=1,2,\dots$. Thus all eigenpairs arewith arbitrary nonzero scalar multiples of each eigenfunction. These are the complete list since all solutions of the constant-coefficient differential equation have been considered; the transformed sine system is also complete in the corresponding weighted Hilbert space.

Parametrize a streamline by $\mathbf{r}(s)$ with $d\mathbf{r}/ds=\mathbf{u}(\mathbf{r}(s))$. The chain rule givesTherefore $H$ is constant along every streamline. A different regular parametrization merely multiplies this derivative by a scalar and gives the same conclusion. For time-dependent flow, the claim concerns streamlines at a fixed time, rather than Lagrangian trajectories.