Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Cm Math Central

1 · Prerequisites

2 · Summary

The central-force layer of the non-relativistic classical mechanics collection, all of it mathematics: the reduction of the non-radial central-force equation to polar data, the Kepler problem up to the determination of its conic, period and one-dimensional action variable, and the finite-range and Coulomb scattering layer with its cross-section densities. Every item on this page carries domain: mathematics and depends only on mathematical suppliers; nothing here asserts that the equations describe a physical interaction, that incident flux is realized by counted events, or that any empirical conclusion follows. The physical reading of masses, forces, energies and scattering measurements remains with the consuming classical-mechanics pages.

The reduction layer begins with Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential, which defines the central-force equation μr′′=−(U′(∣r∣)/∣r∣)r for a curve r:I→R3 with values in J0⊆(0,∞), together with the angular momentum vector j=μ r×r′, the energy and the effective potential UJeff. Local polar coordinates for a curve with nonvanishing position in the plane: a differentiable angle function with prescribed derivative supplies local polar coordinates ρ,σ of a planar curve as a twice differentiable rotation, and Conservation of the angular momentum vector along a central-force trajectory and its orthogonality to position and velocity proves that j is constant and orthogonal to r and r′. Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation then derives the polar identities ρ2θ′=J/μ, the radial equation μρ′′=−(UJeff)′(ρ) and the conserved energy E=12μρ′2+UJeff(ρ). Circular orbits of a non-radial central-force trajectory: the critical-point condition and the squared frequency of the linearised radial equation characterises circular orbits by (UJeff)′(r0)=0, exhibits them explicitly, and reduces the linearised radial equation to the constant-coefficient oscillator, whose solutions are classified by Solutions of v′′+ω2v=0 are exactly the functions Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)). The Kepler specialization proceeds through Binet's equation for the attractive inverse-square field: d2u/dφ2+u=μκ/J2 for u=1/ρ, which derives u′′+u=μκ/J2 for u=1/ρ along the angular variable, and The Kepler orbit equation is a conic: u=c(1+ecos⁡(φ−θ0)) with E=(μκ2/2J2)(e2−1), which proves u=c(1+ecos⁡(φ−θ0)) and the exact relation E=μκ22J2(e2−1) connecting the sign of the energy with the eccentricity. The algebraic prerequisites The vector triple product identity u×(v×w)=⟨u,w⟩v−⟨u,v⟩w and Oriented orthonormal triads in R3: existence over any unit vector, and the orthonormal expansion with no fourth orthogonal direction are proved locally: the first is the vector triple product identity, the second constructs an oriented orthonormal triad over any unit vector and proves that no four orthonormal vectors exist in R3, together with the orthonormal expansion.

The Kepler layer proper is The subcritical Kepler conic is an ellipse with semiaxes a=p/(1−e2), b=a1−e2, focus at the origin and area πab; its area is also 12∫02πϱ2, which proves that the subcritical conic is the ellipse with centre (−ae,0) and semiaxes a=p/(1−e2), b=a1−e2, with area λ2(Eell)=πab=12∫02πϱ(φ)2 dφ; the item declares the Axiom of Countable Choice used by its affine change-of-variables and polar-decomposition arguments. Completeness and periodicity of the negative-energy Kepler orbit, with T=2πμab/J and T2=4π2μa3/κ proves completeness of the negative-energy orbit, the periodicity of the motion and the exact period value T=2πμab/J=2πμ/κ a3/2, hence T2=4π2μa3/κ; it also declares the Axiom of Countable Choice. The one-dimensional action layer is The one-dimensional action variable of a family of two-turning-point energy levels, which defines Ja(E)=1π∫a(E)b(E)2m(E−U(q)) dq for a two-simple-turning-point family, and The action variable of a two-simple-turning-point family is C1 with derivative T(E)/(2π), the period quadrature over 2π, which proves that the action variable is C1 with Ja′(E)=T(E)/(2π) by the inverse function theorem at the turning points, a C3 endpoint rescaling and differentiation under the integral sign.

The measure part of the scattering layer assumes countable choice. The scattering layer is built on Scattering data: the incident flux measure, the solid-angle element, deflection branches and the cross-section density, which defines the finite-range data, the incident flux measure b db dφ, the solid-angle measure sin⁡θ dθ dφ, deflection branches and their cross-section density b(θ)/ ⁣(sin⁡θ ∣θ′(b(θ))∣) as abstract measures and functions; no counting of events and no empirical claim belongs to the definition. Finite-range scattering: the simple outer turning point, the two global branches and the free tails assumes a simple outer turning point and proves the two unique global branches with free tails, Existence of the asymptotic velocities and the impact-parameter identity J=μb∣v−∣ constructs the asymptotic velocities v± and the impact-parameter identity J=μb∣v−∣, and The scattering angle advance Δ and the deflection relation: the total polar advance is 2Δ and the velocity directions make the cosine cos⁡(π−2Δ) proves the advance quadrature Δ=∫r∗∞J dr/(r22μ(E−V(r))) and the deflection relation modulo whole revolutions. One-dimensional substitution on a deflection branch: the outgoing angular measure has density b∣db/dθ∣, hence cross-section density b∣db/dθ∣/sin⁡θ performs the one-dimensional substitution on compact interior intervals and a Borel measure-uniqueness argument, giving the full open-branch pushforward identity for nonnegative or absolutely integrable Borel test functions and the outgoing angular measure density b∣b′∣ and the cross-section density b∣b′∣/sin⁡θ, additively over finitely many disjoint branches. Repulsive Coulomb scattering: the unique outer turning point, global non-colliding motion, the asymptotic impact distance and Δ=π/2−arctan⁡(κ/(2Eb)) separately treats the long-range the repulsive Coulomb potential, where the effective potential is strictly decreasing with a single simple turning point and the motion is global and non-colliding, with genuine asymptotic velocity and impact-distance limits and angle advance Δ=π/2−arctan⁡(κ/(2Eb)); finally The Rutherford deflection map, its inverse, the cross-section density (κ/4E)2csc⁡4(θ/2) and the divergence of the total cross-section proves, for fixed energy and J(b)=2μE b, the deflection map θdef(b)=2arctan⁡(κ/(2Eb)), its inverse b(θ)=κ2Ecot⁡(θ/2), the density (κ/4E)2csc⁡4(θ/2) and the divergence of the total cross-section at forward angles together with finiteness over every compact angular interval. The companion page cm-math-central-examples records exact finite checks of the action derivative, of the Kepler ellipse data and of the Rutherford density, and a zero-angular-momentum collision counterexample that marks the excluded radial degenerate case.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential

Definition

Let μ>0 be a real number, let J0⊆(0,∞) be an open interval, let U:J0→R be of class C2 (Ck Euclidean maps and diffeomorphisms), let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), and let r:I→R3 be twice differentiable with ∣r(t)∣∈J0 for every t∈I (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral). The curve r is a central-force trajectory for (μ,U) when

μ r′′(t)=−U′(∣r(t)∣)∣r(t)∣ r(t)(t∈I),

where ∣x∣:=∥x∥2 is the Euclidean norm (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn) and U′ is the derivative of U on the open interval J0.

Attached objects. The angular momentum vector of such a trajectory is the curve j:I→R3,

j(t):=μ r(t)×r′(t),

with × the cross product of R3 (The cross product in R3); its energy is the function E:I→R,

E(t):=12μ ∥r′(t)∥22+U(∣r(t)∣);

and for every J≥0 the effective potential of momentum magnitude J is the function UJeff:J0→R,

UJeff(ρ):=U(ρ)+J22μ ρ2.

Why the objects are defined. Since ρ>0 throughout J0, the right-hand side of the central-force equation is defined, the map ρ↦J2/(2μρ2) is C2 on (0,∞) (a quotient of C2 functions with nowhere vanishing denominator), and UJeff is again C2; in particular second derivatives of UJeff are available. The vector j(t) is defined for every t∈I because r is differentiable there.

Mathematical status. The data μ, U and J are real numbers and functions; the equation and the attached objects are mathematical statements about curves. Nothing here asserts that the equation describes a physical interaction, and the names are conventional labels for the mathematical objects only.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Local polar coordinates for a curve with nonvanishing position in the plane: a differentiable angle function with prescribed derivative

Statement

Let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let t0∈I, and let ρ,σ:I→R be twice differentiable with (ρ(t),σ(t))≠(0,0) for every t∈I. Here twice differentiable means that f and its derivative f′:I→R are differentiable at every point of I in the domain-relative sense of The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set; f′′=(f′)′, including one-sided limits at included endpoints. Put

r(t):=ρ(t)2+σ(t)2>0(t∈I).

Then there is a twice differentiable θ:I→R with

ρ=rcos⁡θ,σ=rsin⁡θon I,

and every twice differentiable ϑ:I→R with ρ=rcos⁡ϑ and σ=rsin⁡ϑ on I satisfies

ϑ′(t)=ρ(t)σ′(t)−σ(t)ρ′(t)r(t)2(t∈I).

Facts & Assumptions

Given: An interval I⊆R with at least two points, a point t0∈I, and twice differentiable functions ρ,σ:I→R with (ρ,σ)≠(0,0) on I.

[L3]

t↦(cos⁡t,sin⁡t) is a bijection from [0,2π) onto the unit circle S1={(x,y):x2+y2=1}, and sin⁡2+cos⁡2=1 (t↦(cos⁡t,sin⁡t) is a bijection from [0,2π) onto the real unit circle, Parity and the Pythagorean identity for sine and cosine).

Proof

technique · direct
1.1givenL1algebra

By [L1], r is differentiable with r′=(ρρ′+σσ′)/r. Since ρ,σ are twice differentiable and r>0, this quotient is differentiable by the product and quotient rules; hence r is twice differentiable.

1.2givenL3choose

The vector (ρ(t0)/r(t0),σ(t0)/r(t0)) satisfies (ρ/r)2+(σ/r)2=1 by the definition of r, so by [L3] there is a real θ0 with cos⁡θ0=ρ(t0)/r(t0) and sin⁡θ0=σ(t0)/r(t0).

2.1givenstep 1.1step 1.2L2L4construct

Define g:I→R by g:=(ρσ′−σρ′)/r2. The numerator is continuous and r2>0 is continuous, so g is continuous on I by [L4]; hence by [L4] the function θ(t):=θ0+∫t0tg(s) ds is differentiable on I with θ′=g. Moreover g is differentiable by the product and quotient rules, since ρ,σ are twice differentiable and r2>0; therefore θ is twice differentiable with θ′=g.

2.2givenstep 1.2L3algebra

Define A:=ρcos⁡θ+σsin⁡θ and B:=σcos⁡θ−ρsin⁡θ on I. Evaluating at t0 with the choice of θ0 and [L3] gives A(t0)=ρ(t0)2/r(t0)+σ(t0)2/r(t0)=r(t0) and B(t0)=σ(t0)ρ(t0)/r(t0)−ρ(t0)σ(t0)/r(t0)=0.

3.1step 1.1step 2.1L1L2algebra

Put c:=r′/r on I. Differentiating with [L2] and substituting θ′=g gives the two identities ρ′+σθ′=cρ and σ′−ρθ′=cσ: indeed ρ′+σθ′=ρ′+σ(ρσ′−σρ′)/r2=(ρ′ρ2+ρσσ′)/r2=ρ(ρρ′+σσ′)/r2=cρ, and σ′−ρθ′=σ′−ρ(ρσ′−σρ′)/r2=(σ′σ2+ρσρ′)/r2=σ(ρρ′+σσ′)/r2=cσ, using [L1]. Hence A′=(ρ′+σθ′)cos⁡θ+(σ′−ρθ′)sin⁡θ=cA, B′=(σ′−ρθ′)cos⁡θ−(ρ′+σθ′)sin⁡θ=cB.

4.1step 2.2step 3.1L1L2L5algebra

By the quotient rule of [L2] and r′=cr from [L1], (A/r)′=A′r−Ar′r2=cAr−Acrr2=0,(B/r)′=cBr−Bcrr2=0on I. The functions A/r and B/r are continuous on I; by [L5] and step 2.2 they are the constants 1 and 0: A=r and B=0 on I.

5.1step 4.1L3algebra

Solving the two equations A=ρcos⁡θ+σsin⁡θ=r and B=σcos⁡θ−ρsin⁡θ=0 for (ρ,σ) gives ρ=Acos⁡θ−Bsin⁡θ=rcos⁡θ and σ=Asin⁡θ+Bcos⁡θ=rsin⁡θ on I: the displayed linear system has the invertible coefficient matrix (cos⁡θsin⁡θ−sin⁡θcos⁡θ) of determinant cos⁡2θ+sin⁡2θ=1.

6.1givenL1L2algebra∎

Let ϑ:I→R be twice differentiable with ρ=rcos⁡ϑ, σ=rsin⁡ϑ. Differentiating the two identities and eliminating r′ gives ρσ′−σρ′=rcos⁡ϑ (r′sin⁡ϑ+rϑ′cos⁡ϑ)−rsin⁡ϑ (r′cos⁡ϑ−rϑ′sin⁡ϑ)=r2ϑ′, so ϑ′=(ρσ′−σρ′)/r2 on I.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Solutions of v′′+ω2v=0 are exactly the functions Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0))

Statement

Let ω>0, let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let t0∈I, and let v:I→R be twice differentiable with

v′′(t)+ω2v(t)=0(t∈I).

Here twice differentiable means that v and v′:I→R are differentiable at every point of I in the domain-relative sense of The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set; v′′=(v′)′, with one-sided limits at included endpoints.

Then there are reals A,B with

v(t)=Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0))(t∈I).

Conversely, for all reals A,B the function t↦Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)) is twice differentiable on R and satisfies v′′+ω2v=0 there.

Facts & Assumptions

Given: A real ω>0, an interval I⊆R with at least two points, a point t0∈I, and a twice differentiable v:I→R with v′′+ω2v=0 on I.

[L1]

Sine and cosine are differentiable with (sin⁡u)′=cos⁡u and (cos⁡u)′=−sin⁡u; the chain rule differentiates the compositions t↦cos⁡(ω(t−t0)), t↦sin⁡(ω(t−t0)) with derivatives −ωsin⁡(ω(t−t0)) and ωcos⁡(ω(t−t0)) (The derivatives of sine and cosine are cosine and minus sine, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c)).

[L4]

sin⁡2+cos⁡2=1 identically (Parity and the Pythagorean identity for sine and cosine).

Proof

technique · direct
1.1givenL1L2construct

Write φ(t):=ω(t−t0) and define w:=vcos⁡φ−v′ωsin⁡φ,z:=vsin⁡φ+v′ωcos⁡φ on I. Every factor is differentiable (including v′, by the twice differentiable hypothesis), so w and z are differentiable by [L2].

