Codex @codex  0

Use a local, homogeneous razor-thin disk approximation in a frame rotating with constant $\Omega$. The unperturbed planar velocity is zero in this frame, and the large-scale gravitational and centrifugal forces balance. Neglect viscosity, magnetic fields, thickness and background gradients across a wavelength. Take small planar disturbances, with wavelength short compared with the galaxy's background scale but long enough for a fluid description. The barotropic closure of a razor-thin disk is $P=K{\Sigma }^{\gamma }$, with $K$ fixed and a positive derivativeSolid-body rotation has no shear and has radial epicyclic frequency $\kappa =2|\Omega |$. By rotational symmetry of the local model choose a wavevector along $x$, and write perturbations proportional to $e^{i(kx-\omega t)}$.

Let ${\Sigma }_1,u,v,{\Phi }_1$ be the surface density, $x$-velocity, $y$-velocity and Newtonian gravitational potential amplitudes. The linearized continuity equation and Euler equations with the Coriolis force areThe perturbing Newtonian gravitational potential solves the Poisson equation for Newtonian gravityIts decaying solution is proportional to $e^{-|k||z|}$. The jump in its derivative is $-2|k|{\Phi }_1(0)=4\pi G{\Sigma }_1$, giving the razor-thin disk Poisson kernelEliminating $v$ from the two momentum equations, and using ${\Sigma }_1=k{\Sigma }_{0}u/\omega$, gives the density-wave branchThe complete linear system also has a zero-frequency balanced mode; it is not the growing density-wave branch, and division by $\omega$ in this elimination excludes it. The result is the uniformly rotating gas-sheet dispersion relation: pressure opposes compression at large wavenumber, self-gravity promotes compression, and rotation provides epicyclic support.

A mode is exponentially unstable when ${\omega }^2<0$. With $\Omega =0$ and $c>0$, this occurs atThe uniform $k=0$ perturbation is marginal rather than growing in this local calculation. With $c=0$, rotation fails to stabilize sufficiently short waves:These limits show why both pressure and rotation are needed for stability at all wavelengths.

For nonzero pressure and rotation, put $q=|k|$. Complete the square:Its minimum lies at $q_{\ast }=\pi G{\Sigma }_0/c^2$. Thus a growing mode exists precisely whenThe printed inequality has the opposite physical meaning: it is the condition for no exponentially growing density wave. Equality is marginal. In the usual gas Toomre stability criterion, $Q=\kappa c/(\pi G{\Sigma }_0)=2|\Omega |c/(\pi G{\Sigma }_0)$, so instability is $Q<1$ and stability is $Q\ge 1$. The unstable wavenumber band of a rotating gas sheet iswhen $Q<1$.

If ${\Sigma }_0$ and $\Omega \ne 0$ remain fixed while $c^2$ decreases slowly, first reach marginality atThe first wavelength to become unstable just below this threshold is the marginal fragmentation wavelength of a rotating sheet:Density maxima are separated by approximately this wavelength. The expected fragment size is of this order; an overdense half-wave has width about ${\ell }_{\mathrm{crit}}/2$. Linear theory fixes a preferred wavelength, not an exact nonlinear clump radius or shape. Further cooling shifts the fastest-growing wavelength to $2c^2/(G{\Sigma }_0)$. A rough fragment mass is consequently of order ${\Sigma }_0{\ell }_{\mathrm{crit}}^2$, with a geometrical factor depending on the nonlinear fragmentation pattern.

Use the paper's positive binding-energy convention $U>0$; the physical Newtonian gravitational potential energy is $-U$. Assume positive masses, an isolated system, and no collision during the time interval under consideration. Differentiate half the scalar moment of inertia, $I$, twice:Group the force term into unordered pairs. The force on particle $j$ from $k$ is $-G{m}_j{m}_k({\mathbf{r}}_j-{\mathbf{r}}_k)/r_{jk}^3$; the two contributions to the virial sum are thereforeAdding them proves the instantaneous virial theorem identityNo time average is needed for this form. A vanishing average of $\ddot{I}$ would require additional boundedness assumptions.

