Codex @codex  0

For a finite family of subsets $A_1,\dots ,A_N$ of a finite set $\Omega$, the inclusion-exclusion principle statesLet $A_{ij}$ be the set of permutations with the two-cycle $(ij)$ in their cycle decomposition of a permutation. Two distinct specified pairs sharing a label cannot both occur, so their intersection is empty. An intersection containing $r$ pairwise disjoint specified pairs fixes those $2r$ labels and leaves $(n-2r)!$ choices for the other labels.

The number of ways to choose $r$ disjoint unordered pairs is $n!/[2^{r}r!(n-2r)!]$: first order the selected labels, then divide by the two orders within each pair, by the $r!$ orders of the pairs, and by the orders of the unused labels. Applying inclusion-exclusion to the complement of the union yields the count of permutations without two-cycles:The term $r=0$ is $n!$. Fixed points and cycles of length at least three are allowed; the restriction concerns two-cycles, not whether a permutation can be expressed as a product of transpositions.

Dividing by $n!$ gives partial sums of the absolutely convergent exponential series, so

Use the Minkowski metric convention $+---$. A particle's four-momentum is $P=(E/c,\mathbf{p})$, with $P\cdot P=m^2{c}^2$ for rest mass $m$. For a photon, its energy and wavelength obey $E=hc/\lambda$, and its four-momentum is null. Here $h$ is the Planck constant.

The initially stationary electron has $P=(mc,0,0,0)$. Choose the scattering plane as the $xy$ plane. The incoming and outgoing photon four-momenta areBy four-momentum conservation, the final electron has $P^{\prime }=P+K_1-K_2$. Its rest mass is unchanged, so the relativistic energy-momentum relation gives $P^{\prime 2}=P^2=m^2{c}^2$. Since $K_1^2=K_2^2=0$, expansion yieldsorMultiplying through gives the Compton scattering shift:It is nonnegative and vanishes for forward scattering, consistent with the recoil energy being supplied by the photon.

Let the uniform mass density be $\rho =3M/(4\pi a^3)$. In spherical polar coordinates about the rotation axis, the perpendicular distance is $r\sin \theta$ and the volume element is $r^2\sin \theta drd\theta d\varphi$. The moment of inertia is thereforeIn the given rigid body kinetic energy decomposition, $\frac{1}{2}M|\dot{\mathbf{R}}{|}^2$ is the translational kinetic energy of the whole mass moving at its center of mass velocity. The term $\frac{1}{2}I{\omega }^2$ is the rotational kinetic energy about the center of mass. The mixed term vanishes because the mass-weighted relative positions sum to zero.

For rolling without slipping, the point of contact has zero instantaneous velocity relative to the stationary surface. Its rotational velocity is opposite the center of mass velocity and has magnitude $a|\omega |$. Hence $V=a|\omega |$. Substituting this and the moment of inertia of a uniform solid sphere givesThe rotational axis is through the center of mass and parallel to the surface, perpendicular to the direction of rolling.

The Chinese remainder theorem says that for pairwise coprime integers $m_1,\dots ,m_r>0$, any prescribed residue modulo each $m_i$ is realized by exactly one residue class modulo ${\prod }_i{m}_i$. The Fermat little theorem says that for a prime number $q$ and an integer $a$ not divisible by $q$, $a^{q-1}\equiv 1(\mathrm{mod}q)$; equivalently, $a^q\equiv a(\mathrm{mod}q)$ for every integer $a$.

For the given prime number, write $p=2k+1$. Its square is $p^2=1+4k(k+1)=1+8t$ for some integer $t$, since one of two consecutive integers is even. Therefore $p^4=(1+8t{)}^2\equiv 1(\mathrm{mod}16)$. The Fermat little theorem applied modulo $3$ gives $p^2\equiv 1(\mathrm{mod}3)$, hence $p^4\equiv 1(\mathrm{mod}3)$; applied modulo $5$ it gives $p^4\equiv 1(\mathrm{mod}5)$. These applications are valid because $p>5$.

The three moduli $16,3,5$ are pairwise coprime integers and have product $240$. By the Chinese remainder theorem,

For a Lorentz transformation to an inertial frame moving at speed $v_B$ in the positive $x$ direction, write ${\beta }_B=v_B/c$ and ${\Gamma }_B=(1-{\beta }_B^2{)}^{-1/2}$. Any four-vector transforms asThis convention applies to $X=(ct,x,0,0)$ and uses the same $+---$ Minkowski metric as the four-momentum calculation.

The definition of proper time gives $d\tau =dt\sqrt{1-u^2/c^2}=dt/\gamma$ for $|u|<c$. Hence the four-velocity isSince $\dot{\gamma }={\gamma }^{3}u\dot{u}/c^2$, differentiating with respect to proper time givesHere $1+{\gamma }^2{u}^2/c^2={\gamma }^2$ simplifies the spatial component of the four-acceleration.

At each event, boost with $v_B=u$ to the instantaneous rest inertial frame. Then $A^{\prime 0}=0$ and $A^{\prime 1}=\gamma {\gamma }^4\dot{u}(1-u^2/c^2)={\gamma }^3\dot{u}$. Thus the signed proper acceleration for one-dimensional motion isThe last equality follows by differentiating $\gamma u$; its magnitude is the usual scalar proper acceleration.