1.2givenL1L2algebra

Differentiating and using v′′=−ω2v: w′=v′cos⁡φ−ωvsin⁡φ−v′′ωsin⁡φ−v′cos⁡φ=−ω2v+v′′ωsin⁡φ=0, and likewise z′=v′sin⁡φ+ωvcos⁡φ+v′′ωcos⁡φ−v′sin⁡φ=ω2v+v′′ωcos⁡φ=0 on I.

2.1step 1.1step 1.2L3

By [L3] and step 1.2 there are reals A,B with w=A and z=B on I.

3.1step 1.1step 2.1L4algebra

Multiplying out and using [L4], wcos⁡φ+zsin⁡φ=v(cos⁡2φ+sin⁡2φ)=v on I; hence v=Acos⁡φ+Bsin⁡φ, which is the asserted representation.

4.1L1L2algebra∎

Conversely, if u(t):=Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)), then [L1] and [L2] give u′=−ωAsin⁡φ+ωBcos⁡φ and u′′=−ω2Acos⁡φ−ω2Bsin⁡φ=−ω2u on I, so u′′+ω2u=0.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

The vector triple product identity u×(v×w)=⟨u,w⟩v−⟨u,v⟩w

Statement

For all u,v,w∈R3, u×(v×w)=⟨u,w⟩ v−⟨u,v⟩ w.

Facts & Assumptions

Given: Vectors u,v,w∈R3 with coordinates u=(u1,u2,u3), v=(v1,v2,v3), w=(w1,w2,w3).

[L1]

The cross product is defined coordinatewise by u×v=(u2v3−u3v2, u3v1−u1v3, u1v2−u2v1) (The cross product in R3).

[L2]

The Euclidean inner product is ⟨u,w⟩=∑k=13ukwk (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L3]

Finite sums may be expanded, reordered and combined (Laws of finite sums and finite products).

Proof

technique · direct
1.1L1L2L3algebra

First component. By [L1], (u×(v×w))1=u2(v×w)3−u3(v×w)2=u2(v1w2−v2w1)−u3(v3w1−v1w3), which expands to v1(u2w2+u3w3)−w1(u2v2+u3v3); since u2w2+u3w3=⟨u,w⟩−u1w1 and u2v2+u3v3=⟨u,v⟩−u1v1 by [L2], this equals v1(⟨u,w⟩−u1w1)−w1(⟨u,v⟩−u1v1)=(⟨u,w⟩v−⟨u,v⟩w)1, which is the 1-component of the asserted identity.

2.1L1L2L3algebra

Second and third components. The same computation, with the coordinate labels cyclically relabelled and the components of the cross product read off in the corresponding cyclic order, gives (u×(v×w))2=(⟨u,w⟩v−⟨u,v⟩w)2 and (u×(v×w))3=(⟨u,w⟩v−⟨u,v⟩w)3: each is the same finite expansion of [L1] as in step 1.1 with (1,2,3) cyclically replaced, and [L2], [L3] are unaffected by the relabelling.

3.1step 1.1step 2.1L1algebra∎

Two vectors of R3 with equal components are equal, so steps 1.1 and 2.1 prove the identity in all three coordinates.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Oriented orthonormal triads in R3: existence over any unit vector, and the orthonormal expansion with no fourth orthogonal direction

Statement

Call a triple (e1,e2,e3) of vectors of R3 an oriented orthonormal triad when ⟨ei,ej⟩=δij for all i,j∈{1,2,3} and e1×e2=e3. Then:

  1. Existence. For every e3∈R3 with ∣e3∣=1 there are e1,e2∈R3 such that (e1,e2,e3) is an oriented orthonormal triad.
  2. Dimension bound. Every orthonormal list f1,…,fm∈R3 (pairwise orthogonal unit vectors) has m≤3; consequently, for every oriented orthonormal triad (e1,e2,e3) and every x∈R3, x=⟨x,e1⟩e1+⟨x,e2⟩e2+⟨x,e3⟩e3,∣x∣2=⟨x,e1⟩2+⟨x,e2⟩2+⟨x,e3⟩2.

Throughout, ⟨⋅,⋅⟩ and ∣⋅∣ are the Euclidean inner product and norm (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), × is the cross product (The cross product in R3), and s1,s2,s3 is the standard basis of R3 (The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0).

Facts & Assumptions

Given: A unit vector e3∈R3, standard basis vectors s1,s2,s3, and arbitrary orthonormal lists as in the statement.

[L1]

The triple product identity u×(v×w)=⟨u,w⟩v−⟨u,v⟩w holds for all u,v,w∈R3 (The vector triple product identity u×(v×w)=⟨u,w⟩v−⟨u,v⟩w).

[L2]

The cross product is bilinear and alternating and is orthogonal to both factors; and ∣u×v∣2=∣u∣2∣v∣2−⟨u,v⟩2 (The cross product is bilinear, alternating, and orthogonal to both factors, The squared cross-product norm is the Gram determinant of two vectors, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L3]

Relabel the standard basis indexed 0,1,2 in the supplier as s1,s2,s3. This basis is orthonormal: ⟨si,sj⟩=δij, and every x∈R3 is x=∑j=13⟨x,sj⟩sj with ∣x∣2=∑j=13⟨x,sj⟩2, by the coordinate definition of the Euclidean inner product and the standard basis (The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L4]

Bessel's inequality: for an orthonormal family (f1,…,fm) in an inner product space (Real and complex inner product spaces, with the inner product linear in the first argument, Orthogonal vectors and subspaces, orthogonal and orthonormal sets, and orthonormal bases) and every vector x, ∑i=1m⟨x,fi⟩2≤∣x∣2 (The Bessel inequality for an arbitrary orthonormal family).

[L5]

Finite sums may be reordered and combined (Laws of finite sums and finite products); and two distinct standard basis vectors are not parallel, since they have distinct, nonvanishing coordinates.

Proof

technique · direct
1.1givenL3L5

At most one standard basis vector is parallel to e3: if si=λe3 and sj=μe3 with i≠j and λ,μ≠0, then si=(λ/μ)sj, which is false because si and sj have 1 in different coordinates. Hence some k with sk not parallel to e3 exists.

1.2givenL3L4L5

Let f1,…,fm∈R3 be orthonormal, and suppose m≥4. Applying Bessel's inequality [L4] to each standard basis vector sj gives ∑i=1m⟨sj,fi⟩2≤∣sj∣2=1. Summing over j=1,2,3 and exchanging the finite sums [L5] gives ∑i=1m∑j=13⟨sj,fi⟩2≤3.

2.1step 1.1L3algebra

Put v:=sk−⟨sk,e3⟩e3. Then v≠0: otherwise sk=⟨sk,e3⟩e3 is parallel to e3, contrary to step 1.1. Set e1:=v/∣v∣. Then ⟨e1,e3⟩=0 and ∣e1∣=1, because ⟨v,e3⟩=⟨sk,e3⟩−⟨sk,e3⟩∣e3∣2=0 and ∣v∣≠0.

2.2step 1.2L3L5

For each i, the standard-basis expansion of [L3] applied to fi gives ∑j=13⟨sj,fi⟩2=∣fi∣2=1. Hence the left-hand side of step 1.2 equals ∑i=1m1=m≥4, contradicting the bound ≤3 of step 1.2. Therefore no orthonormal list in R3 has four or more members.

3.1step 2.1L1L2algebra

Set e2:=e3×e1. By [L2], ⟨e2,e3⟩=0=⟨e2,e1⟩ and ∣e2∣2=∣e3∣2∣e1∣2−⟨e3,e1⟩2=1; and e1×e2=e1×(e3×e1)=⟨e1,e1⟩e3−⟨e1,e3⟩e1=e3 by the triple product identity [L1]. So (e1,e2,e3) is an oriented orthonormal triad, which is clause 1.

3.2givenstep 2.2L3algebra

Let (e1,e2,e3) be an oriented orthonormal triad and x∈R3. Put w:=x−∑i=13⟨x,ei⟩ei. Then ⟨w,ej⟩=⟨x,ej⟩−∑i⟨x,ei⟩δij=0 for each j, by bilinearity of the inner product. If w≠0, then e1,e2,e3,w/∣w∣ would be four pairwise orthogonal unit vectors, contradicting step 2.2; hence w=0, which is the expansion.

4.1step 3.2L3algebra∎

With the expansion in hand, ∣x∣2=⟨∑i⟨x,ei⟩ei,∑k⟨x,ek⟩ek⟩=∑i,k⟨x,ei⟩⟨x,ek⟩δik=∑i⟨x,ei⟩2, by bilinearity and orthonormality. This completes clause 2.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Conservation of the angular momentum vector along a central-force trajectory and its orthogonality to position and velocity

Statement

Let μ>0, let J0⊆(0,∞) be an open interval, let U:J0→R be of class C2, let I⊆R be an interval with at least two points, and let r:I→R3 be a central-force trajectory for (μ,U) with angular momentum vector j=μ r×r′ (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential). Then:

  1. j is differentiable on I with j′(t)=0 for every t∈I; hence j is constant on I;
  2. for every t∈I, ⟨r(t),j(t)⟩=0 and ⟨r′(t),j(t)⟩=0.

Facts & Assumptions

Given: A central-force trajectory r for (μ,U) on an interval I with at least two points.

[L1]

For differentiable curves u,v:I→R3 the cross product u×v is differentiable with (u×v)′=u′×v+u×v′; moreover sums, and scalar multiples of curves, differentiate by the usual rules (Differentiation of sums, scalar multiples, inner products, cross products and squared norms of differentiable curves in R3).

[L2]

The trajectory satisfies r′′=−(U′(∣r∣)/(μ∣r∣))r on I, where ∣x∣=∥x∥2 (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L3]