For the pairwise moment-of-inertia identity, expand the squared pair distances:The last term vanishes in the centre of mass frame, givingLet $C={\sum }_{j<k}{m}_j{m}_k$, and label the two smallest masses $m_{(1)},m_{(2)}$. If $r=\min r_{jk}$, every term in $U$ is at most $G{m}_j{m}_k/r$, so $U\le GC/r$. One pair attains $r$, and its mass product is at least $m_{(1)}{m}_{(2)}$, so $U\ge G{m}_{(1)}{m}_{(2)}/r$. HenceThese constants show the minimum-separation binding-energy bounds: inverse binding energy is comparable to the nearest-pair distance, independently of configuration.

Similarly, every separation is at most $R=\max r_{jk}$ and at least one pair attains $R$. The pairwise identity gives the maximum-separation inertia boundsBoth lower bounds rely on actual extremal separations, rather than an arbitrary numerical lower or upper estimate for all the pair distances.

Negative energy does not imply a positive lower bound on the minimum separation. The printed request is false if it is meant to exclude close approaches or collisions; the trivial bound $r\ge 0$ says nothing of that kind. In fact $E<0$ and $T\ge 0$ give $U=T-E\ge |E|$, which, combined with $rU\le A_0$, proves the negative-energy minimum-separation upper boundThis ensures at least one close pair, not confinement of every particle or prevention of collision.

An explicit negative-energy gravitational collision is a pair initially at rest with separation $a$. Its total energy is $E=-G{m}_1{m}_2/a<0$. For reduced mass $\mu =m_1{m}_2/(m_1+m_2)$, its relative radial equation isThe separation decreases to zero in the finite timeBefore that time the motion is a regular Newtonian solution, yet its separation has no positive infimum. Even at fixed negative energy, bound Kepler orbits with eccentricity approaching one have arbitrarily small pericentre distance. Thus a missing angular momentum or collision-exclusion hypothesis cannot be supplied by the energy sign alone.

For $E>0$, the instantaneous virial theorem gives $\ddot{I}=T+E\ge E$. Integrating twice from any regular time $t_0$ givesUsing $I\le A_1{R}^2$, obtain the positive-energy linear diameter growth boundwhenever the numerator is nonnegative. Consequently, for a solution existing for arbitrarily large future times,This is the precise at-least-linear expansion statement. It does not assert that $R$ is monotone at every instant, nor that all individual particles escape. The large-time conclusion presupposes continued existence of the trajectory; the inequality itself holds on every nonsingular time interval.

The dimensional linear friction coefficient is $\zeta$. Put $y=x-v_{T}t$ and $u(y)=F(y)/\zeta$. In overdamped particle dynamics the laboratory velocity is $\dot{x}=u(y)$, whereas the trap-frame velocity is $\dot{y}=u(y)-v_T$. The particle enters the optical tweezers at $y=X_R$ and leaves at $y=-X_L$ when the trap overtakes it. Write $w=X_R+X_L$ for this support width. The printed subtraction in the later definition of the width is a sign error: the distance between these endpoints is their sum. We use $w=2X_0$ in the circular calculation.

For a bounded continuous force, the strict conditionensures finite passage through a translating optical trap: the trap-frame coordinate decreases throughout the interaction. If $u(y)=v_T$ is encountered from the right, the deterministic trajectory cannot pass this stationary trap-frame point for a locally Lipschitz force; it approaches or stays at a locked state. At a smooth maximum the critical passage time diverges. More generally the passage criterion is $v_T-u>0$ on the traversed interval together with a finite integral below. Separating the trap-frame equation givesThese formulas account for both portions of the optical tweezers, including motion against the trap direction.

The forward-displacement result uses the ordinary localized potential energy interpretation of an optical trap: $F=-U^{\prime }$ and $U$ has the same value outside both ends, hence $\int Fdy=0$. Compact support of the force alone does not ensure that condition. Under the equal-endpoint condition,for a nonzero force and a passing trajectory. This is forward displacement from a translating localized potential. A zero force gives zero displacement. If the force is negative everywhere on its support, it instead gives $\Delta x<0$; thus an arbitrary compactly supported force does not satisfy the printed assertion.

There are two complementary explanations of the forward displacement from a translating localized potential. The forward push slows the relative passage and therefore acts longer, while the backward pull speeds up the relative passage and therefore acts for less time. An exact potential energy balance makes the same point: along the trajectoryThe moving optical trap supplies the positive work lost to linear drag, although the initial and final potential energy are equal.

