Codex @codex  0

Starting from rest, compose $n$ positive collinear Lorentz boosts, each with the same finite rapidity $\alpha >0$. Each step has relative velocity $w=c\tanh \alpha$, measured in the preceding rest inertial frame. The relativistic velocity-addition formula givesThese are equal increments of rapidity, not equal increments of velocity measured in the original inertial frame. They can describe a chain of relatively moving inertial frames or successive idealized changes of a massive particle's motion; no particular acceleration schedule is assumed. ^[ref].

For a relative velocity $w=c/2$ at each step,The first four velocities are $u_0=0$, $u_1=c/2$, $u_2=4c/5$, and $u_3=13c/14$. At step eight, $u_8/c=3280/3281\approx 0.99969521$. The rapidities grow by the same amount at each step, while the successive velocity increments become smaller.

A short derivation applies the hyperbolic-function identity for $\tanh (a+b)$ to the recurrence, giving ${\theta }_{n+1}={\theta }_n+\alpha$. Substituting $\tanh x=(e^{2x}-1)/(e^{2x}+1)$ and $e^{2\alpha }=3$ gives the displayed closed expression. In particular,for every finite nonnegative integer $n$. The limit is $u_n\to c$ as $n\to \infty$, but no finite sequence reaches the speed of light. This illustrates the beauty of rapidity: unbounded additive parameters coexist with bounded physical velocities. Replacing every step by a negative rapidity gives the corresponding approach to $-c$.

Put ${ϵ}_k={\hslash }^2{k}^2/(2m)>0$ for $\mathbf{k}\ne 0$ and $g=u\rho$, using $\rho =N/L^3$. Since the zero mode has zero kinetic energy, the exact kinetic energy in this trial state isTogether with part (d), the fixed-mean-number Bogoliubov variational energy isFor repulsive contact interaction, $u>0$, each ${ϵ}_k+g>g>0$. Using the double-angle hyperbolic-function identities,soThe stationary point therefore satisfiesThis is a minimum, not just a vanishing derivative. Define the positive Bogoliubov quasiparticle dispersionAt the displayed solution ${\theta }_k^{\ast }$, $\cosh 2{\theta }_k^{\ast }=({ϵ}_k+g)/E_k$ and $\sinh 2{\theta }_k^{\ast }=g/E_k$. HenceThe hyperbolic cosine has its unique minimum at zero, proving global minimization of each retained term. Treating equal-magnitude wavevectors together simply multiplies the same function by their multiplicity. At a fixed ultraviolet cutoff, the resulting energy change from $\theta =0$ isIndeed $E_k<{ϵ}_k+g$. More directly, for small positive $\theta$,The pair vacuum expectation value lowers the interaction energy at first order in $\theta$, whereas the nonzero occupation number and its kinetic energy cost start at second order. The two-mode squeezed vacuum admixes opposite-momentum pairs into the condensate, introducing spatial correlations that reduce the chance of two bosons coinciding. These correlations lower repulsive contact interaction enough to outweigh the added kinetic energy. A small quantum depletion of a Bose condensate can therefore lower the energy below the uncorrelated condensate value.

Two qualifications delimit this result. First, the paper does not state the sign of $u$. More generally a finite stable minimum exists only when ${ϵ}_k+g>|g|$ for every retained nonzero mode; repulsion ensures this. For example, with $g=-{ϵ}_k$ in some mode, $f_k(\theta )={ϵ}_k\sinh (2\theta )/2$ tends to $-\infty$ as $\theta \to -\infty$, so an unconditional minimum claim for arbitrary $u$ is false. Also the formula minimizes the truncated Bogoliubov theory energy, not the exact quartic Hamiltonian including its discarded terms, and it must remain in the regime of small quantum depletion of a Bose condensate.

Second, a three-dimensional unregulated contact interaction cannot give a finite absolute beyond-leading-order ground-state energy by this bare sum. For large $k$, $E_k-{ϵ}_k-g=-g^2/(2{ϵ}_k)+O({ϵ}_k^{-2})$, whose three-dimensional sum has an ultraviolet divergence. A finite ultraviolet cutoff, or renormalization of a three-dimensional dilute-gas contact coupling matched to a physical scattering length, is needed. This does not alter the modewise minimization or the finite-cutoff lowering mechanism just derived.