The cross product is bilinear and alternating, so u×u=0 and (u×v)⊥u, (u×v)⊥v; orthogonality means vanishing Euclidean inner product (The cross product is bilinear, alternating, and orthogonal to both factors, The cross product in R3, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

Proof

technique · direct
1.1givenL1algebra

By [L1] the curve j=μ(r×r′) is differentiable with j′=μ(r′×r′+r×r′′)=μ r×r′′.

2.1step 1.1L2L3algebra

Substituting [L2] and using alternation [L3], j′=μ r×(−U′(∣r∣)μ∣r∣ r)=−U′(∣r∣)∣r∣ r×r=0on I.

3.1step 2.1L4

Let a,b∈I with a<b. Then [L4] applied to j gives ∫abj′=j(b)−j(a); since j′ is the zero curve by step 2.1, the integral vanishes and j(b)=j(a). As a,b were arbitrary, j is constant on I.

4.1L3algebra∎

For every t∈I, [L3] gives ⟨r(t),j(t)⟩=μ⟨r(t),r(t)×r′(t)⟩=0 and ⟨r′(t),j(t)⟩=μ⟨r′(t),r(t)×r′(t)⟩=0: the cross product is orthogonal to each factor.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation

Statement

Let μ>0, let J0⊆(0,∞) be an open interval, let U:J0→R be of class C2, let I⊆R be an interval with at least two points, and let r:I→R3 be a central-force trajectory for (μ,U) whose angular momentum j=μ r×r′ is a nonzero constant vector (Conservation of the angular momentum vector along a central-force trajectory and its orthogonality to position and velocity); write J:=∣j∣>0. Let (e1,e2,e3) be an oriented orthonormal triad with e3=j/J (Oriented orthonormal triads in R3: existence over any unit vector, and the orthonormal expansion with no fourth orthogonal direction).

Then there are twice differentiable ρ:I→(0,∞) and θ:I→R with

r(t)=ρ(t)(cos⁡θ(t) e1+sin⁡θ(t) e2)(t∈I),

and for every such pair (ρ,θ):

  1. ρ(t)=∣r(t)∣ for all t;
  2. ρ2θ′=J/μ on I, so θ′>0;
  3. the effective-potential equation μρ′′=−(UJeff)′(ρ) holds on I, with UJeff=U+J2/(2μρ2) (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential);
  4. E:=12μ∥r′∥22+U(∣r∣) equals 12μρ′2+UJeff(ρ) and is constant on I.

Facts & Assumptions

Given: A central-force trajectory r for (μ,U) with constant nonzero angular momentum j=Je3, J>0, an oriented orthonormal triad (e1,e2,e3), and twice differentiable ρ,θ with r=ρ(cos⁡θ e1+sin⁡θ e2).

[L1]

For the oriented orthonormal triad: every x∈R3 is x=⟨x,e1⟩e1+⟨x,e2⟩e2+⟨x,e3⟩e3 and ∣x∣2=∑i⟨x,ei⟩2 (Oriented orthonormal triads in R3: existence over any unit vector, and the orthonormal expansion with no fourth orthogonal direction).

[L2]

The angular momentum is orthogonal to r and to r′: ⟨r,e3⟩J=⟨r,j⟩=0 and ⟨r′,j⟩=0 (Conservation of the angular momentum vector along a central-force trajectory and its orthogonality to position and velocity).

[L3]

Polar lift: for twice differentiable ρ1,ρ2:I→R with (ρ1,ρ2)≠(0,0) there is a twice differentiable angle function such that ρ1=r0cos⁡θ, ρ2=r0sin⁡θ with r0=ρ12+ρ22 (Local polar coordinates for a curve with nonvanishing position in the plane: a differentiable angle function with prescribed derivative).

[L5]

Cross product rules: e1×e2=e3, e2×e1=−e3, the product is bilinear, and ei are unit vectors (The cross product in R3, The cross product is bilinear, alternating, and orthogonal to both factors); cos⁡2+sin⁡2=1 (Parity and the Pythagorean identity for sine and cosine).

Proof

technique · direct
1.1givenL1L2L4algebra

Put ρi:=⟨r,ei⟩ for i=1,2; these are twice differentiable by [L4]. They are never both zero: by [L1] and [L2], r=ρ1e1+ρ2e2+⟨r,e3⟩e3=ρ1e1+ρ2e2, so ρ1=ρ2=0 would give r=0, contradicting ∣r∣∈J0⊆(0,∞).

1.2givenL1L5algebra

Clause 1 for an arbitrary pair: if r=ρ(cos⁡θ e1+sin⁡θ e2), then ∣r∣2=∣r∣2 is computed by [L1] and [L5] as ρ2(cos⁡2θ+sin⁡2θ)=ρ2, hence ρ=∣r∣.

1.3givenL4L5algebra

Differentiating [L5]'s expansion r=ρ n with n:=cos⁡θ e1+sin⁡θ e2 and n⊥:=−sin⁡θ e1+cos⁡θ e2 by [L4] gives r′=ρ′n+ρθ′ n⊥,r′′=(ρ′′−ρθ′2)n+(2ρ′θ′+ρθ′′)n⊥. Here n×n⊥=e3 and n⊥×n=−e3 by bilinearity of the cross product and e1×e2=e3.

2.1step 1.1L1L2L3algebra

Apply the polar lift [L3] to (ρ1,ρ2): there is a twice differentiable θ with ρ1=r0cos⁡θ, ρ2=r0sin⁡θ for r0=ρ12+ρ22; by Parseval in [L1] and ⟨r,e3⟩=0, r02=ρ12+ρ22=∣r∣2, so r0=∣r∣=:ρ is twice differentiable with values in J0⊆(0,∞). Substituting into the expansion of step 1.1 gives r=ρ(cos⁡θ e1+sin⁡θ e2), so such a pair exists.

2.2givenstep 1.3L5algebra

Clause 2: from step 1.3, r×r′=ρ n×(ρ′n+ρθ′n⊥)=ρ2θ′ e3, so j=Je3=μr×r′=μρ2θ′e3, giving J=μρ2θ′ and hence ρ2θ′=J/μ>0.

3.1givenstep 1.2step 1.3step 2.2L1L5algebra

Clause 3: the trajectory equation is μr′′=−(U′(∣r∣)/∣r∣)r=−U′(ρ) n by clause 1. Comparing with the decomposition of step 1.3 and using the linear independence of n and n⊥ — their components in the orthonormal triad are distinct, so αn+βn⊥=0 forces α=β=0 by taking inner products with n and n⊥ — the radial component gives μ(ρ′′−ρθ′2)=−U′(ρ). Substituting ρ2θ′=J/μ from step 2.2 gives μρ′′=−U′(ρ)+μρθ′2=−U′(ρ)+J2μρ3=−(UJeff)′(ρ), because (UJeff)′=U′(ρ)−J2μρ3.

3.2step 1.3step 2.2L1algebra

Clause 4, the energy identity: from step 1.3 and orthonormality, ∥r′∥22=ρ′2+ρ2θ′2, so E=12μρ′2+12μρ2θ′2+U(ρ)=12μρ′2+J22μρ2+U(ρ)=12μρ′2+UJeff(ρ) by step 2.2.

4.1step 3.1step 3.2L4L6algebra∎

The right-hand side is differentiable on I and its derivative is μρ′ρ′′+(UJeff)′(ρ) ρ′=ρ′(μρ′′+(UJeff)′(ρ))=0 by step 3.1 (the equation holds at every point of I); hence E is constant on I by [L6].

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Circular orbits of a non-radial central-force trajectory: the critical-point condition and the squared frequency of the linearised radial equation

Statement

Let μ>0, let J0⊆(0,∞) be an open interval, let U:J0→R be of class C2, let I be an interval with at least two points, and let r:I→R3 be a central-force trajectory for (μ,U) with constant nonzero angular momentum j=Je3, J>0, and polar data (ρ,θ) as in Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation, so that ρ2θ′=J/μ and μρ′′=−(UJeff)′(ρ) on I. Fix r0∈J0. Then:

  1. (Necessity.) If ρ≡r0 is constant on I, then (UJeff)′(r0)=0 and θ′≡J/(μr02) is constant; the trajectory is the uniform rotation r(t)=r0(cos⁡θ(t) e1+sin⁡θ(t) e2).
  2. (Existence.) Conversely, if (UJeff)′(r0)=0, then for every θ0∈R the curve r∘(t):=r0(cos⁡(θ0+ωt) e1+sin⁡(θ0+ωt) e2),ω:=Jμr02, is a central-force trajectory for (μ,U) defined on R, with angular momentum j=Je3 and constant energy E=12μr02ω2+U(r0)=UJeff(r0).
  3. (Linearised radial equation.) Let ρ be a solution of the reduced equation μρ′′=−(UJeff)′(ρ) on I with (UJeff)′(r0)=0, and put δ:=ρ−r0. Then μδ′′=−(UJeff)′′(r0) δ+g(δ),g(x):=(UJeff)′′(r0)x−(UJeff)′(r0+x), where g(x)/x→0 as x→0; if moreover (UJeff)′′(r0)>0, then the associated linearised equation μδ′′=−(UJeff)′′(r0)δ has exactly the solutions δ(t)=Acos⁡(ωlin(t−t0))+Bsin⁡(ωlin(t−t0)) with A,B∈R and ωlin:=(UJeff)′′(r0)/μ>0, for any t0∈I. This is a statement about the quadratic radial approximation; no nonlinear stability of the circular orbit is asserted.

Facts & Assumptions

Given: The central-force data with constant nonzero angular momentum j=Je3 and polar data (ρ,θ) from Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation, a radius r0∈J0, and the functions n(θ):=cos⁡θ e1+sin⁡θ e2, n⊥(θ):=−sin⁡θ e1+cos⁡θ e2.

[L1]

Polar reduction: r=ρ n(θ) with ρ=∣r∣, ρ2θ′=J/μ, μρ′′=−(UJeff)′(ρ), and the energy is constant (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation).

[L3]

sin⁡2+cos⁡2=1; and for the linearised equation: the solutions of v′′+ω2v=0 on an interval with ω>0 are exactly Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)) (Parity and the Pythagorean identity for sine and cosine, Solutions of v′′+ω2v=0 are exactly the functions Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)), Ck Euclidean maps and diffeomorphisms).

[L4]

The derivative of (UJeff)′ at r0 is (UJeff)′′(r0): for every ε>0 there is η>0 with ∣(UJeff)′(r0+x)−(UJeff)′(r0)−(UJeff)′′(r0)x∣≤ε∣x∣ for ∣x∣<η (The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set, Limit point, isolated point, adherent point, derived set, and dense subset of R, Intervals of R: the nine order-convex forms, nondegeneracy, and length).

Proof

technique · direct
1.1givenL1algebra

Clause 1: if ρ≡r0 then ρ′′=0, so the reduced equation of [L1] gives (UJeff)′(r0)=0; and θ′=J/(μρ2)=J/(μr02) is constant by [L1].

1.2givenL2L3algebraconstruct

Clause 2, the kinematics: for r∘(t)=r0 n(θ0+ωt) with ω=J/(μr02), [L2] gives (r∘)′=r0ω n⊥(θ0+ωt) and (r∘)′′=−r0ω2 n(θ0+ωt); hence j=μ r∘×(r∘)′=μr02ω e3=Je3 and ∥(r∘)′∥22=r02ω2, using [L2] and [L3].

1.3givenL1L4algebra

Clause 3, the remainder: by [L4] applied to (UJeff)′ at r0, for every ε>0 there is η>0 with ∣g(x)∣≤ε∣x∣ for ∣x∣<η, where g(x)=(UJeff)′′(r0)x−(UJeff)′(r0+x) and (UJeff)′(r0)=0; that is, g(x)/x→0. Substituting ρ=r0+δ into the reduced equation of [L1] and expanding gives μδ′′=−(UJeff)′(r0+δ)=−(UJeff)′′(r0)δ+g(δ).

2.1givenstep 1.2L1algebra

Clause 2, the equation of motion: from step 1.2 and ω2=J2/(μ2r04), μ(r∘)′′=−μr0ω2n=−J2μr03 n, while −(U′(r0)/r0)r∘=−U′(r0) n; the two agree because (UJeff)′(r0)=0 means U′(r0)=J2/(μr03). The energy is 12μr02ω2+U(r0)=J22μr02+U(r0)=UJeff(r0), and it is constant because the equation is satisfied and the trajectory is periodic; more directly it equals the constant right-hand side of clause 4 of [L1] along this motion.

3.1step 1.3L3algebra∎

Clause 3, the linearised solutions: if (UJeff)′′(r0)>0 then ωlin2=(UJeff)′′(r0)/μ>0 and the linearised equation reads δ′′+ωlin2δ=0; by [L3] its solutions on I are exactly the displayed trigonometric functions for any t0∈I.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Binet's equation for the attractive inverse-square field: d2u/dφ2+u=μκ/J2 for u=1/ρ

Statement

Let μ>0, κ>0, let U:(0,∞)→R be U(ρ):=−κ/ρ (of class C2), let I be an interval with at least two points, and let r:I→R3 be a central-force trajectory for (μ,U) with constant angular momentum j=Je3, J>0. Let (ρ,θ) be polar data as in Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation, let I∘ be the interior of I (a nonempty interval with at least two points, or empty if I is a singleton, which is excluded), and put W:=θ(I∘), an interval. Then:

  1. θ:I∘→W is a C1 bijection with inverse σ:W→I∘ of class C1, and the function u:=1/(ρ∘σ):W→(0,∞) is of class C2;
  2. on W, u′′+u=μκ/J2 (derivatives with respect to the variable of W);
  3. the auxiliary identity ρ′=−Jμ u′∘θ holds on I∘;
  4. the energy satisfies E=J22μ((u′)2+u2)−κu on W, where E is the constant energy of the trajectory.

Facts & Assumptions

Given: The inverse-square central-force trajectory with j=Je3, J>0, and polar data (ρ,θ), the interval W=θ(I∘), and the associated functions σ and u.

[L1]

Polar reduction: on all of I, r=ρ n(θ), ρ2θ′=J/μ (so θ′>0 and θ is strictly increasing), μρ′′=−(UJeff)′(ρ) with UJeff(ρ)=−κ/ρ+J2/(2μρ2), and the energy E is constant (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation).

Proof

technique · direct
1.1givenL1L2L3algebra

By [L1], ρ′′=−(UJeff)′(ρ)/μ on I∘. Since UJeff is C2 and ρ is continuous, the right-hand side is continuous, so ρ is C2 there. Thus θ′=J/(μρ2) is C1 and positive. By [L2], θ:I∘→W has a C1 inverse σ with σ′=1/(θ′∘σ); the right-hand side is C1 by [L3], so σ is C2. Therefore u:=1/(ρ∘σ) is a quotient of C2 functions with nonvanishing denominator and is C2 on W.

2.1step 1.1L1L3

Chain rule along the trajectory: for t∈I∘, u(θ(t))=1/ρ(t), so [L3] gives u′(θ(t)) θ′(t)=−ρ′(t)/ρ(t)2; multiplying by ρ2 and using ρ2θ′=J/μ from [L1], ρ′(t)=−ρ(t)2θ′(t)u′(θ(t))=−Jμ u′(θ(t)), which is clause 3.

3.1step 1.1step 2.1L1L3algebra

Differentiating clause 3 with [L3] and using θ′=J/(μρ2) again gives ρ′′(t)=−Jμu′′(θ(t)) θ′(t)=−J2μ2ρ(t)2 u′′(θ(t))(t∈I∘). Since W=θ(I∘), every φ∈W equals θ(t) for some t∈I∘, so this identity supplies the corresponding value of u′′ at every point of W.

4.1step 1.1step 3.1L1algebra

The reduction equation of [L1] for U(ρ)=−κ/ρ reads μρ′′=−κ/ρ2+J2/(μρ3). Substituting steps 2.1 and 3.1, with ρ=1/u, gives −J2μρ2u′′=−κρ2+J2μρ3  ⟹  u′′=μκJ2−u, that is, u′′+u=μκ/J2 on W, which is clause 2.

5.1step 1.1step 2.1L1L3algebra∎

Clause 4: the energy identity of [L1] gives E=12μρ′2+UJeff(ρ)=12μρ′2−κρ+J22μρ2. Substituting ρ=1/u and clause 3, 12μρ′2=J22μ(u′∘θ)2 and J22μρ2=J22μu2; hence E=J22μ((u′∘θ)2+u(θ)2)−κ u(θ) along the trajectory. Since θ ranges over W, the displayed identity between the constant E and the continuous function J22μ((u′)2+u2)−κu on W holds on W.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

The Kepler orbit equation is a conic: u=c(1+ecos⁡(φ−θ0)) with E=(μκ2/2J2)(e2−1)

Statement

Let μ>0, κ>0, U(ρ)=−κ/ρ on (0,∞), let r:I→R3 be a central-force trajectory for (μ,U) with constant angular momentum j=Je3, J>0 and constant energy E (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential), and let u:W→(0,∞), θ and σ be as in Binet's equation for the attractive inverse-square field: d2u/dφ2+u=μκ/J2 for u=1/ρ, with W=θ(I∘). Write c:=μκ/J2>0 and p:=1/c=J2/(μκ). Then:

  1. there are e≥0 and θ0∈R with u(φ)=c(1+ecos⁡(φ−θ0))(φ∈W), and for t∈I∘ accordingly ρ(t)=p/(1+ecos⁡(θ(t)−θ0));
  2. the energy and the eccentricity are related by E=μκ22J2(e2−1); consequently E<0  ⟺  e<1, E=0  ⟺  e=1 and E>0  ⟺  e>1.

Facts & Assumptions

Given: The inverse-square central-force trajectory with constant angular momentum j=Je3, J>0, constant energy E, the solution u of Binet's equation on the interval W, and the constants c=μκ/J2, p=1/c.

[L1]

Binet's equation and the auxiliary identities: u′′+u=μκ/J2=c on W, W is an interval, and E=J22μ((u′)2+u2)−κu on W (Binet's equation for the attractive inverse-square field: d2u/dφ2+u=μκ/J2 for u=1/ρ).

[L2]

Solutions of v′′+v=0 on an interval with at least two points are Acos⁡(φ−φ0)+Bsin⁡(φ−φ0) for fixed φ0 in that interval (Solutions of v′′+ω2v=0 are exactly the functions Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)), Intervals of R: the nine order-convex forms, nondegeneracy, and length).

[L3]

Subtraction formula: cos⁡(α−ψ)=cos⁡αcos⁡ψ+sin⁡αsin⁡ψ; (t↦(cos⁡t,sin⁡t)) is a bijection [0,2π)→S1; sin⁡2+cos⁡2=1 (The subtraction formulas for sine and cosine, t↦(cos⁡t,sin⁡t) is a bijection from [0,2π) onto the real unit circle, Parity and the Pythagorean identity for sine and cosine).

Proof

technique · direct
1.1L1L2algebra

Note u′′+u=c is equivalent to (u−c)′′+(u−c)=0 on W by linearity of the second derivative. Since W is an interval and Binet's equation holds on it (it contains at least two points because I∘ has at least two points), [L2] gives reals A,B with u(φ)−c=Acos⁡(φ−φ0)+Bsin⁡(φ−φ0)(φ∈W) for a fixed φ0∈W.

2.1step 1.1L3algebra