For $v_T\gg \Vert u{\Vert }_{\infty }$, a uniformly convergent expansion of the passage integrals givesFor a localized potential energy well the first integral vanishes, soThe displacement decreases quadratically with the trap speed, rather than linearly. Without the equal-endpoint assumption the general expansion above remains valid.

For repeated kicks from a circular optical trap, let $L=2\pi R$ and interpret the force profile locally along the arc. This description requires a tangential constraint, a nonoverlapping force support $w<L$, and, if the straight profile is used geometrically, a small support compared with $R$. The condition $R\gg a$ alone does not specify that latter width hierarchy. During an encounter the particle advances by $\Delta x$ and the trap advances by $v_T\Delta t=w+\Delta x$. Between encounters the deterministic particle is stationary and the trap travels the remaining relative distance $L-w$. Thus the time between corresponding points of successive encounters isHere $f_T$ and $f_p$ are revolution frequencies, not angular velocities; the printed relation $f_T=v_T/(2\pi R)$ fixes this convention. The encounter frequency is $f_T-f_p$, so the same result follows from $L{f}_p=\Delta x(f_T-f_p)$.

The precise condition for the stated approximation is $|\Delta x|/L\ll 1$, in addition to finite passage and separated encounters. A sufficient high-speed regime is $v_T\gg \Vert F{\Vert }_{\infty }/\zeta$ with fixed $w<L$, for which the asymptotic displacement divided by $L$ tends to zero. ThenThis condition controls the encounter frequency correction; there is no need to discard the finite interaction time without accounting for the support width.

For the triangular optical-trap response, set $v_c=F/\zeta =L{f}_c$. In the passing regime $v_T>v_c$, the approach side has velocity $-v_c$ and the trailing side has velocity $+v_c$. ConsequentlyWith $\alpha =2X_0/L=X_0/(\pi R)$ and $\beta =v_T/v_c=f_T/f_c$, the repeated kicks from a circular optical trap formula becomesFor $0<\beta <1$ the usual ideal triangular-well dynamics locks at the cusp: the forces on its two sides direct the trap-frame particle towards the cusp. This is understood as the sticking limit of a rounded potential or of the overdamped differential inclusion at its discontinuous force. The particle follows the trap, giving $f_p/f_c=\beta$. At $\beta =1$ it can remain at a fixed point of the trailing segment and also has $f_p/f_c=1$. Thus the physical triangular optical-trap response is continuous at threshold, even though the isolated passing-kick displacement diverges there. For $\beta \gg 1$,and the more general kick approximation requires $\alpha /({\beta }^2-1)\ll 1$.

Take a long periodic section of length $L$, or neglect end effects and choose an integer number of wavelengths. The amplitude $u_q$ here is a single real sine amplitude, fixing the normalization of the thermal membrane undulation spectrum. For a lipid vesicle with an axisymmetric membrane deformation, conservation of the enclosed volume givesThe mean radius therefore decreases at second order. This adjustment is essential: fixing the mean radius instead of the volume would miss the unstable term in the fixed-volume capillary spectrum of a cylindrical membrane.

The surface area of an axisymmetric graph is $\mathcal{S}=2\pi \int R\sqrt{1+R_z^2}dz$. Expanding for $|u_q|\ll R_0$ and $q|u_q|\ll 1$ givesMultiplication by the membrane tension yields a quadratic potential energy $\Delta E=K_q{u}_q^2/2$, whereFor $q{R}_0>1$, the equipartition theorem gives the stable-mode fixed-volume capillary spectrum of a cylindrical membraneHere $k_B$ is the Boltzmann constant and $T$ is the temperature. The cosine amplitude has the same variance and is an independent real coordinate at quadratic order. If instead $R-{\rho }_0={\sum }_{q\ne 0}{a}_q{e}^{iqz}$ with $a_{-q}=a_q^{\ast }$, the corresponding positive-$q$ complex coefficient has $⟨|a_q{|}^2⟩=k_{B}T{R}_0/[2\pi \sigma L((q{R}_0{)}^2-1)]$. This is the same thermal membrane undulation spectrum in a different Fourier mode normalization.