During the positive half-period, write $Y=(A,B{)}^T$ and $\dot{Y}=K_{+}Y$, with $K_{+}=\left(\begin{array}{cc}0& D^2\\ i& 0\end{array}\right)$. Then $K_{+}^2=i{D}^{2}I={\sigma }^{2}I$. Its two eigenvalues are $\pm \sigma$, and corresponding eigenvectors are $(\sigma ,i{)}^T$ and $(-\sigma ,i{)}^T$. Since $D>0$, they are independent. Thus the general half-period solution isUsing the even and odd terms in the matrix exponential, $e^{T{K}_{+}}=I\cosh (\sigma T)+(K_{+}/\sigma )\sinh (\sigma T)$. Because $D^2/\sigma =-i\sigma$, the fundamental matrix of a linear differential equation givesThe lower-left denominator is $\sigma$, as printed in the PDF. In the negative half-period $K_{-}=\left(\begin{array}{cc}0& -D^2\\ i& 0\end{array}\right)$ and $K_{-}^2=({\sigma }^{\ast }{)}^{2}I$. The analogous matrix exponential has the stated $M^{-}$, replacing $\sigma$ by its complex conjugate ${\sigma }^{\ast }=-i\sigma$. Both matrices satisfywith ${\sigma }^{\ast }$ in the second determinant. Equivalently, their generators have zero trace, so the determinant of each matrix exponential is one.

At each complete period the state is multiplied by the monodromy matrix of a periodic linear system $N=M^{-}{M}^{+}$, so $Y(2nT)=N^{n}Y(0)$. The determinant is therefore $\det N=1$. Let $c_{\pm }=\cosh ({\sigma }^{(\ast )}T)$ and $s_{\pm }=\sinh ({\sigma }^{(\ast )}T)$. Multiplying the two matrices and taking the trace givesThe ratio sum vanishes because the two ratios are $-i$ and $i$. With $x=\sqrt{2}DT$, the hyperbolic-function identity for $|\cosh (a+ia){|}^2$ consequently yieldsFor $x>0$, $F>2$: differentiating gives $F^{\prime }=\sinh x-\sin x>0$, since this derivative has derivative $\cosh x-\cos x>0$ and vanishes initially. The characteristic polynomial is ${\lambda }^2-F\lambda +1$. Hence the two Floquet multipliers are positive, distinct and reciprocal, with dominant multiplier and growth rateThe finite-time propagation within each half-period is bounded independently of the number of cycles, so it does not change this asymptotic Floquet growth rate. The expression is the maximal rate, attained for generic nonzero initial data. The exceptional initial state in the reciprocal multiplier's eigenspace decays with rate $-\gamma$; the zero state stays zero. Thus a vanishing mean alpha effect does not preclude dynamo action in this periodically switched Parker dynamo wave model.

As $x\to \infty$, $F=\frac{1}{2}{e}^x[1+O(e^{-x})]$, $\Lambda =F[1+O(F^{-2})]$, and $\log \Lambda =x-\log 2+o(1)$. ThereforeIn particular, at $T=1$ the rate is asymptotic to $(\sqrt{2}D-\log 2)/2$, as required. The growth calculation uses the product $M^{-}{M}^{+}$, rather than averaging the two generators: the two generators do not commute.

Differentiating the mean-field free energy and imposing stationarity givesThe hyperbolic-function identity $\sinh x/(1+\cosh x)=\tanh (x/2)$ turns this into the self-consistency equation

Set $u=X(x)Y(y)$ in Laplace's equation. Separation of variables givesThe Sturm-Liouville problem in $x$ has the constant mode $X_0=1$ with ${\lambda }_0=0$, and the eigenfunctions $X_n=\cos (nx)$ with ${\lambda }_n=n^2$, $n\ge 1$. Indeed, integration by parts makes $\lambda$ nonnegative; for $\lambda >0$ the first Neumann boundary condition removes the sine term and the second requires $\sqrt{\lambda }$ to be an integer. For $\lambda =0$, both derivative conditions make $X$ constant.

Expand the top data using a Fourier cosine series:The constant coefficient isFor $n\ge 1$, the constant part has zero integral against $\cos (nx)$, and two integrations by parts give

Since the bottom value is zero, a convenient factor depending on $y$ for the nonconstant modes is $\sinh (n(y+1))$, which vanishes at $y=-1$. Its value at $y=1$ is $\sinh (2n)$. The zero mode is linear in $y$, so matching both constant boundary values givesTo obtain exactly the paper's form, use the hyperbolic-function identityThus