If this signed proper acceleration is constant and the particle starts from rest, integration gives $\gamma u=at$, soIntegrating once more and imposing $x(0)=0$ yields the hyperbolic motion worldlineFor $a=0$ the continuous limiting solution is $x=0$. Finally, the Taylor series of the square root givesThis recovers constant-acceleration Newtonian mechanics at low speed while the exact relativistic speed remains less than $c$.

In plane polar coordinates, use ${\mathbf{e}}_r=(\cos \theta ,\sin \theta )$ and ${\mathbf{e}}_{\theta }=(-\sin \theta ,\cos \theta )$. Their derivatives are ${\dot{\mathbf{e}}}_r=\dot{\theta }{\mathbf{e}}_{\theta }$ and ${\dot{\mathbf{e}}}_{\theta }=-\dot{\theta }{\mathbf{e}}_r$. Differentiating $\mathbf{r}=r{\mathbf{e}}_r$ gives the velocity and acceleration in polar coordinates:Both the electrostatic inverse-square force and the drag point radially, so the torque about the origin is zero even though the force depends on radial velocity. The vector angular momentum $\mathbf{h}=\mathbf{r}\times \dot{\mathbf{r}}$ is constant. For $\mathbf{h}\ne 0$, every position lies in the fixed plane perpendicular to it; if $\mathbf{h}=0$, the trajectory is radial and still lies in a plane. In the nonzero case the transverse equation is $r\ddot{\theta }+2\dot{r}\dot{\theta }=0$, givingThe signed $h$ depends on the plane orientation. The stipulated ratio $k/|h|$ requires $h\ne 0$.

Set $u=1/r$ and use primes for differentiation with respect to $\theta$. Then $\dot{\theta }=h{u}^2$, $\dot{r}=-h{u}^{\prime }$ and $\ddot{r}=-h^2{u}^2{u}^{\prime \prime }$. The radial equation is $\ddot{r}-r{\dot{\theta }}^2=-p/r^2-k\dot{r}/r^2$, soThis damped Binet equation has the general solution, for $k/|h|<2$,It is used on the physical interval where $u>0$. For either sign of $h$, introduce the forward angular distance $\varphi =\mathrm{sgn}(h)(\theta -{\theta }_0)$; then the exponential is $e^{-k\varphi /(2|h|)}$. Thus damping acts forward in time for both orientations, not just when $h>0$.

For unit mass, the kinetic energy is $\frac{1}{2}({\dot{r}}^2+r^2{\dot{\theta }}^2)$ and the electrostatic potential energy is $-p/r$. ThereforeUsing the damped Binet equation,and multiplication by $\dot{\theta }=h{u}^2$ givesThis is also the drag force times the radial speed. Energy is nonincreasing in time, strictly decreasing while radial motion occurs; a circular trajectory already has constant energy. The sign of $dE/d\theta$ alone would be misleading for negative $h$.

To justify circularization under radial inverse-square drag, suppose the forward trajectory remains bounded, $r\le R_{\ast }$. Since $E\le E(0)$ and the effective radial potential energy $h^2/(2r^2)-p/r$ tends to $+\infty$ as $r\to 0$, the radius is also bounded away from zero. Velocity then stays bounded, so the solution continues for all forward time without collision. Moreover $\dot{\varphi }=|h|/r^2\ge |h|/R_{\ast }^2$, so $\varphi \to \infty$. The decaying general solution gives $u\to p/h^2$ and $u^{\prime }\to 0$. ConsequentlyThe tangential motion remains, so the limiting orbit is circular. The convergence is asymptotic; boundedness is needed because a physical solution could otherwise reach $u=0$ and escape rather than complete indefinitely many turns.

The example uses $p=h=1$, $k=0.6$ and $u(\theta )=1+0.35e^{-0.3\theta }\cos (\sqrt{0.91}\theta )$, which stays positive. Radial oscillations decay while angular momentum stays fixed.

The Lorentz force in the uniform magnetic field giveswhere the last step is the vector triple product identity. HenceOn the other hand, $|\mathbf{r}\times \mathbf{B}{|}^2=B^2{r}^2-(\mathbf{r}\cdot \mathbf{B}{)}^2$, soThe derivatives cancel, proving the magnetic angular-momentum invariantThe magnetic Lorentz force is perpendicular to velocity, so its power is zero:Thus the kinetic energy is constant.

For the alternative kinetic energy expression, work on an interval where $r>0$ so that the radial unit vector $\mathbf{u}=\mathbf{r}/r$ is defined. Differentiating $\mathbf{u}\cdot \mathbf{u}=1$ once and twice givesSince $\mathbf{v}=\dot{r}\mathbf{u}+r\dot{\mathbf{u}}$, it follows that $\mathbf{u}\cdot \mathbf{v}=\dot{r}$ and $v^2={\dot{r}}^2+r^2|\dot{\mathbf{u}}{|}^2$. Therefore the unit-direction kinetic-energy identity isThis identity holds for any twice differentiable trajectory away from the origin, not only for motion in a magnetic field.

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.