Put R:=A2+B2≥0. If R>0, then (A/R,B/R) is a unit vector, so by [L3] there is ψ∈[0,2π) with cos⁡ψ=A/R, sin⁡ψ=B/R; if R=0 take ψ:=0. Then for every φ, [L3] gives Acos⁡(φ−φ0)+Bsin⁡(φ−φ0)=Rcos⁡(φ−φ0−ψ). Setting θ0:=φ0+ψ and e:=R/c=R J2/(μκ)≥0, step 1.1 becomes u(φ)=c+e ccos⁡(φ−θ0)=c(1+ecos⁡(φ−θ0)) on W, which is the first assertion; the polar form ρ=p/(1+ecos⁡(θ−θ0)) follows by inverting u=1/ρ on I∘.

3.1step 2.1L1L4algebra

Energy: by [L1] and step 2.1, with ψ(φ):=φ−θ0, u=c(1+ecos⁡ψ) and u′=−cesin⁡ψ (differentiate with [L4]); hence (u′)2+u2=c2(e2sin⁡2ψ+1+2ecos⁡ψ+e2cos⁡2ψ)=c2(1+2ecos⁡ψ+e2), and κu=κc(1+ecos⁡ψ). Substituting into [L1]: E=J2c22μ(1+2ecos⁡ψ+e2)−κc (1+ecos⁡ψ). With c=μκ/J2 one has J2c22μ=μκ22J2 and κc=μκ2J2, so E=μκ22J2[(1+2ecos⁡ψ+e2)−2(1+ecos⁡ψ)]=μκ22J2(e2−1).

4.1step 3.1algebra∎

Classification: since μκ2/(2J2)>0, the sign of E equals the sign of e2−1; with e≥0 this is negative exactly for e<1, zero exactly for e=1 and positive exactly for e>1.

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

The one-dimensional action variable of a family of two-turning-point energy levels

Definition

Let J0⊆R be an open interval, let U:J0→R be continuous (Continuity of f:A→R at a point of A and on A: the ε-δ condition, its agreement with lim⁡x→cf(x)=f(c) at a limit point, and continuity at an isolated point), let m>0, let E0<E1 be reals, and let a,b:(E0,E1)→J0 be functions with

a(E)<b(E),U(a(E))=U(b(E))=E,U(q)<E (a(E)<q<b(E))(E0<E<E1).

The action variable of this datum is the function Ja:(E0,E1)→R,

Ja(E):=1π∫a(E)b(E)2m(E−U(q))  dq,

the integral being the Riemann integral on the compact interval [a(E),b(E)] (The lower and upper Darboux integrals of a bounded f on [a,b] as sup⁡PL(f,P) and inf⁡PU(f,P), Darboux integrability as their equality, and the notation ∫abf). It is defined because q↦2m(E−U(q)) is continuous on [a(E),b(E)]: the radicand is a continuous nonnegative function there (it vanishes exactly at the two endpoints). Existence and uniqueness of nonnegative roots follow from Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}. For x,y≥0, if s=x≥t=y, then (s−t)2≤(s−t)(s+t)=x−y; interchanging x,y gives ∣x−y∣2≤∣x−y∣. Thus ∣x−y∣<ε2 implies ∣x−y∣<ε, proving continuity, including at zero. The continuous integrand is Riemann integrable by A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion.

Reading of the notation. The interval [a(E),b(E)] is the classically allowed segment between the two turning points of the energy E, and the closed curve p=±2m(E−U(q)) encloses the area A(E)=2∫a(E)b(E)2m(E−U(q)) dq (the upper graph minus the lower graph), so Ja(E)=A(E)/(2π); no physical interpretation of m, U or E is a mathematical premise here, and no period is claimed by this definition.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

The subcritical Kepler conic is an ellipse with semiaxes a=p/(1−e2), b=a1−e2, focus at the origin and area πab; its area is also 12∫02πϱ2

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let p>0 and 0≤e<1, and put a:=p1−e2>0,b:=a1−e2,ϱ(φ):=p1+ecos⁡φ(φ∈R).

  1. ϱ is C1, 2π-periodic, and p1+e≤ϱ≤p1−e.
  2. Let Q(x,y):=(1−e2)x2+y2+2pex−p2 and Eell:={(x,y)∈R2:Q(x,y)≤0}. Then Eell is Lebesgue measurable, and for the closed unit disc B‾⊆R2 and T(u,v):=(au−ae,bv) one has Eell=T[B‾]; consequently λ2(Eell)=πab=πa21−e2.
  3. For every φ, the ray {s(cos⁡φ,sin⁡φ):s>0} meets ∂Eell in the unique point ϱ(φ)(cos⁡φ,sin⁡φ); that point satisfies Q=0, and λ2(Eell)=12∫02πϱ(φ)2 dφ.

Thus the conic Q=0 is the ellipse with centre (−ae,0), semiaxes a,b, its right focus at the origin, and enclosed area πab.

Facts & Assumptions

Given: The constants p>0, 0≤e<1, a=p/(1−e2), b=a1−e2, the profile ϱ, and the region Eell.

[L1]

Substitution and completion of the square: Q(T(u,v))=Q(au−ae,bv)=a2(1−e2)(u2+v2−1) for all (u,v)∈R2; equivalently Q(x,y)=(1−e2)(x+pe/(1−e2))2+y2−p2/(1−e2). [algebra]

[L2]

Assuming countable choice, a C1 diffeomorphism satisfies the change-of-variables formula for L1 functions; the polar map Φ(φ,s)=(scos⁡φ,ssin⁡φ) is a C1 diffeomorphism of (0,2π)×(0,∞) onto its image (injective there, with nonsingular derivative, hence a diffeomorphism onto its open image by the injective-regular theorem); and a continuous function on a product with a nonnegative area element may be integrated iteratively by Tonelli (A C^1 diffeomorphism satisfies the change-of-variables formula for L^1 functions, An injective regular C1 map is a diffeomorphism onto its image, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, The Axiom of Countable Choice (ACω)).

Proof

technique · direct
1.1givenL4algebra

Clause 1: ϱ is a quotient of C1 functions with denominator 1+ecos⁡φ≥1−e>0, hence C1 by [L4]; it is 2π-periodic because cos⁡ is; and p/(1+e)≤ϱ≤p/(1−e) because cos⁡φ∈[−1,1].

1.2givenL1algebra

Clause 2, the identity Eell=T[B‾]: by [L1], for (u,v)∈R2, Q(T(u,v))=a2(1−e2)(u2+v2−1), so Q(T(u,v))≤0 if and only if u2+v2≤1, that is, (u,v)∈B‾; hence (x,y)∈Eell exactly when (x,y)=T(u,v) for some (u,v)∈B‾, since T is a bijection of R2 onto itself (its linear part is invertible, det⁡diag(a,b)=ab>0).

1.3givenL1L4algebra

Clause 3, the radial description: for s>0 put (x,y):=s(cos⁡φ,sin⁡φ). Using [L4], Q(x,y)=s2(1−e2cos⁡2φ)+2pescos⁡φ−p2, a quadratic in s with positive leading coefficient (since 1−e2cos⁡2φ≥1−e2>0) and negative constant term −p2<0; hence the set of s>0 with Q≤0 is the interval (0,ϱ(φ)]: indeed s=ϱ(φ) is the unique positive root of Q=0, because substituting ϱ=p/(1+ecos⁡φ) and multiplying by (1+ecos⁡φ)2 gives p2(1−e2cos⁡2φ)+2p2ecos⁡φ (1+ecos⁡φ)−p2(1+ecos⁡φ)2=0, a direct expansion — and Q(0)=−p2<0 while Q(s)→+∞ as s→∞, so the two roots have opposite signs and the positive one is ϱ(φ). Therefore Eell={s(cos⁡φ,sin⁡φ):0≤φ<2π, 0≤s≤ϱ(φ)}, the boundary is the displayed profile, and each ray meets it exactly once.

2.1step 1.2L2L3algebra

Clause 2, the area: T is an affine map with invertible linear part, hence a C1 diffeomorphism of R2 with ∣det⁡DT(u,v)∣=ab everywhere; Eell={Q≤0} is closed (so Lebesgue measurable with indicator in L1), and T[B‾]=Eell by step 1.2. The change-of-variables formula of [L2] applied to the indicator of Eell gives λ2(Eell)=∫B‾ab dλ2=ab λ2(B‾); the closed unit disc has λ2(B‾)=π by [L3], so λ2(Eell)=πab=πa21−e2.

3.1step 1.3step 2.1L2L3algebra∎

Clause 3, the polar area formula: let W:={(φ,s):0<φ<2π, 0<s≤ϱ(φ)}⊆D:=(0,2π)×(0,∞) and S:=Φ(W), where Φ(φ,s)=(scos⁡φ,ssin⁡φ) is the restriction of the polar map of [L2]. The set W is Borel because ϱ is continuous, and S is Borel because Φ is a homeomorphism from D onto its image. Since S⊆Eell by step 1.3 and Eell has finite measure by step 2.1, 1S is integrable. The change-of-variables formula [L2] applied on D gives λ2(S)=∫Ws dλ2(φ,s)=∫02π ⁣ ⁣∫0ϱ(φ)s ds dφ=12∫02πϱ(φ)2 dφ, the middle equality by Tonelli [L2] applied to the nonnegative function (φ,s)↦s1W. By step 1.3, Eell∖S is contained in the closed positive x-axis, including the origin, a subset of a line and hence null by [L3]. Therefore λ2(Eell)=λ2(S)=12∫02πϱ(φ)2 dφ, which is clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Completeness and periodicity of the negative-energy Kepler orbit, with T=2πμab/J and T2=4π2μa3/κ

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let μ>0, κ>0, U(ρ)=−κ/ρ on (0,∞), and let r:I→R3 be a central-force trajectory for (μ,U) with constant angular momentum j=Je3, J>0, and constant energy E<0. Let e∈[0,1), θ0∈R, p=J2/(μκ), a=p/(1−e2), b=a1−e2, and ϱ(φ):=p/(1+ecos⁡(φ−θ0)) be the conic data of The Kepler orbit equation is a conic: u=c(1+ecos⁡(φ−θ0)) with E=(μκ2/2J2)(e2−1), and let (ρ,θ) be polar data as in Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation. Then:

  1. (Completeness.) r extends uniquely to a central-force trajectory for (μ,U) defined on all of R; along it ρ(t)∈[p/(1+e), p/(1−e)] for every t, and ρ and ρ′ are bounded. Consequently every maximal central-force trajectory with these data has domain R.
  2. (Periodicity.) θ:R→R is a strictly increasing C1 bijection with Jμ(1−ep)2≤θ′≤Jμ(1+ep)2; if t0∈R and T>0 is determined by θ(t0+T)=θ(t0)+2π, then r(t+T)=r(t) and r′(t+T)=r′(t) for all t∈R.
  3. (Period value.) With rmax⁡:=p/(1−e) and rmin⁡:=p/(1+e), T=2πμabJ=2πμκ a3/2,T2=4π2μa3κ.

Facts & Assumptions

Given: The inverse-square central-force trajectory with constant nonzero angular momentum and negative energy, its conic data from The Kepler orbit equation is a conic: u=c(1+ecos⁡(φ−θ0)) with E=(μκ2/2J2)(e2−1), and its polar data.

[L1]

Polar reduction: r=ρ n(θ) with ρ=∣r∣, ρ2θ′=J/μ, μρ′′=−(UJeff)′(ρ) with UJeff(ρ)=−κ/ρ+J2/(2μρ2), so (UJeff)′(ρ)=κ/ρ2−J2/(μρ3) is C1 on (0,∞); the energy is constant (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation).

[L3]

Compact continuation: for z′=F(z) with F continuous and locally Lipschitz on an open U⊆Rn, a solution whose values near a finite endpoint lie in a compact subset of U extends past that endpoint; solutions agreeing at one time agree on every common interval containing it (Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it, Picard-Lindelöf local existence and uniqueness for first-order systems, First-order systems, initial value problems, and solutions on intervals, Local Lipschitz continuity in the state variable, locally uniform in time and parameters).

Proof

technique · direct
1.1givenL1L2L4algebra

By [L1] and [L2], for t in the interior of I the function φ(t):=θ(t) satisfies ρ(t)=ϱ(φ(t))∈[rmin⁡,rmax⁡] and, by [L4], ρ′(t)=ϱ′(φ(t))θ′(t) with ∣θ′∣=J/(μρ2)≤J/(μrmin⁡2); since ϱ′ is continuous and 2π-periodic, it is bounded on R by M1:=max⁡∣ϱ′∣, so ∣ρ′∣≤M:=M1J/(μrmin⁡2) on the interior of I; by continuity the bounds for ρ and ρ′ hold on all of I.

1.2givenL1L3algebra

Consider the state z:=(ρ,ρ′,θ) in the open set U0:=(0,∞)×R×R with right-hand side F(x):=(x2, −μ−1(UJeff)′(x1), J/(μx12)); the radial equation of [L1] and θ′=J/(μρ2) give z′=F(z). Here (UJeff)′(s)=κ/s2−J2/(μs3), so F is locally Lipschitz: on a state box with x1≥δ>0, each x1−k for k=2,3 satisfies ∣x1−k−y1−k∣≤kδ−k−1∣x1−y1∣ by factoring the difference of reciprocal powers, while the first component is linear. Along every solution of this system, H(ρ,ρ′):=12μρ′2+UJeff(ρ) is constant, since H′=ρ′(μρ′′+(UJeff)′(ρ))=0; its value is the given E<0. With p=J2/(μκ), direct algebra gives E−UJeff(s)=κ(1−e2)2ps2(s−p1+e)(p1−e−s). Thus every such solution has rmin⁡=p/(1+e)≤ρ≤rmax⁡=p/(1−e); its energy identity bounds ρ′, and θ′=J/(μρ2) is bounded as well. On a bounded time interval, θ is consequently bounded, so near every finite endpoint z lies in a compact subset of U0. Compact continuation [L3] extends it past every finite right endpoint. For a finite left endpoint apply [L3] to z^(s):=z(−s), which solves the time-reversed locally Lipschitz system z^′=−F(z^). Hence the maximal extension of z is defined on all of R. For n(θ):=cos⁡θ e1+sin⁡θ e2 and n⊥(θ):=−sin⁡θ e1+cos⁡θ e2, differentiating r=ρn(θ) gives r′′=(ρ′′−ρθ′2)n(θ)+(2ρ′θ′+ρθ′′)n⊥(θ). Since θ′=J/(μρ2), θ′′=−2ρ′θ′/ρ and the second coefficient vanishes; the radial equation then gives μ(ρ′′−ρθ′2)=−U′(ρ). Thus the reconstructed r satisfies the central-force equation. For any other continuation, the product rule and the central-force equation give j′=μ(r′×r′+r×r′′)=0, so its angular momentum vector agrees with the one on the overlap; cross-product orthogonality places it in the same plane. By [L1] it has polar data there, and shifting its angle lift by an integer multiple of 2π makes the angle agree at a common time. The resulting state solves the same locally Lipschitz system with the same initial state, so uniqueness [L3] gives the same ρ,θ and hence the same trajectory. This is clause 1.