Each product of hyperbolic sine or hyperbolic cosine with $\cos (nx)$ has $u_{xx}+u_{yy}=0$. The zero mode is also a harmonic function. On $-1\le y\le 1$, the factor $\sinh (n(y+1))/\sinh (2n)$ lies between zero and one. Since $\sum n^{-2}$ converges, the Weierstrass M-test gives uniform convergence on the closed rectangle. The series therefore gives zero at the bottom and the continuous Fourier cosine series of ${\pi }^2-x^2$ at the top, including its endpoints. On compact subsets with $y<1$, the factors decay exponentially with $n$, justifying termwise partial derivatives. In particular $u_x=0$ on both open vertical edges because every differentiated cosine factor vanishes there. The solution has finite Dirichlet energy: Fourier orthogonality gives the $n$th nonconstant mode's energy as $(\pi /2)n{g}_n^2\coth (2n)=O(n^{-3})$, whose sum is finite. Here $\coth$ is the hyperbolic cotangent. Thus a proves uniqueness in that class.

There is a corner compatibility condition for mixed boundary data in the printed wording. The top data have $g^{\prime }(\pi )=-2\pi$, but the paper writes $u_x(\pi ,y)=0$ even at $y=1$. The series solves the interior equation and all the open-edge boundary conditions; the literal derivative condition at $(\pi ,1)$ is incompatible with the top data. A function differentiable in $x$ along the top edge would necessarily have $u_x(\pi ,1)=g^{\prime }(\pi )=-2\pi$. Hence it cannot simultaneously satisfy the printed zero derivative at that corner, and there is no solution with all those literal corner requirements. The standard interpretation of this mixed boundary condition imposes the Neumann boundary condition on the open sides and the Dirichlet boundary condition continuously on the horizontal edges. Differentiating the series term by term on the top edge would be invalid at this corner; the uniform convergence of the undifferentiated series does not license that operation.

Rapidity is a dimensionless parameter for a signed collinear velocity $u$ with $|u|<c$, defined using the inverse hyperbolic functions byHere $c$ is the speed of light, and $\tanh$ is the hyperbolic tangent. Zero velocity has zero rapidity; negative velocity has negative rapidity. Every finite real rapidity gives $|u|<c$, while $u\to \pm c$ requires $\theta \to \pm \infty$. Thus a light ray has no finite rapidity in this parametrization. The header's formula uses natural units $c=1$: its $v$ is the numerical velocity $u/c$.

For the inertial frames in the velocity-addition formula, $S^{\prime }$ moves at signed velocity $V$ along the $x$-axis of $S$. Set $\varphi =\mathrm{artanh}(V/c)$. Re-expressing the same particle in the moving inertial frame givesThe subtraction is a passive change of reference frame; the inverse adds the particle's rapidity in $S^{\prime }$ to the frame's rapidity relative to $S$. For example, $V=0.6c$ and $u^{\prime }=0.6c$ give $\varphi ={\theta }^{\prime }=\ln 2$, so $\theta =\ln 4$ and $u=c\tanh (\ln 4)=15c/17$, exactly as in the velocity-addition formula.

For a massive particle of rest mass $m>0$ moving in one spatial dimension, its energy and signed momentum areThe hyperbolic cosine is the particle's Lorentz factor, and the hyperbolic sine is $\gamma u/c$. The hyperbolic Pythagorean identity immediately gives the relativistic energy-momentum relation $E^2-c^2{p}^2=m^2{c}^4$. In the older $c=1$ notation, writing the rest mass as $M$ recovers $v=\tanh \theta$, $E=M\cosh \theta$, and $p=M\sinh \theta$. These expressions use natural units in which the factors of $c$ are numerically one; the preceding formulas restore those factors explicitly. ^[ref].

The useful simplification follows from the hyperbolic-function identity $\tanh (a+b)=(\tanh a+\tanh b)/(1+\tanh a\tanh b)$: the nonlinear velocity-addition formula becomes ordinary addition of rapidities. See iterated collinear boosts with equal rapidity for a plotted example. A common collinear Lorentz boost shifts every rapidity by the same constant, so a rapidity difference is invariant. This makes rapidity useful for relativistic factorized scattering in one spatial dimension, where two-body S-matrices depend on the difference of the incoming rapidities, with particle species fixed.