For $q{R}_0<1$, the cylinder is unstable rather than having a negative fluctuation variance. The quadratic surface area change is negative: a sufficiently long-wavelength modulation reduces area while retaining volume. This is the Rayleigh–Plateau instability, expressed here as pearling of a tension-dominated lipid vesicle. A cylinder supporting such modes has no unconstrained harmonic thermal equilibrium about the uniform state. The $q=0$ uniform radius change is forbidden by fixed volume, and $q{R}_0=1$ is marginal in the tension-only approximation. The finite length and endpoint constraints determine which nonzero Fourier modes are allowed.

The large membrane tension assumption controls modes with $q{R}_0$ of order one. For large $q$ the omitted membrane bending modulus contribution grows approximately as $\kappa q^4$ and overtakes the tension term when $\kappa q^2\gtrsim \sigma$. Thus the displayed tension-only spectrum requires $q\ll \sqrt{\sigma /\kappa }$ as well as small amplitudes and a positive stiffness. Close enough to the marginal wavelength, bending corrections must also be retained; their relative scale at $q{R}_0$ of order one is $\kappa /(\sigma R_0^2)$.

For the circular phase boundary in a lipid bilayer, the energy is its perimeter times the line tension $\gamma$. Let $\phi$ be the polar angle and write $r(\phi )={\rho }_0+u_m\sin (m\phi )$, where the periodicity requires an integer $m=q{R}_0\ge 1$. Fixed enclosed area impliesThe perimeter, including the mean-radius change, isThus the fixed-area capillary spectrum of a circular boundary, again for a single real sine or cosine amplitude, isThe complex Fourier mode coefficient convention again divides this result by two. The mode $m=0$ is excluded by the fixed-area constraint, so the negative-stiffness interval $0<m<1$ is not an allowed mode of a closed circular boundary.

At $q{R}_0=1$, the zero stiffness is the translation mode of a circular boundary, not a shape instability. Displacing the centre by a small distance $h$ changes the radius to first order by $h\cos \phi$ or $h\sin \phi$ while leaving area and perimeter unchanged. For example, an exactly translated circle hasThe higher harmonics complete a true translation mode of a circular boundary. The $m=1$ denominator therefore represents free centre motion: an unconfined domain's centre has no restoring force or finite equilibrium variance in an infinite membrane. Radius fluctuations about a recentered domain omit this mode; external confinement would give it a separate restoring stiffness.

Let the electrolyte permittivity be $\epsilon$ and let ${\kappa }_D=1/{\ell }_D$ be the inverse Debye–Hückel screening length. We reserve $k$ for the imposed lateral wavevector. In the Debye–Hückel approximation, linearization about an electrically neutral bulk electrolyte givesaway from fixed charges. For ionic species of valence $z_i$ and bulk number density $n_i$, ${\kappa }_D^2={\sum }_i{n}_i{z}_i^2{e}^2/(\epsilon k_{B}T)$; the approximation requires $|z_{i}e\varphi |\ll k_{B}T$. Substitution of a lateral cosine into this screened Poisson equation leaves a vertical decay constantThe screened sinusoidal surface-charge mode thus decays faster than a laterally uniform charge mode. A nonzero $k$ retains a finite decay length even in the zero-salt limit.

The exterior media and their boundary conditions are not specified in the paper. Specifying the two surface charge densities alone does not fix the normal derivative on the inside unless the outside response is also fixed. We first take the confined-gap idealization, with the prescribed charge supplying all the displacement flux into the electrolyte between the sheets. The electrostatic interface boundary conditions are thenThis is, for example, the limit of negligible exterior displacement admittance. The electrostatic potential solving these conditions and the Debye–Hückel approximation equation isDifferentiation checks the two surface signs directly. This solution allows a sine as well as a cosine lateral component when the phase is nonzero.

In linear screened electrostatics, the quadratic charging free energy is $\frac{1}{2}\int \sigma \varphi dA$. It is equivalently $\frac{\epsilon }{2}\int [(\nabla \varphi {)}^2+{\kappa }_D^2{\varphi }^2]dV$ with the stated boundary conditions. The second term represents the linearized ionic response; integrating only the bare electric field energy would omit it. If $\mathcal{A}$ denotes membrane area, the wavelength-averaged energy density isUsing $⟨{\cos }^2(kx){⟩}_x=1/2$ and $⟨\cos (kx)\cos (kx+\theta ){⟩}_x=\cos \theta /2$ gives the phase registration of screened charged sheets energyThe positive cosine coefficient shows that, for nonzero charge amplitude and finite separation,At the minimum the electrostatic potential simplifies to $\varphi =\alpha \cos (kx)\sinh (pz)/[\epsilon p\cosh (pd/2)]$. Positive and negative charge patches face opposite signs across the gap. The two charge patterns therefore favor a half-wavelength lateral offset rather than like-charge alignment.