2.1step 1.1step 1.2L1L2algebra

On R, ρ∈[rmin⁡,rmax⁡] and θ′=J/(μρ2)∈[J(1−e)2/(μp2),J(1+e)2/(μp2)], since rmin⁡=p/(1+e) and rmax⁡=p/(1−e). Applying the conic supplier [L2] to the global extension gives the same profile ϱ: its sine and cosine coefficients agree with the original ones on the nonempty open angle interval θ(I∘), so their differences vanish identically (evaluate the difference and its derivative at one point). Thus ρ=ϱ∘θ on R. Hence θ is strictly increasing with θ′(t)≥c1>0, so θ(t)→±∞ as t→±∞ and θ(R)=R.

3.1step 2.1L4L6algebra

Let t0∈R and let T>0 satisfy φ(t0+T)=φ(t0)+2π, i.e. θ(t0+T)=θ(t0)+2π (existence and uniqueness of T come from step 2.1). Then ρ(t0+T)=ϱ(φ(t0)+2π)=ϱ(φ(t0))=ρ(t0) by the 2π-periodicity of ϱ, and ρ′(t0+T)=ϱ′(φ(t0)+2π)θ′(t0+T)=ϱ′(φ(t0)) Jμρ(t0+T)2=ϱ′(φ(t0)) Jμρ(t0)2=ρ′(t0), using [L4], the 2π-periodicity of ϱ′, and step 2.1. Also n(θ(t0+T))=n(θ(t0)) by the periodicity of sine and cosine [L6], so r(t0+T)=r(t0) and r′(t0+T)=r′(t0) (differentiate r=ρ n(θ) and use the equalities just proved).

4.1step 2.1step 3.1L2L4L5algebra

Change of variables in the defining relation of T: θ:[t0,t0+T]→[θ(t0),θ(t0)+2π] is a continuous bijection with θ′>0, so its inverse θ−1 is differentiable with (θ−1)′=1/θ′ composed with θ−1 [L5]. Hence, using θ′=J/(μρ2) and ρ=ϱ(θ) on the interval, T=∫t0t0+Tdt=∫θ(t0)θ(t0)+2π1θ′ dθ=μJ∫θ(t0)θ(t0)+2πϱ(θ)2 dθ=μJ∫02πϱ(ψ)2 dψ, the last equality by the 2π-periodicity of ϱ.

4.2step 1.2step 3.1L3L6algebra

Define zT(t):=(ρ(t+T),ρ′(t+T),θ(t+T)−2π). Since the state field in step 1.2 is independent of θ, zT solves the same equation as z(t):=(ρ(t),ρ′(t),θ(t)). By step 3.1 and the definition of T, zT(t0)=z(t0); uniqueness [L3] gives zT=z on R. Periodicity of sine and cosine then yields r(t+T)=r(t), and differentiation yields r′(t+T)=r′(t). This is clause 2.

5.1step 4.1L2algebra

By clause 3 of The subcritical Kepler conic is an ellipse with semiaxes a=p/(1−e2), b=a1−e2, focus at the origin and area πab; its area is also 12∫02πϱ2, applied to the unshifted profile ϱ~(ψ):=p/(1+ecos⁡ψ), with a translation of the integration variable by θ0 and invariance of a full-period integral under translation, 12∫02πϱ(ψ)2 dψ=λ2(Eell)=πab, so step 4.1 gives T=2πμabJ.

6.1step 5.1algebra∎

Algebra: p=J2/(μκ) and p=a(1−e2) by the definitions, so J2=μκa(1−e2) and J=μκa(1−e2); with b=a1−e2, 2πμabJ=2πμa21−e2μκa(1−e2)=2πμκ a3/2, and squaring gives T2=4π2μa3/κ. This is clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

The action variable of a two-simple-turning-point family is C1 with derivative T(E)/(2π), the period quadrature over 2π

Statement

Let J0⊆R be an open interval, let U:J0→R be of class C3 (Ck Euclidean maps and diffeomorphisms), let m>0, let E0<E1, and let a,b:(E0,E1)→J0 be functions with a(E)<b(E), U(a(E))=U(b(E))=E, U(q)<E for a(E)<q<b(E), and

U′(a(E))<0<U′(b(E))(E0<E<E1)

(the energy E has exactly the two simple turning points a(E), b(E) in J0). Then:

  1. a and b are of class C1 on (E0,E1);
  2. the action variable Ja of The one-dimensional action variable of a family of two-turning-point energy levels is of class C1 on (E0,E1), and for every E∈(E0,E1) Ja′(E)=12πT(E),T(E):=2m∫a(E)b(E)dqE−U(q)∈(0,∞), where the integral is the improper integral over the open interval (a(E),b(E)), convergent at both endpoints.

Facts & Assumptions

Given: The data J0,U,m,E0,E1,a,b of the statement, with U of class C3, exactly two simple turning points, and U<E in between.

[L1]

The Euclidean inverse function theorem: if F is C1 near c and DF(c) is invertible, then F restricts to a C1 bijection from a neighbourhood of c onto a neighbourhood of F(c) with C1 inverse (The Euclidean inverse function theorem). A Ck map with invertible derivative has a Ck local inverse by the same theorem applied to the derivative and induction on smoothness; in the one-variable case U′(ξ)≠0 makes ξ↦U(ξ) a C3 diffeomorphism near ξ onto its image.

[L2]

Leibniz's rule on a compact rectangle: for continuous g,h on [a,b]×[c,d] with x↦g(x,t) differentiable on (a,b) with derivative h(x,t), the function x↦∫cdg(x,t) dt is differentiable with derivative ∫cdh(x,t) dt; with the chain rule this gives the variable-limit rule ddE∫cz(E)G(E,z) dz=G(E,z(E))z′(E)+∫cz(E)∂EG(E,z) dz (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c)).

[L3]

Substitution: for φ differentiable on [c,d] with integrable φ′ and f continuous on an interval containing φ([c,d]), ∫φ(c)φ(d)f=∫cd(f∘φ)φ′ (Substitution: if φ is differentiable on [c,d] with φ′ integrable and f is continuous on an interval containing φ([c,d]), then ∫φ(c)φ(d)f=∫cd(f∘φ) φ′).

Proof

technique · direct
1.1givenL1choosealgebra

Fix E∗∈(E0,E1). Since U′(a(E∗))≠0 and U′(b(E∗))≠0, the inverse function theorem [L1] applied to U at a(E∗) and at b(E∗) gives δ>0 and C3 (hence C1) functions A,B:(E∗−δ,E∗+δ)→J0 with U(A(E))=U(B(E))=E, A(E∗)=a(E∗) and B(E∗)=b(E∗). Shrinking δ, continuity gives A<B on the smaller interval. For every such E, the two real numbers A(E) and B(E) are distinct solutions of U(q)=E; by hypothesis the only two solutions in J0 are a(E)<b(E), so {A(E),B(E)}={a(E),b(E)}; since A<B and a<b, this forces a=A and b=B on that neighbourhood. Hence a and b are C1 near the arbitrary point E∗, which is clause 1.

2.1givenstep 1.1constructalgebra

Fix E∗ again. Since U′ is continuous and U′(a(E∗))<0<U′(b(E∗)), choose disjoint neighborhoods of a(E∗) and b(E∗) on which U′<0 and U′>0, respectively. Choose q1,q2 in these neighborhoods with a(E∗)<q1<q2<b(E∗). By continuity of a and b from step 1.1, shrink δ>0 so that for ∣E−E∗∣<δ, a(E) and q1 remain in the first neighborhood, q2 and b(E) remain in the second, and a(E)<q1<q2<b(E). Thus U is strictly decreasing on [a(E),q1] and strictly increasing on [q2,b(E)]. On this energy interval write M(E):=∫q1q22m(E−U(q)) dq,L(E):=∫a(E)q12m(E−U(q)) dq,R(E):=∫q2b(E)2m(E−U(q)) dq, so that Ja=(L+M+R)/π there.

3.1step 2.1L2L4L5algebra

The function u(E,q):=2m(E−U(q)) is continuous on the compact rectangle [E∗−δ,E∗+δ]×[q1,q2] and ∂Eu(E,q)=m/2m(E−U(q)) is continuous there, because U(q)<E∗ for every q∈[q1,q2] and U is continuous on the compact interval, so that E−U(q)≥ϵ>0 on the rectangle after shrinking δ [L4, L5]. Hence by [L2] the middle term is C1 with M′(E)=∫q1q2m dq2m(E−U(q))(∣E−E∗∣<δ).

3.2step 2.1L1L3algebraconstruct

For the left leg, use the open interval about a(E∗) chosen in step 2.1, containing every segment [a(E),q1] and having U′<0. On this interval U is strictly monotone; its local C3 inverses from [L1] agree on overlaps by injectivity, giving a C3 inverse ψ on the whole image interval. For E near E∗ the substitution s=U(q) [L3] gives L(E)=∫U(q1)E2m(E−s) w(s) ds,w(s):=−1U′(ψ(s))>0, and w is C2 on the image interval; here U(q1)<E for E near E∗ since U(q1)<E∗.

3.3step 2.1L1L2L3L5algebra

The right leg is treated identically: near b(E∗) one has U′>0, the C3 local inverse ψ of U gives R(E)=∫U(q2)E2m(E−s) w(s) ds with w:=1/U′(ψ(s))>0 of class C2, the substitution s=E−z2 and the variable-limit Leibniz rule give R′(E)=∫q2b(E)m dq2m(E−U(q)).

4.1step 3.2L2L3L5algebra

Substituting s=E−z2 in the last integral [L3] turns the left leg into L(E)=22m∫0z(E)z2w(E−z2) dz,z(E):=E−U(q1)>0, with z of class C1 and the integrand G(E,z):=22mz2w(E−z2) of class C1 on a compact rectangle. By the variable-limit Leibniz rule [L2], L′(E)=G(E,z(E))z′(E)+∫0z(E)∂EG(E,z) dz=2m z(E)w(E−z(E)2)+22m∫0z(E)z2w′(E−z2) dz.

5.1step 4.1L3L5algebra

Integration by parts, equivalently the product rule applied to z↦z w(E−z2), gives ∫0z(E)w(E−z2) dz=z(E) w(E−z(E)2)+2∫0z(E)z2w′(E−z2) dz; substituting this into step 4.1 yields L′(E)=2m∫0z(E)w(E−z2) dz=∫a(E)q1m dq2m(E−U(q)), the last equality by the same substitutions read backwards [L3]: 2m∫0z(E)w(E−z2) dz=∫a(E)q1u(E,q)−1m dq.

6.1step 3.1step 5.1step 3.3L3L4algebra

Adding the three derivatives of steps 3.1, 5.1 and 3.3 gives Ja′(E)=1π∫a(E)b(E)m dq2m(E−U(q))=2m2π∫a(E)b(E)dqE−U(q)=T(E)2π for ∣E−E∗∣<δ, with the middle integral understood improperly at a(E) and b(E). The endpoint pieces converge by steps 3.2–5.1 and 3.3; additivity on interior truncations and passage to the endpoint limits justify splitting this improper integral at q1,q2. The second equality uses m/2m=2m/2, and the last is the definition of T(E).

7.1step 6.1L4algebra∎

The integrand q↦1/E−U(q) is continuous and strictly positive on the nonempty open interval (a(E),b(E)), so T(E)>0 by [L4]; and the derivative computed in step 6.1 is independent of the auxiliary E∗, so the formula holds on all of (E0,E1) and Ja is C1 there. This is clause 2.

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

Scattering data: the incident flux measure, the solid-angle element, deflection branches and the cross-section density

Definition

Assume ACω (The Axiom of Countable Choice (ACω)) for the Borel/Lebesgue measure foundation below.

Data. Let μ>0, let U:(0,∞)→R be of class C2 (Ck Euclidean maps and diffeomorphisms), let R>0 satisfy U(r)=0 for every r≥R (a finite-range potential), and let E>0 and J>0 be reals. Write

V(r):=U(r)+J22μr2(r>0).

Incident flux measure. Let ν be the Borel measure on the open rectangle (0,∞)×(0,2π) given by the density b against Lebesgue measure λ2,

ν(A):=∫Ab dλ2(b,φ)(A⊆(0,∞)×(0,2π) Borel),

the integral being that of Integral over a measurable subset; the nonnegative Borel function (b,φ)↦b defines a measure by The indefinite integral of a nonnegative measurable function is a measure (Borel measurable and Lebesgue measurable functions on Rn). In the model this measure records a uniform incident beam: an annulus of impact parameters [b1,b2) contributes π(b22−b12) over the full azimuth range, or 12(b22−b12) per unit azimuth.

Solid-angle element. Let Ω be the Borel measure on (0,π)×(0,2π) with density sin⁡ against λ2,

Ω(B):=∫Bsin⁡θ dλ2(θ,φ)(B⊆(0,π)×(0,2π) Borel),

where sin⁡ is the function of Sine and cosine defined by their real power series. It is continuous, being differentiable by The derivatives of sine and cosine are cosine and minus sine. Also sin⁡0=sin⁡π=0 by The zero sets of sine and cosine and the least positive common period 2 pi; strict increase on [0,π/2] and strict decrease on [π/2,π] (Signs, monotonicity intervals, and ranges of sine and cosine) give sin⁡θ>0 for 0<θ<π. Thus (θ,φ)↦sin⁡θ is a nonnegative Borel function (continuous preimages of Borel sets are Borel, A continuous map has Borel preimages of Borel sets), and The indefinite integral of a nonnegative measurable function is a measure establishes that Ω is a measure. This is the standard solid-angle element written sin⁡θ dθ dφ (the surface measure of the unit sphere in spherical coordinates, restricted to the open chart).

Deflection branch. A deflection branch is an interval (a1,a2)⊆(0,∞) with at least two points together with a function θ:(a1,a2)→(0,π) that is of class C1 and satisfies θ′(b)≠0 for every b∈(a1,a2). Such a θ is strictly monotone (the continuous nonvanishing derivative has constant sign, and The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a) gives strict monotonicity) and its image is an interval with endpoints θ1<θ2 (Intervals of R: the nine order-convex forms, nondegeneracy, and length); its inverse b:(θ1,θ2)→(a1,a2) is differentiable with b′(θ)=1/θ′(b(θ)), by the derivative-of-an-inverse theorem Derivative of an inverse: if f is continuous and injective on a nondegenerate interval I and differentiable at c∈I with f′(c)≠0, then the inverse g is differentiable at f(c) with g′(f(c))=1/f′(c); and if f′(c)=0 then g is not differentiable at f(c) applied to the continuous injective θ (every point of (a1,a2) is a limit point of it, Limit point, isolated point, adherent point, derived set, and dense subset of R, The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set).

