Write , , and retain the frozen axial field . The radial magnetostatic equilibrium equation now includes both the gradient of the toroidal magnetic pressure and the inward magnetic tension:
For , cylindrical gravity gives
Multiplication by the integrating factor therefore produces
The integration constant must vanish for regularity on the axis. Direct integration gives
Its Taylor expansion at the axis is
so the regular field behaves as . At large radius,
Nonnegativity at infinity requires
This condition is also sufficient: at the limiting value, the braces divided by reduce to , and decreasing only increases them. Thus the toroidal magnetic field is real and regular at every radius exactly in the stated range.
Let
Flux conservation gives . Substitution into the extensional-force equation yields the third-order equation
Since
one integration, using and its derivatives tending to zero on the flat-film side, gives
Put . After division by , this becomes
The integrating factor gives
A second integration and therefore give
On the bubble side, matching to its cylindrical shape gives , so . Taking in the integrated equation yields . Hence
With the convention in the question, the constant velocity-anisotropy parameter is . Since , the Spherical Jeans equation becomes
For and , its general integrating factor solution is
The homogeneous term represents a boundary pressure. For an extended scale-free system the physical boundary condition normally removes it, leaving
Positivity requires . A complete physical model must also have a nonnegative galactic distribution function and sensible inner and outer boundary behaviour; for example, strong radial anisotropy is restricted by density-slope--anisotropy inequalities.
Observationally, the tracer density can be estimated from star counts only after correcting distances, extinction, survey selection, and incompleteness. Spectroscopy supplies mainly line-of-sight velocities; proper motions add transverse information but become less precise for distant halo stars. The equation shows the mass--anisotropy--density degeneracy directly: the same measured can result from a larger , a steeper tracer slope , or a different . Even globally constant power laws therefore do not determine the galactic mass profile unless some of these quantities are independently constrained.
If , , or changes near a break radius, the solution at one radius also depends on the outer boundary integral. A break in observed dispersion may be attributed to a mass-profile feature, a tracer-density break, or a change in orbital anisotropy. Separate tracer populations, full three-dimensional velocities, higher velocity moments, and measurements over a wide radial range help break this degeneracy.
More flexible alternatives model a nonnegative solution of the Collisionless Boltzmann equation itself. An action-based galactic distribution function gives an analytic or parametrized ; a Schwarzschild orbit-superposition model assigns nonnegative weights to an orbit library; and a made-to-measure stellar-dynamical model adjusts particle weights to reproduce observations. These methods retain more phase-space information than Jeans moments, although their flexibility introduces model choices and regularization.
For a Keplerian orbit, write , so . Substituting the stated torque laws into
gives
With
the integrating factor is . The boundary condition then gives
Thus
and
Let the source's hydrogen-ionizing photon production rate be
where is the hydrogen ionization energy. The idealized Strömgren sphere has an almost fully ionized interior and a thin ionization front. In photoionization equilibrium, every ionization is balanced by a Case B recombination, so spherical symmetry gives
For pure hydrogen of constant number density , the interior has , and the Strömgren radius is therefore
Now let the effective number of dust grains per hydrogen nucleus be , so the dust absorption coefficient is . If is the ionizing-photon rate crossing the sphere of radius , recombinations and dust absorption give the linear ordinary differential equation
If denotes a dust mass fraction instead, the grain mass and gas mean particle mass are simply absorbed into the effective product . Define the dust optical depth
Multiplication by the integrating factor and integration to the dusty front gives its governing equation
The exponential function weights recombinations at large optical depth by the extra source photons that dust must remove before those photons reach that radius.
For constant , put and . The elementary integral
reduces the equation to
Equivalently, with and ,
Its dust-free limit is , while absorption makes for nonzero dust abundance.
The odd copy number decreases by two until it reaches one, so . The exact first-moment equation is
Because copy numbers are nonnegative integers, , and the last mixed moment is nonnegative. Hence
The integrating factor , together with , gives the comparison bound
Subtracting the exact first-moment equation from the equation cancels reaction 5:
The assumed inequality is precisely
Choose a positive smaller than this gap divided by . The limit superior bound implies that, for all sufficiently large ,
Thus grows at least linearly. Since ,
so the required species index is .
The critical-point equations are and , so the points are and . The Hessian is
Thus is a strict minimum, a strict maximum, and are saddles. The local contour patterns are respectively ellipses, inverted ellipses, and crossing hyperbolas.
The gradient flow has and . Hence
Putting gives . Multiplication by the integrating factor yields
and therefore the general trajectory is
Let increase downward. The slender viscous thread is locally in uniaxial extension. The transverse stress equals the ambient pressure, so the Trouton ratio gives the excess axial stress . Mass conservation and axial force balance therefore give
In steady flow . Dividing the momentum equation by gives
With and , choose
The dimensionless equation becomes
Treating as a function of and using an integrating factor yields
Because , the area sketches are the reciprocals of the functions found below. The dimensional vertical deviatoric stress is
Substitute part (ii) into . Multiplication by the integrating factor gives
Since , integration and rearrangement yield
Complexifying the directional transport equation yields spectral solutions analytic inside and outside the unit circle. An integrating factor built from known attenuation removes its multiplication term. The limiting Green kernels encode weighted line measurements and transverse Hilbert transforms. Their jump gives a Riemann-Hilbert problem; its solution recovers the source from an asymptotic coefficient. A general primary treatment of this spectral construction is www.stat.uchicago.edu/~guillaumebal/PAPERS/AttRadTh.pdf .