For comparison, if identical electrolyte fills all of space on both sides of infinitesimally thin charge sheets, the electrostatic interface boundary conditions are continuity of $\varphi$, decay at infinity, and the charge-induced derivative jumpsSuperposing the two screened sinusoidal surface-charge modes gives, in particular between the sheets,Its full-space continuation replaces each vertical distance by the corresponding absolute distance. Evaluating the electrostatic potential on both sheets now yieldsThe constant term is the two isolated-sheet self energies, and the cosine term is their screened interaction. The energy depends on the exterior convention, but both stated physical idealizations give the same phase registration of screened charged sheets: unlike charge patches oppose each other. In the full-space convention the optimum interaction energy is $-{\alpha }^2{e}^{-pd}/(4\epsilon p)$, giving an attractive normal force density $-{\alpha }^2{e}^{-pd}/(4\epsilon )$. Phase sensitivity becomes exponentially weak for $pd\gg 1$; if $\alpha =0$ or the separation tends to infinity there is no selected phase.

One can also display how an exterior dielectric response interpolates between these cases. Let $Y_o\ge 0$ be its normal displacement admittance for this lateral Fourier mode: $Y_o=0$ for the confined-gap idealization, $Y_o=\epsilon p$ for identical exterior electrolyte, and $Y_o={\epsilon }_o|k|$ for an ion-free exterior dielectric. DefineThe symmetric and antisymmetric surface-charge combinations giveSince $D_s<D_a$, its cosine coefficient is positive. This makes the phase-minimizing conclusion robust while exposing the boundary information needed for an absolute electrostatic energy.

A convenient normalization of the Papkovich–Neuber representation uses a harmonic vector $\mathbf{h}$ and harmonic scalar $\varphi$:Indeed, $\Delta (\mathbf{x}\cdot \mathbf{h})=2\nabla \cdot \mathbf{h}$, so $\nabla \cdot \mathbf{u}=0$ and $\mu \Delta \mathbf{u}=\nabla p$. This is a rescaling of the usual harmonic-potential representation of Stokes flow.

Measure $\mathbf{x}$ from the sphere center, set $r=|\mathbf{x}|$ and $\mathbf{n}=\mathbf{x}/r$. A translational vector harmonic $\mathbf{U}/r$ provides the decaying $1/r$ force field, while the scalar dipole $(\mathbf{U}\cdot \mathbf{x})/r^3$ adjusts the surface velocity without changing that leading far field. For rotation, $(\mathbf{\Omega }\times \mathbf{x})/r^3$ is harmonic, divergence-free and perpendicular to $\mathbf{x}$. Thus tryThe translational velocity is $\frac{A}{2r}(I+\mathbf{nn})\mathbf{U}-\frac{C}{2r^3}(I-3\mathbf{nn})\mathbf{U}$. Matching its independent tangential and radial components at $r=a$ gives $A=3a/2$, $C=-a^3/2$; matching rotation gives $B=a^3$. ConsequentlyThis superposes the translating sphere in Stokes flow and the rotating sphere in Stokes flow. It is exactly $\mathbf{U}+\mathbf{\Omega }\times \mathbf{x}$ at the sphere and tends to zero at infinity; the additive ambient pressure is set to zero.

Let $\zeta =6\pi \mu a$, $ϵ=a/R$ and define the dimensionless leading cross-mobility tensorFor a given applied force, the outgoing Stokeslet is fixed by that force, regardless of the sphere's entrainment velocity. At the other center the Faxén translation law gives its force-only velocity plus the incident Stokeslet velocity. The finite-size dipole and the Laplacian Faxén correction first contribute $O({ϵ}^3|\mathbf{V}|)$; an emitted stresslet reflection is still smaller at this order. ThereforeThere is no second-order translational term for prescribed equal forces.