Cross-section density. The cross-section density of the branch is the function

dσdΩ(θ):=b(θ)sin⁡θ ∣θ′(b(θ))∣(θ1<θ<θ2).

Mathematical status. The measures, branches and density are abstract data independent of the finite-range trajectory data (μ,U,R,E,J). A branch is not asserted to arise from any trajectory; the scattering theorems of this page identify the deflection angle computed from the central-force equation with the values of such a function under their hypotheses, and they compute the pushforward of ν under the induced map. No counting of detected particles, and no empirical claim, is part of this definition.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Finite-range scattering: the simple outer turning point, the two global branches and the free tails

Statement

Let μ>0, let U:(0,∞)→R be of class C2, let R>0 satisfy U(r)=0 for every r≥R, let E>0 and J>0, and put V(r):=U(r)+J2/(2μr2) (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential). Suppose r∗>0 satisfies

V(r∗)=E,V′(r∗)<0,V(r)<E(r>r∗).

Let r:I→R3 be a central-force trajectory for (μ,U) with constant nonzero angular momentum j=Je3, J>0, energy E, and ∣r(t)∣≥r∗ for all t∈I, whose radial part attains the value r∗ at some t∗∈I (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation). Then:

  1. (ρ′)2=2μ(E−V(ρ)) on I, where ρ=∣r∣; ρ′(t∗)=0 and ρ′′(t∗)=−V′(r∗)/μ>0, so t∗ is a strict minimum of ρ;
  2. r extends uniquely to a central-force trajectory on all of R; along the extension ρ≥r∗, ρ(t)>ρ(t∗)=r∗ for t≠t∗, and ρ is strictly decreasing on (−∞,t∗] and strictly increasing on [t∗,∞);
  3. the time quadrature τ(r):=∫r∗rds(2/μ)(E−V(s))(r≥r∗) is finite for every r and satisfies ρ(t∗+τ(r))=r for all r≥r∗;
  4. for every time t1≥t∗ with ρ(t1)>R the motion is exactly free from t1 on: r(t)=r(t1)+(t−t1)r′(t1) for all t≥t1, and ρ is strictly increasing there; symmetrically, for t2≤t∗ with ρ(t2)>R, r(t)=r(t2)+(t−t2)r′(t2) for all t≤t2.

Facts & Assumptions

Given: The finite-range central-force data (μ,U,R), the energy E>0, the angular momentum j=Je3 with J>0, the simple outer turning point r∗ of V, and a trajectory with ∣r∣≥r∗ attaining r∗ at t∗.

[L1]

Polar reduction and energy: r=ρ n(θ), ρ=∣r∣, ρ2θ′=J/μ, μρ′′=−V′(ρ), and the constant energy is E=12μ(ρ′)2+V(ρ) (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation; here V=UJeff by Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential).

[L2]

Compact continuation and uniqueness for locally Lipschitz autonomous systems: solutions trapped in a compact subset of the open state domain have no finite endpoints, and two solutions with the same state at one time agree on every common interval (Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it, Picard-Lindelöf local existence and uniqueness for first-order systems, First-order systems, initial value problems, and solutions on intervals, Local Lipschitz continuity in the state variable, locally uniform in time and parameters).

[L4]

Behaviour of the quadrature: s↦E−V(s) is continuous and strictly positive on (r∗,∞), comparable to s−r∗ near r∗ (so 1/E−V(s) is integrable at r∗), and tends to E>0 as s→∞ (so τ(r)→∞, ∣ρ′∣≤2E/μ whenever ρ≥R, since there V(ρ)=J2/(2μρ2)≥0); continuous functions on intervals have primitives and the integral is additive over adjacent intervals (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G, Continuity of f:A→R at a point of A and on A: the ε-δ condition, its agreement with lim⁡x→cf(x)=f(c) at a limit point, and continuity at an isolated point, Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0).

Proof

technique · direct
1.1givenL1algebra

Square the energy identity of [L1]: 12μ(ρ′)2=E−V(ρ), which is clause 1's first display. At a time with ρ=r∗ the right side is E−V(r∗)=0, so ρ′=0; then μρ′′(t∗)=−V′(r∗)>0, so t∗ is a strict local minimum of ρ, and in particular ρ(t)>ρ(t∗) for t near but different from t∗.

2.1step 1.1L1algebra

On the open set where ρ>r∗, the right side of the energy identity is strictly positive, so ρ′≠0; hence the sign of ρ′ is constant on every connected component of {t:ρ(t)>r∗} (a union of open intervals on which the continuous nonvanishing function ρ′ has a sign).

2.2step 1.1L1L2L4algebra

Global extension: take the state z:=(ρ,ρ′,θ) in (0,∞)×R×R with right-hand side F(x):=(x2,−μ−1V′(x1),J/(μx12)). This field is locally Lipschitz: near any x1>0, continuity of V′′ bounds it by some M on a small closed interval; the fundamental theorem for continuous V′′ gives V′(a)−V′(b)=∫baV′′, and the Darboux definition with ∣V′′∣≤M gives ∣V′(a)−V′(b)∣≤M∣a−b∣. Also x1−2 is locally Lipschitz by factoring the difference of reciprocal powers. Along every solution, H(ρ,ρ′):=12μ(ρ′)2+V(ρ) is constant by the radial equation. Also ρ≥r∗ on every continuation: if it entered ρ<r∗, continuity from t∗ would give a boundary time with ρ=r∗ and values below r∗ arbitrarily nearby on one side; there H=E gives ρ′=0 and the equation gives ρ′′=−V′(r∗)/μ>0, contradicting the strict local-minimum property. Thus V(ρ)≤E and ρ′ 2=(2/μ)(E−V(ρ)). The function V is bounded below on [r∗,∞) because U is continuous and vanishes beyond R, so ρ′ is bounded; also 0<θ′=J/(μρ2)≤J/(μr∗2). On every bounded time interval, ρ and θ are therefore bounded, so near every finite right endpoint z lies in a compact subset of the state domain and [L2] continues it. For a finite left endpoint apply [L2] to z^(s):=z(−s), which solves z^′=−F(z^). Hence the maximal extension of z is defined on R. Set n(θ):=cos⁡θ e1+sin⁡θ e2 and n⊥(θ):=−sin⁡θ e1+cos⁡θ e2. The relation θ′=J/(μρ2) gives θ′′=−2ρ′θ′/ρ; differentiating gives r′′=(ρ′′−ρθ′2)n(θ)+(2ρ′θ′+ρθ′′)n⊥(θ), whose second coefficient vanishes. Using μρ′′=−V′(ρ) then yields μ(ρ′′−ρθ′2)=−U′(ρ), so the reconstructed curve satisfies the central-force equation. Any other continuation has the same angular momentum on the overlap: along it j′=μ(r′×r′+r×r′′)=0 by the central-force equation, so it remains in the same plane. By [L1] it has polar data there; shift its angle lift by an integer multiple of 2π to match θ at a common time. The resulting state solves the same locally Lipschitz system with the same initial state, so uniqueness [L2] gives the same ρ,θ and hence the same trajectory.

3.1step 1.1step 2.1algebra

Shape of the set {ρ>r∗}. This set is open in R, and on each of its connected components ρ′ has a constant sign (step 2.1). Every bounded component (t0,t1) satisfies ρ(t0)=ρ(t1)=r∗ (endpoints of a component of the open set {ρ>r∗} lie in its complement, where ρ≤r∗, i.e. ρ=r∗); at both endpoints ρ′′=−V′(r∗)/μ>0 by step 1.1, so ρ′ changes sign from negative to positive at each of them, whence just inside the component the sign is positive near t0 and negative near t1, contradicting constancy of the sign. Hence {ρ>r∗} has no bounded component. Since t∗ is an isolated zero of ρ−r∗ (a strict local minimum by step 1.1) and points just left and just right of t∗ belong to {ρ>r∗}, the two components adjacent to t∗ are (−∞,t∗) and (t∗,+∞): an open subset of R has at most two unbounded components, and each of these two is unbounded on its outer side (a bounded component is impossible, and it cannot end on the inner side at a point other than t∗, since its inner endpoint would be an endpoint of a component at which r∗ is attained, producing a second strict minimum and the same sign contradiction as above). Therefore ρ>r∗ exactly for t≠t∗, ρ′<0 on (−∞,t∗) and ρ′>0 on (t∗,∞).

4.1step 2.2step 3.1L3L4algebra

The quadrature: set q(s):=(2/μ)(E−V(s)) for s>r∗. The integral defining τ is finite and continuous at r∗ by [L4], and τ′(s)=1/q(s) for s>r∗. On the outgoing branch, the chain rule and ρ′=q(ρ) give ddtτ(ρ(t))=1 for t>t∗; taking the limit as t↓t∗ yields τ(ρ(t))=t−t∗. On the incoming branch, ρ′=−q(ρ) gives ddtτ(ρ(t))=−1, hence τ(ρ(t))=t∗−t. The function τ is strictly increasing on [r∗,∞) and tends to +∞ by [L4], so it maps this interval onto [0,∞). Inverting the outgoing identity gives ρ(t∗+τ(r))=r for every r≥r∗; the incoming identity gives the corresponding incoming branch. In particular the unique time with ρ=r∗ is t∗, and ρ(t)>r∗ for t≠t∗. These are clauses 2 and 3.

5.1step 2.2step 3.1L1algebra∎

Free tails: suppose t1≥t∗ with ρ(t1)>R. By step 3.1, ρ is strictly increasing on [t∗,∞), so ρ(t)≥ρ(t1)>R for all t≥t1. On that regime the potential is identically zero and the Cartesian central-force equation reduces to μr′′=−(U′(ρ)/ρ)r=0, since U′(ρ)=0 for ρ>R; hence r is affine on [t1,∞) and r(t)=r(t1)+(t−t1)r′(t1) there. The incoming statement is the same argument on (−∞,t2]. This is clause 4.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Existence of the asymptotic velocities and the impact-parameter identity J=μb∣v−∣

Statement

Let the hypotheses of Finite-range scattering: the simple outer turning point, the two global branches and the free tails hold, with finite-range potential U, energy E>0, angular momentum j=Je3, J>0, outer turning point r∗, and let r:R→R3 be the extension of the trajectory. Then:

  1. the limits v±:=lim⁡t→±∞r′(t) exist and satisfy ∣v+∣=∣v−∣=2E/μ;
  2. for every t2≤t∗ with ρ(t2)>R, the incoming free tail is the half-line {r(t2)+s v−:s≤0}, whose affine extension is the straight line {r(t2)+s v−:s∈R}; its perpendicular distance from the origin is b:=∣r(t2)×v−∣∣v−∣=Jμ∣v−∣,i.e.J=μb∣v−∣.

Facts & Assumptions

Given: The finite-range scattering orbit with its simple outer turning point r∗ and its two global branches, and times t1≥t∗, t2≤t∗ with ρ(t1),ρ(t2)>R.

[L1]

Free tails: for t1≥t∗ with ρ(t1)>R, r(t)=r(t1)+(t−t1)r′(t1) for all t≥t1; and symmetrically for t2≤t∗ with ρ(t2)>R, r(t)=r(t2)+(t−t2)r′(t2) for all t≤t2 (Finite-range scattering: the simple outer turning point, the two global branches and the free tails).

[L2]

On the free regime ρ>R the potential vanishes identically, so the constant energy is 12μ∥r′∥22=E there (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential, Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential).

[L4]

For nonzero u,v∈R3, the perpendicular distance from the origin to the line {u+sv:s∈R} equals ∣u×v∣/∣v∣: the closest point is u−⟨u,v⟩∣v∣2v, whose squared norm is ∣u∣2−⟨u,v⟩2∣v∣2=∣u×v∣2∣v∣2 by the Gram identity ∣u×v∣2=∣u∣2∣v∣2−⟨u,v⟩2 (The squared cross-product norm is the Gram determinant of two vectors, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The cross product in R3).

Proof

technique · direct
1.1givenL1

By [L1] the velocity is constant on the outgoing tail: r′(t)=r′(t1) for all t≥t1, so v+:=lim⁡t→∞r′(t)=r′(t1) exists; symmetrically v−:=r′(t2) is the constant velocity on the incoming tail t≤t2.

2.1step 1.1L2algebra

By [L2] applied at t1 and t2, ∥r′(t1)∥22=∥r′(t2)∥22=2E/μ; hence ∣v+∣=∣v−∣=2E/μ. This is clause 1.

2.2step 1.1givenL1algebra

The incoming line: by [L1] all points r(t) with t≤t2 lie on the set {r(t2)+s v−:s∈R}, a line through r(t2) with direction v−≠0 (it has norm 2E/μ>0).

3.1step 1.1step 2.2L3L4algebra∎

Its perpendicular distance from the origin is, by [L4] applied to u=r(t2), v=v−, ∣r(t2)×v−∣∣v−∣=∣r(t2)×r′(t2)∣∣v−∣=∣j∣/μ∣v−∣=Jμ∣v−∣, using that j=μ r(t2)×r′(t2) and ∣j∣=J from [L3]. Defining b to be this distance gives J=μb∣v−∣, which is clause 2.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

The scattering angle advance Δ and the deflection relation: the total polar advance is 2Δ and the velocity directions make the cosine cos⁡(π−2Δ)

Statement

Let the hypotheses of Finite-range scattering: the simple outer turning point, the two global branches and the free tails hold, let r:R→R3 be the extension of the trajectory, and let (ρ,θ) be polar data with r=ρ n(θ), n(φ)=cos⁡φ e1+sin⁡φ e2 (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation). Write t∗ for the apsis time, and for α∈R let Rα be the rotation of the plane span⁡(e1,e2), Rα(x):=cos⁡α x+sin⁡α (e3×x). Then:

  1. the half-angle Δ:=∫r∗∞J drr22μ (E−V(r)) is finite and positive, and θ(∞)−θ(t∗)=Δ=θ(t∗)−θ(−∞),θ(∞)−θ(−∞)=2Δ;
  2. lim⁡t→±∞n(θ(t)) exist and equal v+/∣v+∣ (as t→+∞) and −v−/∣v−∣ (as t→−∞);
  3. consequently v+∣v+∣=−R2Δ(v−∣v−∣),⟨v−,v+⟩∣v−∣ ∣v+∣=−cos⁡2Δ=cos⁡(π−2Δ); in particular, if 0<Δ<π/2 then π−2Δ∈(0,π) is the (undirected) angle between the asymptotic velocity directions, and the signed velocity rotation is −(π−2Δ) in the orientation defined by e3.

Facts & Assumptions

Given: The finite-range scattering orbit with its apsis time t∗, its polar data, its asymptotic velocities v± (Existence of the asymptotic velocities and the impact-parameter identity J=μb∣v−∣) and the frozen angular momentum j=Je3, J>0.

[L1]

Branch structure and quadrature: ρ(t∗)=r∗, ρ is strictly decreasing on (−∞,t∗] and strictly increasing on [t∗,∞), ρ′ has the corresponding signs, and (ρ′)2=(2/μ)(E−V(ρ)) (Finite-range scattering: the simple outer turning point, the two global branches and the free tails).

[L4]

The addition formulas give n(φ+α)=Rα(n(φ)) for all φ,α, and the cross product identity ⟨u,e3×u⟩=0 for every u: the mixed term of ⟨u,Rαu⟩ vanishes because ⟨u,e3×u⟩=det⁡[u,e3,u]=0 under the determinant form of the scalar triple product (The addition formulas for sine and cosine, Parity and the Pythagorean identity for sine and cosine, The cross product is bilinear, alternating, and orthogonal to both factors, The cross product in R3).

[L5]

Free tails: for t1≥t∗ with ρ(t1)>R, r(t)=r(t1)+(t−t1)v+ for t≥t1; symmetrically on the incoming tail with v− (Finite-range scattering: the simple outer turning point, the two global branches and the free tails, Existence of the asymptotic velocities and the impact-parameter identity J=μb∣v−∣).

Proof

technique · direct
1.1givenL1L2L3algebra

Outgoing half-angle: on [t∗,∞), θ is C1 and strictly increasing with θ′=J/(μρ2)>0, and ρ is a C1 strictly increasing bijection of [t∗,∞) onto [r∗,∞); substituting ρ for t first on compact subintervals of (t∗,∞), where ρ′>0, and then taking endpoint limits [L3] gives θ(∞)−θ(t∗)=∫t∗∞θ′(t) dt=∫r∗∞Jμs2dsρ′(ρ−1(s))=∫r∗∞J dss22μ (E−V(s))=Δ, using ρ′=(2/μ)(E−V) from [L1]; the integral converges at r∗ because E−V(s) is comparable to s−r∗ there and diverges nowhere else, and it is positive.

2.1step 1.1L1L2L3algebra

Incoming half-angle: the same computation on the decreasing branch (−∞,t∗], where ρ′=−(2/μ)(E−V(ρ)) and the substitution reverses the sign, gives θ(t∗)−θ(−∞)=Δ; adding the two halves gives θ(∞)−θ(−∞)=2Δ, which is clause 1.

3.1step 1.1step 2.1L5algebra

Direction limits: on the outgoing tail [L5], r(t)/∥r(t)∥2→v+/∣v+∣ as t→∞ (an affine curve with nonzero velocity direction tends to its velocity direction), and r(t)/ρ(t)=n(θ(t))→n(θ(∞)) because θ has a finite limit and n is continuous; hence n(θ(∞))=v+/∣v+∣. On the incoming tail, t−t2<0 tends to −∞, so r(t)/∥r(t)∥2→−v−/∣v−∣, while r(t)/ρ(t)→n(θ(−∞)); hence n(θ(−∞))=−v−/∣v−∣. This is clause 2.

4.1step 2.1step 3.1L4algebra

Rotation form: by [L4], v+/∣v+∣=n(θ(∞))=n(θ(−∞)+2Δ)=R2Δ(n(θ(−∞)))=−R2Δ(v−/∣v−∣), the first equality of clause 3.

5.1step 4.1L4algebra∎

Cosine: writing u:=v−/∣v−∣, the last identity gives ⟨v−,v+⟩/(∣v−∣∣v+∣)=⟨u,−R2Δu⟩=−cos⁡2Δ, because ⟨u,R2Δu⟩=cos⁡2Δ ∣u∣2+sin⁡2Δ ⟨u,e3×u⟩=cos⁡2Δ by [L4]. Since cos⁡(π−2Δ)=−cos⁡2Δ, this gives the displayed cosine; if 0<Δ<π/2 then π−2Δ∈(0,π) and the identity says that the angle between the two nonzero velocity vectors (the unique angle in [0,π] whose cosine is the displayed value, cosine being strictly decreasing there) is π−2Δ.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

One-dimensional substitution on a deflection branch: the outgoing angular measure has density b∣db/dθ∣, hence cross-section density b∣db/dθ∣/sin⁡θ

Statement

Assume ACω (The Axiom of Countable Choice (ACω)). Let θ:(a1,a2)→(θ1,θ2)⊆(0,π) be a deflection branch of Scattering data: the incident flux measure, the solid-angle element, deflection branches and the cross-section density, with inverse b(θ), and set D=(a1,a2)×(0,2π), H=(θ1,θ2)×(0,2π) and T(b,φ)=(θ(b),φ). Restrict the incident measure ν to D. Its pushforward, defined by T∗ν(B)=ν(T−1(B)) for Borel B⊆H, satisfies T∗ν(B)=∫Bb(θ)∣b′(θ)∣ dλ2(θ,φ). For every nonnegative Borel g:H→[0,+∞], and for every real Borel g with either side's integral of ∣g∣ finite, consequently ∫Dg(θ(b),φ)b dλ2(b,φ)=∫Hg(θ,φ)b(θ)∣b′(θ)∣ dλ2(θ,φ). The nonnegative equality allows +∞; in the real case both sides are finite and absolute integrability on one side is equivalent to that on the other. In particular it applies to continuous periodic test functions that are nonnegative or absolutely integrable, and to compactly supported continuous test functions. No equality of undefined signed improper integrals is asserted. Nonnegative double-integral notation denotes the corresponding Lebesgue integrals, with iterated integrals justified by Tonelli.

Thus the outgoing angular measure has density b(θ)∣b′(θ)∣ against dθ dφ and density dσdΩ(θ)=b(θ)∣b′(θ)∣sin⁡θ=b(θ)sin⁡θ ∣θ′(b(θ))∣ against the solid-angle measure. For finitely many branches on disjoint incident intervals the measures add, including on overlapping image intervals; in particular contributions with pairwise disjoint images add.

Facts & Assumptions

Given: Countable choice, the branch, its inverse, D,H,T,ν and the test function as stated.

[L2]

Compact substitution with absolute derivative holds on every closed interval lying strictly inside the branch; continuous functions there are Riemann integrable and their Riemann and Lebesgue integrals agree under countable choice (In one dimension the compact-Jordan formula is substitution over the unoriented image interval with the absolute derivative, A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, A bounded Riemann integrable function on a closed bounded interval is Lebesgue measurable and has the same integral).

[L3]

On Borel subsets of the plane λ2=λ1×λ1, and Tonelli applies to nonnegative Borel integrands. Measures agreeing on a generating pi-system with an increasing common finite-measure exhaustion agree everywhere (On Borel subsets of R^{m+n}, the product lambda_m times lambda_n agrees with lambda_{m+n}, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Measures agreeing on a generating pi-system are equal under an increasing finite-measure exhaustion from that pi-system).

[L5]

sin⁡θ>0 on (0,π): the solid-angle definition proves this from sin⁡0=sin⁡π=0 and strict increase/decrease on the two half intervals; its continuous density is Borel measurable (Scattering data: the incident flux measure, the solid-angle element, deflection branches and the cross-section density).

Proof

technique · direct
1.1givenL1L4

By [L1], b is C1 and h(θ)=b(θ)∣b′(θ)∣ is continuous and positive. The map T is a homeomorphism D→H, so T∗ν is a Borel measure: preimages are Borel and preserve disjoint countable unions. The nonnegative density h defines another Borel measure η(B)=∫Bh dλ2 by [L4].

2.1step 1.1L2L3algebra

Let c<d lie in (θ1,θ2) and 0<p<q<2π. Apply compact substitution [L2] to the inverse map b on [c,d], with integrand x↦x: writing u=min⁡(b(c),b(d)), v=max⁡(b(c),b(d)), it gives ∫cdh(θ) dθ=∫uvx dx. Endpoints have Lebesgue measure zero, so Tonelli [L3] yields η((c,d)×(p,q))=(q−p)∫uvx dx=ν(T−1((c,d)×(p,q))). These quantities are finite.

3.1step 1.1step 2.1L3

Interior open rectangles, together with the empty set, form a pi-system generating the Borel sets of H: rational rectangles form a countable base, and every open set is the union of those base rectangles it contains. Choose explicit nested interior intervals exhausting each open coordinate interval; their product rectangles exhaust H and have the common finite measure of step 2.1. Thus measure uniqueness [L3] gives T∗ν=η on every Borel set of H, proving the asserted pushforward identity even for unbounded open branches.

4.1step 3.1L3L4algebra

The test-function identity first holds for indicator functions by step 3.1 and the definition of pushforward, then for nonnegative simple functions by finite addition, and finally for every nonnegative Borel g by increasing simple approximation and monotone convergence [L4] on both sides. Apply this to ∣g∣ to obtain equivalence of absolute integrability; applying it to the positive and negative parts gives the real identity whenever that absolute integral is finite. Tonelli [L3] supplies the nonnegative iterated form.

5.1step 1.1step 3.1step 4.1L1L4L5algebra∎

Since sin⁡θ>0 by [L5], rewriting h dλ2=(h/sin⁡θ) dΩ gives the solid-angle density; [L1] gives its reciprocal-derivative expression. Pushforward commutes with finite addition of measures by its definition, so branches on disjoint incident intervals contribute the sum of these measures, and the sum of the densities wherever their images overlap. This proves all assertions.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Repulsive Coulomb scattering: the unique outer turning point, global non-colliding motion, the asymptotic impact distance and Δ=π/2−arctan⁡(κ/(2Eb))

Statement

Let μ>0, let κ>0 and U(r)=κ/r on (0,∞) (the repulsive Coulomb potential, of class C2), let E>0 and J>0, and put V(r):=κr+J22μr2,b:=J2μE>0,α:=κEb>0. Then:

  1. V is strictly decreasing on (0,∞), with V(r)→∞ as r↓0 and V(r)→0 as r→∞; hence there is exactly one r∗>0 with V(r∗)=E, it is simple (V′(r∗)<0), and V(r)<E for all r>r∗.

  2. Let r:I→R3 be a central-force trajectory for (μ,U) with constant angular momentum j=Je3, J>0, energy E, with ∣r(t)∣≥r∗ for all t∈I and with ∣r(t∗)∣=r∗ for some t∗∈I. Then r extends uniquely to all of R; along the extension ρ=∣r∣≥r∗>0 (no collision), ρ is strictly decreasing on (−∞,t∗] and strictly increasing on [t∗,∞), and the asymptotic velocities v±=lim⁡t→±∞r′(t) exist with ∣v+∣=∣v−∣=2E/μ. Moreover the asymptotic perpendicular distance of each branch to its limiting direction exists and equals b=J2μE=Jμ∣v±∣,so thatJ=μb∣v±∣. The distance in clause 2 means the limit of ∣r(t)×v±∣/∣v±∣ as t→±∞, the distance to the line through the origin in the limiting velocity direction; it does not assert a free tail or a finite longitudinal intercept.

  3. For the polar angle, with u∗:=1/r∗, Δ:=∫r∗∞J drr22μ (E−V(r))=∫0u∗b du1−κEu−b2u2=π2−arctan⁡κ2Eb.

  4. Write n±=lim⁡t→±∞r(t)/∣r(t)∣. Then n+=v+/∣v+∣, n−=−v−/∣v−∣, and in an oriented polar plane the outgoing velocity direction is the incoming direction rotated by 2Δ−π∈(−π,0). In particular ⟨v−,v+⟩∣v−∣∣v+∣=cos⁡(π−2Δ),θs=π−2Δ=2arctan⁡κ2Eb. For fixed μ,κ,E this gives an actual family for every b>0, with angular momentum J(b)=2μE b: choose its unique turning radius and initial data r(0)=r∗e1, r′(0)=J(b)e2/(μr∗). Local existence and the extension in clause 2 give a member of the family for every b, unique up to apsis time and oriented plane placement.

Facts & Assumptions

Given: The repulsive Coulomb data μ,κ,E,J, the effective potential V, the outer turning point r∗, and a scattering trajectory with apsis time t∗.

[L1]

Polar reduction and energy: r=ρ n(θ), ρ=∣r∣, ρ2θ′=J/μ, μρ′′=−V′(ρ), E=12μ(ρ′)2+V(ρ) (Polar reduction of a non-radial central-force trajectory: ρ2θ′=J/μ, the effective-potential equation and energy conservation); here V=UJeff (Central-force trajectories: the central-force equation, angular momentum, energy and the effective potential).

[L3]

The Cartesian field G(x,v):=(v,κx/(μ∣x∣3)) is defined on the open set D:={(x,v):∣x∣>0} and is C1 there. Its acceleration derivative is Dx(κx/(μ∣x∣3))[h]=κμ(∣x∣−3h−3∣x∣−5x⟨x,h⟩). Thus DG is bounded on a sufficiently small convex ball about each point of D; along a segment in such a ball, the chain rule and the vector-valued fundamental theorem of calculus give G(z)−G(w)=∫01DG(w+s(z−w))(z−w) ds, and the norm inequality for vector integrals gives a finite Lipschitz bound. Hence G is locally Lipschitz, so Picard–Lindelöf uniqueness and compact continuation apply on D (Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it, Picard-Lindelöf local existence and uniqueness for first-order systems, First-order systems, initial value problems, and solutions on intervals, Local Lipschitz continuity in the state variable, locally uniform in time and parameters, If f:[a,b]→Rm is differentiable with integrable f′ then ∫abf′=f(b)−f(a); and a bounded derivative makes f Lipschitz, For a≤b and f:[a,b]→Rm integrable when a<b, ∥∫abf∥2≤∫ab∥f∥2; for a<b, ∥f∥2 is integrable, Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, Ck Euclidean maps and diffeomorphisms).

[L5]

Trigonometry (The derivatives of sine and cosine are cosine and minus sine gives sin⁡x/x→1 from the derivative at zero): arctan⁡:R→(−π/2,π/2) is the inverse of tan⁡; arcsin⁡, arccos⁡ are the principal inverses; the addition and subtraction formulas hold; sin⁡2+cos⁡2=1; cos⁡ is strictly decreasing on [0,π] and sin⁡ is strictly increasing on [−π/2,π/2] (The principal inverse tangent arctan⁡:R→(−π/2,π/2), Principal inverse sine and inverse cosine, The addition formulas for sine and cosine, The subtraction formulas for sine and cosine, Parity and the Pythagorean identity for sine and cosine, Signs, monotonicity intervals, and ranges of sine and cosine, Tangent, cotangent, secant, and cosecant on their exact natural domains).

Proof

technique · direct
1.1givenL2algebra