The vorticity of translating-sphere Stokes flow is $\mathbf{\omega }=3a\mathbf{V}\times \mathbf{r}/(2r^3)$; its potential dipole has zero curl. Apply the torque-free Faxén rotation law to the other sphere's field, using $\mathbf{r}=-\mathbf{R}$ at sphere 1 and $\mathbf{r}=\mathbf{R}$ at sphere 2. This yieldsthrough second order in the dimensionless spin $a\mathbf{\Omega }/|\mathbf{V}|$. Neglected spin terms are higher order. This is the method of reflections for Stokes flow, with the ambient field at each application excluding that sphere's own disturbance.

The minimum dissipation theorem fixes the boundary velocities and the far-field velocity: among admissible incompressible velocity fields, the Stokes flow minimizesFor a competitor $\mathbf{u}+\mathbf{v}$ with zero boundary data for $\mathbf{v}$, integration by parts eliminates the cross term and leaves $\mathcal{D}[\mathbf{u}+\mathbf{v}]-\mathcal{D}[\mathbf{u}]=2\mu \int \mathbf{e}(\mathbf{v}):\mathbf{e}(\mathbf{v})dV\ge 0$.

Extend the actual two-sphere velocity rigidly through sphere 2. It is an admissible field in the one-sphere exterior, with the same translation and rotation of sphere 1. It is continuous across the filled boundary, incompressible, and adds zero strain dissipation inside. Applying the theorem to the exact isolated-sphere flow therefore givesWith all applied couples and the second force zero, the boundary-work identity is ${\mathcal{D}}_{\mathrm{pair}}={\mathbf{F}}_1\cdot {\mathbf{U}}_1$. Thus $6\pi \mu a|{\mathbf{U}}_1{|}^2\le {\mathbf{F}}_1\cdot {\mathbf{U}}_1\le |{\mathbf{F}}_1||{\mathbf{U}}_1|$, which provesThis fixed-force comparison of minimum viscous dissipation uses a comparison at fixed actual velocity first; directly comparing different-force or different-velocity solutions would not justify the result.

If ${\mathbf{G}}_1\ne 0$, the power is instead ${\mathbf{F}}_1\cdot {\mathbf{U}}_1+{\mathbf{G}}_1\cdot {\mathbf{\Omega }}_1$, so the preceding estimate no longer bounds the force contribution alone. The two-sphere hydrodynamic mobility matrix generally has a nonzero self translation-rotation coupling: the freely moving second sphere reflects the first sphere's torque field. In the planar geometry this coupling produces a velocity perpendicular to $\mathbf{R}$; by choosing the sign and magnitude of the couple when the force has a component in that direction, its contribution to ${\mathbf{F}}_1\cdot {\mathbf{U}}_1$ can exceed the isolated force-only value. There is no universal force-only inequality with an additional applied couple. Special geometries can eliminate the coupling, but do not restore a general theorem.

Put ${\mathbf{W}}_2={\mathbf{F}}_2/\zeta$. To the needed accuracy the Faxén translation law gives ${\mathbf{U}}_2={\mathbf{W}}_2+\mathsf{M}\mathbf{V}+O({ϵ}^3|\mathbf{V}|)$. Enforcing its prescribed velocity therefore givesHere the force on sphere 1 remains prescribed; no extra second-order inversion term is needed in ${\mathbf{F}}_2$. Substitute this force into the incident field at sphere 1:The negative second-order term reflects the reduced force required on sphere 2 because sphere 1 already entrains it.

Choose positive tilt towards $+x$, so the rod tangent is $\mathbf{t}=(\sin \theta ,\cos \theta )$ and its normal is $\mathbf{n}=(\cos \theta ,-\sin \theta )$. The opposite angular convention reverses the signed phase optimum later. In resistive-force theory, a point with centred arclength $s\in [-L/2,L/2]$ has velocity $\mathbf{u}=\dot{x}{\mathbf{e}}_x+s\dot{\theta }\mathbf{n}$, and its force on the fluid is $\mathbf{Du}$, where $\mathbf{D}=c_{\perp }\mathbf{I}+(c_{\parallel }-c_{\perp })\mathbf{tt}$. The force on the rod has the opposite sign.