The derivative V′(r)=−κ/r2−J2/(μr3) is negative on (0,∞), so V is strictly decreasing by [L2]; κ/r→∞ and J2/(2μr2)→∞ as r↓0, while both terms tend to 0 as r→∞. Hence V is a continuous strictly decreasing bijection (0,∞)→(0,∞); since E>0 lies in the range, r∗=V−1(E) is the unique solution of V=E, it is simple because V′<0 everywhere, and V(r)<E for r>r∗. This is clause 1.

2.1step 1.1L1L3algebra

Global continuation: energy and angular momentum are conserved by [L1], and V(ρ)≤E forces ρ≥r∗ since V is strictly decreasing. Also ∣r′∣≤c:=2E/μ because U>0. At any finite endpoint, ∣r∣ is bounded above by its value at a fixed time plus c times the elapsed time, and r stays at least r∗ from the origin. Thus the Cartesian state (r,r′) stays in a compact subset of {(x,v):∣x∣>0} near that endpoint. The Cartesian field (v,κx/(μ∣x∣3)) is smooth there; compact continuation [L3] extends the solution past every finite right endpoint. For a finite left endpoint apply the same argument to the time-reversed state (r(−s),−r′(−s)), which solves the same autonomous Cartesian system; uniqueness [L3] then makes the global extension unique. The maximal extension therefore has domain R.

3.1step 1.1step 2.1L1L2L4algebra

Branch shape: ρ′′=−V′(ρ)/μ>0 for all times by step 1.1. Hence ρ′ is strictly increasing, vanishes at t∗ by the energy identity, is negative before t∗ and positive after it. Consequently the apsis is unique and ρ>r∗ elsewhere. Fixing any t1>t∗ gives ρ′(t)≥ρ′(t1)>0 for t≥t1, so ρ→∞ forward; the same argument backward gives ρ→∞ there. On either open branch, compact substitution [L4] and then limits give ∣t−t∗∣=∫r∗ρ(t)ds/(2/μ)(E−V(s)). Its lower endpoint is integrable since E−V(s) is comparable to s−r∗ by V′(r∗)<0; its upper endpoint diverges since the integrand tends to 1/c>0. This proves the branch and quadrature claims without a finite-range hypothesis.

4.1step 3.1L1L4L6algebra

Asymptotic velocities and their speed: since V(ρ(t))→0 along each branch (as ρ(t)→∞), ∣ρ′∣→2E/μ. Put c:=2E/μ and c0:=c/2>0. The radial-speed limit implies ∣ρ′(t)∣≥c0 for sufficiently large ∣t∣; integrating on each monotone branch gives ρ(t)≥c0∣t∣−C0 for a finite constant C0. In particular ρ(t)≥(c0/2)∣t∣ after increasing the tail cutoff. The acceleration of the Cartesian curve is μr′′=−(U′(ρ)/ρ)r, so ∥r′′∥2=∣U′(ρ)∣/μ=κ/(μρ2)≤C/(1+∣t∣)2 for large ∣t∣, an integrable bound along each half-line; the Cauchy criterion for the vector integral (the integral of r′′ over [s,t] is bounded by ∫st∥r′′∥2) therefore gives existence of v±=lim⁡t→±∞r′(t). Passing to the limit in the energy identity E=12μ∥r′∥22+U(ρ(t)) with U(ρ(t))→0 gives ∣v±∣=2E/μ.

4.2step 3.1L1L4algebra

Half-angle quadrature: on compact subintervals of either open branch use θ′=J/(μρ2) and ∣ρ′∣=(2/μ)(E−V) in [L4]. Taking limits gives the same positive half-angle on both sides of the apsis, Δ=∫r∗∞J dss22μ(E−V(s))=∫0u∗b du1−(κ/E)u−b2u2. Each substitution is first on an interior compact interval. The integral converges at r∗ by the square-root bound of step 3.1 and at infinity by the bound C/s2, so it gives finite limits θ± with θ+−θ∗=θ∗−θ−=Δ.

5.1step 4.2L4L5algebra

Evaluation with scaled roots: put w=bu, A=1+α2/4 and w∗=(α2+4−α)/2=bu∗, where α=κ/(Eb). Then 1−αw−w2=A2−(w+α/2)2. Substitute w+α/2=Asin⁡ψ on interior compact intervals, where ψ∈(ψ0,π/2) and ψ0=arcsin⁡(α/(2A))∈(0,π/2). The derivative is Acos⁡ψ>0 and the positive denominator is Acos⁡ψ, so the integral equals π/2−ψ0 after taking limits. Since sin⁡ψ0=α/α2+4 and cos⁡ψ0=2/α2+4, tan⁡ψ0=α/2, giving ψ0=arctan⁡(α/2) by [L5]. Hence Δ=π/2−arctan⁡(κ/(2Eb))∈(0,π/2). This proves clause 3.

6.1step 4.1step 4.2step 5.1L1L5algebra

Direction limits: differentiating r=ρn(θ) gives r′=ρ′n(θ)+(J/(μρ)) n⊥(θ), where n⊥(θ)=−sin⁡θ e1+cos⁡θ e2. By step 4.1, ρ′→c forward and ρ′→−c backward, while the tangential term tends to zero. Step 4.2 therefore gives v+=c n(θ+) and v−=−c n(θ−). The addition formulas [L5] and θ+−θ−=2Δ show the velocity rotation 2Δ−π and the displayed cosine. There is no hidden extra winding: the radial coefficient ρ′ changes sign only once, the tangential coefficient is always positive, so the velocity direction relative to n(θ) varies continuously from π to 0 through the upper half-plane; its net change is 2Δ−π.

7.1step 3.1step 4.1step 4.2step 6.1L1L4L5L6algebra

Impact limit: put q(s)=(2/μ)(E−V(s)) for s>r∗. On the outgoing branch the remaining angle is Ψ(t)=∫ρ(t)∞J/(μs2q(s)) ds, so ρ(t)Ψ(t)=Jμ∫1∞duu2q(ρ(t)u). As q(s)→c>0, for any ε>0 one has ∣1/q(s)−1/c∣<ε for every s≥S for some S. Once ρ(t)≥S, the difference between the integral above and 1/c is at most ε∫1∞u−2 du=ε. Thus ρΨ→J/(μc) and Ψ→0; the elementary limit sin⁡Ψ/Ψ→1 gives ρsin⁡Ψ→J/(μc)=b. By step 6.1 this is exactly ∣r(t)×v+∣/∣v+∣. The incoming computation uses q=∣ρ′∣ and the remaining angle θ(t)−θ− and gives ∣r(t)×v−∣/∣v−∣→b as well.

8.1step 1.1step 2.1step 3.1step 5.1step 6.1step 7.1L1L3algebra∎

For every b>0 set J(b)=2μE b and take the turning radius of clause 1. The initial data stated in clause 4 have angular momentum J(b)e3 and energy V(r∗)=E, radial derivative zero and positive radial second derivative. Local existence [L3] supplies a trajectory and conservation [L1] forces ρ≥r∗; the preceding continuation, angle and impact-limit steps apply to it. Uniqueness [L3] determines it from its apsis state. Thus this actual family realises the deflection formula for all b>0, and all clauses are proved.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

The Rutherford deflection map, its inverse, the cross-section density (κ/4E)2csc⁡4(θ/2) and the divergence of the total cross-section

Statement

Assume ACω (The Axiom of Countable Choice (ACω)) for the cross-section measures. Fix μ>0, κ>0 and E>0. For every b>0 put J(b)=2μE b and use the Coulomb trajectory family of Repulsive Coulomb scattering: the unique outer turning point, global non-colliding motion, the asymptotic impact distance and Δ=π/2−arctan⁡(κ/(2Eb)) clause 4. Its impact distance is b and its deflection map θdef:(0,∞)→(0,π),θdef(b):=2arctan⁡κ2Eb. Then:

  1. θdef is C1 and strictly decreasing, with θdef(b)→π as b↓0 and θdef(b)→0 as b→∞; it is a C1 bijection (0,∞)→(0,π) with inverse b(θ)=κ2Ecot⁡θ2(0<θ<π).
  2. For every interval (θ1,θ2)⊆(0,π) the pair (θ,b) is a deflection branch with cross-section density dσdΩ(θ)=b(θ)∣dbdθ(θ)∣1sin⁡θ=(κ4E)2csc⁡4θ2, finite and positive on (θ1,θ2); and the pushforward of the incident flux measure under the branch is described by One-dimensional substitution on a deflection branch: the outgoing angular measure has density b∣db/dθ∣, hence cross-section density b∣db/dθ∣/sin⁡θ.
  3. For every θ0∈(0,π) the truncated total cross-section ∫0θ0 ⁣ ⁣∫02π(κ4E)2csc⁡4θ2 sin⁡θ dφ dθ=+∞, while for 0<θ1<θ2<π the integral over (θ1,θ2)×(0,2π) is finite; in this sense the total cross-section of the ideal Coulomb model diverges at forward angles and every compact angular acceptance away from 0 has finite cross-section.

Facts & Assumptions

Given: Countable choice, the fixed data μ,κ,E, the family J(b)=2μE b, the deflection map θdef and its asserted inverse b(θ).

[L1]

arctan⁡:R→(−π/2,π/2) is the inverse of tan⁡, arctan⁡ is C1 and strictly increasing, and tan⁡(arctan⁡y)=y for every real y; on 0<θ/2<π/2 one has arctan⁡(tan⁡(θ/2))=θ/2 (The principal inverse tangent arctan⁡:R→(−π/2,π/2), Tangent is a continuous strictly increasing bijection from (−π/2,π/2) onto R, Principal arctangent: derivative, integral, power series, and the Gregory–Leibniz series).

[L3]

The half-angle advance for the Coulomb branch is Δ=π/2−arctan⁡(κ/(2Eb)) and the deflection magnitude is π−2Δ=2arctan⁡(κ/(2Eb)) (Repulsive Coulomb scattering: the unique outer turning point, global non-colliding motion, the asymptotic impact distance and Δ=π/2−arctan⁡(κ/(2Eb)) clauses 3-4).

[L6]

A continuous nonnegative function with a positive value has positive integral, the integral is monotone, and for 0<ε<θ0 the computed integral ∫εθ0θ−3 dθ=12(ε−2−θ0−2) diverges as ε↓0; sin⁡x≤x for x≥0 (mean value theorem applied to sin⁡ on [0,x], using cos⁡≤1) and cos⁡(θ/2)≥cos⁡(θ0/2)>0 for 0<θ≤θ0≤π/2 (The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a), If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a), A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G, The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a), Parity and the Pythagorean identity for sine and cosine, Signs, monotonicity intervals, and ranges of sine and cosine, For a<c<b: f is integrable on [a,b] if and only if it is integrable on [a,c] and on [c,b], and then ∫abf=∫acf+∫cbf; with the oriented form for arbitrary a,b,c, Sums, scalar multiples, products, absolute values, maxima, minima and quotients with nonvanishing denominator of continuous functions are continuous, as are constants, the identity and every polynomial function, Continuity of f:A→R at a point of A and on A: the ε-δ condition, its agreement with lim⁡x→cf(x)=f(c) at a limit point, and continuity at an isolated point).

Proof

technique · direct
1.1givenL1L2L4algebra

Differentiation of the branch: by [L4] and [L2], θdef′(b)=2⋅11+(κ2Eb)2⋅(−κ2Eb2)<0(b>0), so θdef is C1 with strictly negative derivative, hence strictly decreasing; as b↓0, κ/(2Eb)→∞ and arctan⁡ tends to π/2, so θdef→π; as b→∞, κ/(2Eb)→0 and θdef→0. The image is therefore (0,π) and the inverse exists, is C1 by [L4], and satisfies the claimed formula: for 0<θ<π, arctan⁡(tan⁡(θ/2))=θ/2 by [L1] (since 0<θ/2<π/2) and tan⁡(θ/2)=κ/(2Eb(θ)) with b(θ)=κ2Ecot⁡(θ/2), because tan⁡(θ/2)cot⁡(θ/2)=1 by [L2]; substituting shows θdef(b(θ))=θ.

2.1step 1.1L3algebra

The half-angle and deflection relation of [L3] applied for each b>0 to the corresponding trajectory with J=J(b) give exactly the deflection value π−2Δ=2arctan⁡(κ/(2Eb)); thus the abstract branch formula is the computed scattering deflection of the Coulomb problem, and clause 1 is complete.

2.2step 1.1L2L5algebra

Cross-section density: by [L2], b′(θ)=−κ4Ecsc⁡2(θ/2), so ∣b′(θ)∣=κ4Ecsc⁡2(θ/2); combining with [L5], dσdΩ(θ)=b(θ)∣b′(θ)∣sin⁡θ=κ28E2⋅cot⁡(θ/2)csc⁡2(θ/2)sin⁡θ. Using sin⁡θ=2sin⁡(θ/2)cos⁡(θ/2) from [L2], cot⁡θ2csc⁡2θ2/sin⁡θ=cos⁡(θ/2)sin⁡3(θ/2)⋅12sin⁡(θ/2)cos⁡(θ/2)=12csc⁡4θ2, whence dσ/dΩ=κ216E2csc⁡4(θ/2)=(κ/(4E))2csc⁡4(θ/2), finite and positive on (θ1,θ2); this is clause 2.

3.1step 2.2L6algebra

Divergence: for 0<θ≤θ0 with θ0≤π/2, the integrand of clause 3 equals (κ/(4E))2sin⁡θ csc⁡4(θ/2)=(κ/(4E))2⋅2cos⁡(θ/2)/sin⁡3(θ/2) by the same identities; since sin⁡(x)≤x for x≥0 and cos⁡(θ/2)≥cos⁡(θ0/2)>0 by [L6], the integrand is at least Cθ−3 with C=16(κ/(4E))2cos⁡(θ0/2)>0. Hence for 0<ε<θ0, by monotonicity of the integral and the φ-integration of a constant, ∫εθ0 ⁣ ⁣∫02π(κ/(4E))2csc⁡4(θ/2)sin⁡θ dφ dθ≥2πC∫εθ0θ−3 dθ=πC(ε−2−θ0−2)⟶+∞ as ε↓0, by the fundamental theorem and [L6]; hence the truncated integral over (0,θ0) diverges. If instead π/2<θ0<π, cut the interval at π/2 and use the divergence on (0,π/2] together with additivity.

4.1step 3.1L6algebra∎

Finiteness on compact angular intervals: the integrand (κ/(4E))2csc⁡4(θ/2)sin⁡θ is continuous on the compact rectangle [θ1,θ2]×[0,2π] because sin⁡(θ/2)>0 there, so its integral is finite by [L6]; this completes clause 3.

5 · Examples, experiments and counterexamples

None yet.

Sources