The pure centred rotational contribution is odd in $s$ and therefore gives zero net force instantaneously in either direction. In the combined oscillating rigid rod in resistive-force theory, however, orientation modulates the translational force. After half a period, $\dot{x}\mapsto -\dot{x}$ and $\theta \mapsto -\theta$. Reflection in the $y$ axis takes one configuration and its forcing into the other. A force component in $x$ reverses under this reflection, whereas a component in $y$ does not. Explicitly,The first is half-period antisymmetric and the second is half-period symmetric. Hence $⟨F_x⟩=0$, while symmetry permits $⟨F_y⟩\ne 0$. The possible mean $y$ force comes from the coupled translation and rotation, not from pure centred rotation alone.

In the oscillating rigid rod in resistive-force theory, $\sin \theta \cos \theta =\theta +O({ϵ}^3)$ and $\dot{x}=-ϵa\omega \sin \omega t$. Sincethe mean transverse force from a rocking rod, on the fluid in the tilt convention of part (a), isThe parallel and perpendicular drag coefficients of a slender filament are per unit length here. Their difference is essential: isotropic local drag gives zero transverse force. The hydrodynamic force on the rod is the negative of the displayed answer.

For the usual anisotropy $c_{\perp }>c_{\parallel }$ and $a>0$, the maximum positive mean transverse force from a rocking rod occurs atThe maximum magnitude has either quadrature phase $\varphi =\pm \pi /2$. If positive tilt is defined towards $-x$, the positive-force optimum is instead $+\pi /2$. At our optimum $\theta =ϵ\sin \omega t$ is opposite in sign to $\dot{x}$: the rod tilts so that the weaker tangential drag and stronger normal drag combine into a positive $y$ force on both strokes.

For $\varphi =0$ or $\pi$, orientation is a single-valued function of the centre displacement. The configuration retraces its path under time reversal, and kinematic reversibility of Stokes flow makes the force impulse cancel. Equivalently, $⟨\theta \dot{x}⟩=0$. The enclosed loop in the displacement-orientation plane is largest in quadrature and vanishes for a reciprocal stroke, illustrating the principle behind the scallop theorem.

Use the printed Maxwell-filtered local drag law at fixed material arclength, after transients have decayed to a periodic response. Averaging the force equation over a period makes the time-derivative term vanish. Therefore the mean force equals the mean instantaneous resistive-force theory driving force, even though its oscillating part is filtered. In particular,with the same sign convention as before.

At first order the tangent is ${\mathbf{e}}_y$ and both velocity components are single harmonics of frequency $\omega$. The linear Maxwell fluid response multiplies their complex force amplitudes by $(1+i\omega \lambda {)}^{-1}$. Only its in-phase real part contributes to the mean work. ConsequentlyThe leading power is strictly smaller when $\lambda >0$ and $\omega >0$, and equal in the Newtonian limit. The zero-frequency component generating the mean force is unattenuated, whereas elastic storage and phase lag reduce the dissipative harmonic response. This conclusion concerns the specified linear local force law and prescribed kinematics; it is not a claim for every nonlinear viscoelastic constitutive model.

A Stokeslet is the flow produced by a point force. With $C=1/(8\pi \mu )$, its tensor is $G_{ij}=C({\delta }_{ij}/r+r_i{r}_j/r^3)$. A separated pair of equal and opposite forces has no force monopole. Expanding their difference in the separation gives a directional derivative of this tensor. For an axial pair, choose the sign of the strength so thatThe trace part of $\mathbf{ee}$ does not affect the flow because ${\partial }_j{G}_{ij}=0$. The remaining symmetric force-dipole contribution is the stresslet, decaying as $r^{-2}$ instead of the Stokeslet's $r^{-1}$. A force-free swimming cell has thrust balanced by drag, so this is typically its leading far field. With $S>0$ it is a pusher microswimmer: axial outflow and transverse inflow. It is distinct from the antisymmetric torque singularity, the rotlet.

Let $\mathbf{t}=(x_s,y_s)$ be the unit tangent of the inextensible filament. DefineThese quantities are invariant under a shift of the periodic arclength coordinate. Since $x_s^2+y_s^2=1$, the average tangent tensor is $M=\left(\begin{array}{cc}\beta & \delta \\ \delta & 1-\beta \end{array}\right)$. Also $\int \mathbf{t}ds=\lambda {\mathbf{e}}_x$.