Write , and . The polar amplitude equations are
with their reflected counterparts. On the equal-amplitude standing wave branch , equilibrium point requires
Squaring and adding eliminates the phase:
For the candidate squared amplitudes are
For each positive root, set and . The quadratic identity ensures these define a phase. Conversely these phase equations imply the quadratic, so a nonzero phase-locked branch exists exactly when the discriminant is nonnegative and at least one root is positive.
Since , the discriminant condition is the printed for . It is sufficient under the usual above-onset supercritical assumptions , because then the plus root is positive. Without these assumptions the condition is not sufficient: satisfies it but gives . This is the positivity condition for a phase-locked standing wave, not a removable algebraic detail. The zero solution always exists and is distinct from the nonzero branch. For , the undivided equations must be used; a nonzero equilibrium point in this fixed rotating frame requires and . If , nonzero equilibrium points require and a compatible phase, with no cubic amplitude selection.
Strong enough forcing can overcome the nonlinear frequency shift, lock the sum phase to the external signal and produce a standing wave at half the forcing frequency. The difference phase remains free through spatial translation. The amplitude inequality is an existence condition, not a stability theorem: attraction requires the linearized forced amplitude equations, and nonlinear coefficients have not been restricted enough here to promise stability of every root.
A persistent saddle loop must return with the same energy. The Melnikov energy-balance method therefore requires, at first order,
For the normalized orbit of part (c), put . Since , and ,
For the general orbit has and . The energy mismatch is therefore proportional to , whose derivative with respect to is nonzero. This simple zero selects the saddle-loop curve:
The remainder can be sharpened. Reversal of and interchanges and without changing the saddle-loop condition. Smoothness of the homoclinic splitting thus makes the selected an even function of , giving
The leading coefficient alone is enough to separate it from the Hopf curve . The integral hint in the PDF equals for nonnegative even integers ; for odd its integrand is odd and the integral is zero. Only even powers enter the calculation above.
Near the double-zero point, and with the stipulated absence of further bifurcations, the distinct open regions are:
On the homoclinic curve the cycle becomes a saddle loop; its positive saddle trace makes the loop repelling on the cycle side. At the Hopf curve the unstable cycle shrinks into the negative equilibrium point. The node–focus discriminant curves merely change the spiral appearance and are not additional bifurcations. The impossibility of a periodic orbit for can also be seen without a picture: at a local minimum of periodic one would have and , but at its maximum the same strictly positive expression would be required to be nonpositive.
Figure 1.
Local saddle-node, subcritical Hopf and homoclinic curves with the four nearby phase portraits
. Details.
The ring of integers is , since , and its discriminant is . Rationalizing gives
Both are algebraic integers, so is a unit. Set . Since ,
Also . It is the unique ideal of ideal norm two: any quotient of size two is , and the image of must satisfy , hence be one. The kernel of a ring homomorphism onto is uniquely , which is .
The given Minkowski constant means every ideal class has an integral representative of ideal norm at most . Such a ideal norm is one or two. The norm-one ideal is the whole ring, and the sole norm-two ideal is principal. Thus the class number is one.
If , the principal ideal has ideal norm two and hence equals . Its generator differs from by a unit. Conversely multiplication by a norm-one unit preserves the field norm two. With the assumed fundamental unit, all units are ; in fact a norm-minus-one unit is impossible since has no solution. Therefore the complete answer is
For and the positive sign, the product is , giving ; checks the result.
The factorization form of Kummer-Dedekind theorem is as follows. If with integral, minimal polynomial , and does not divide the index , factor into distinct monic irreducibles over . Then
The index condition is essential. All integral bases used below make the index one, so it applies to every prime.
For a squarefree quadratic , take , and field discriminant . At an odd prime, a repeated factor occurs exactly when , in which case and . At , reduction is if is even and if is odd, so is totally ramified in both cases. Thus the totally ramified primes are exactly those dividing .
For , take , with polynomial and discriminant . At an odd prime it has a double root exactly when , giving a single linear factor squared and total ramification. At , its derivative is , so no repeated factor occurs. Hence again
For the prime cyclotomic field, use . Modulo this is , so , proving total ramification. For , the discriminant has no prime factor , so the reduction is squarefree and is unramified. There is an actual exponent error in the printed discriminant: the correct value is , not . To check it, at a nontrivial root ,
and multiplying these derivatives and using gives the exponent . The prime support is still just , so the requested ramification conclusion is unaffected.
For , the integral basis gives and . Modulo , ; modulo , . Dedekind factorization gives and . For other primes the polynomial is squarefree because its discriminant is not divisible by the prime. Therefore the equivalence also holds in this cubic field. This is a conclusion for these specific fields, not the false general assertion that every ramified prime in every higher-degree number field is totally ramified.
For the given Möbius transformation, the fixed point equation is , hence
For the general Möbius transformation, its finite fixed points satisfy
With and nonzero discriminant, the quadratic formula gives the two distinct fixed points of a Möbius transformation:
The fixed-point equation for also gives . Therefore
Dividing yields the conjugation of a Möbius transformation with two distinct fixed points:
Neither factor in this quotient vanishes, since . The formula extends over the Riemann sphere, including the transformation's pole. For the special transformation and the ordering , .
The Torelli theorem concerns a connected smooth projective curve over and its canonical principal polarization of a complex torus. Its conclusion is
The polarization of a complex torus is part of the hypothesis. We will reconstruct the curve intrinsically from the Jacobian variety and the theta divisor, using the Gauss map of a theta divisor, and prove the reconstruction in both the nonhyperelliptic and hyperelliptic curve cases. The needed biduality and ordinary-tangency facts about projective dual varieties are the general projective-duality results permitted here.
First explain why the polarized data determine the theta divisor up to translation. On , take integral one-form classes dual to the symplectic cycle basis. Pullback by the degree-one Abelian sum map identifies them with the corresponding classes on ; the Riemann bilinear relations give
On , the sum map pulls each one-form back to the sum of its pullbacks from the factors. It has generic degree onto : the Abel theorem for divisors and the existence of classes with give a unique effective unordered representative, and ordering gives the factorial. Integrating a wedge of basis one-forms over the product, only terms assigning a matching pair to each factor survive. Dividing by gives integral one when precisely one matching pair is omitted, and zero for all other basis wedges. Hence the Poincaré dual class of is
This is exactly the integral form of the canonical principal polarization of a complex torus constructed by periods.
Its associated line bundle has one-dimensional section space. One can check this without any theta multiplicity assumption: in normalized period coordinates a representative has sections satisfying
Expand the first periodicity in a Fourier series . The second periodicity forces . Positive definiteness of makes this series normally convergent, so the section space has dimension exactly one. Other line bundles giving the same polarization of a complex torus differ by a degree-zero line bundle. The map is an isomorphism from to its dual abelian variety for a principal polarization of a complex torus: on lattices it is the unimodular alternating form. Thus all these line bundles are translates, and so are their unique zero divisors. In particular the given polarized data recover the reduced divisor on an algebraic curve up to translation. Translation does not change tangent hyperplanes after identifying tangent spaces with .
Put and . At a smooth class , the Riemann singularity theorem says . The derivative of the Abelian sum map identifies the translated tangent hyperplane with . Thus the intrinsic Gauss map of a theta divisor is
In particular
We will see that this map is dominant and generically finite. Its function-field extension therefore determines an intrinsic finite cover: normalize in , obtaining . Normalization in a finite extension is finite here. Let be its reduced codimension-one branch locus. Using this finite cover avoids counting omitted points of the rational Gauss map of a theta divisor as ramification. Everything in this construction is determined by the polarized Jacobian variety.
Suppose first that and is nonhyperelliptic. Its canonical map embeds it as a nondegenerate curve . For clarity, this follows directly from the Riemann-Roch theorem: on a positive-genus curve, so has no base point. If a length-two effective divisor had , a nonconstant meromorphic function with poles bounded by would give a map of degree at most two to . Degree one forces genus zero, and degree two is the hyperelliptic curve case. Hence , and the Riemann-Roch theorem gives . The canonical sections therefore separate every pair of points and every tangent direction, proving the embedding.
Let be the open set of hyperplanes transverse to . A hyperplane in cuts out distinct points. These points span the hyperplane: if a second independent canonical section vanished at all of them, its zero divisor would have to equal the same degree- divisor, and the ratio of the two sections would be a constant. Consequently some of the intersection points are independent. For the divisor formed by such a subset, and .
We need to know that every -subset of a general hyperplane section has this property. Here is the monodromy argument, with its essential steps. The incidence space choosing one intersection point is a projective bundle over before restricting to , so it is irreducible. The incidence space choosing two distinct ordered intersection points is likewise a projective bundle over minus its diagonal, with fibre ; it too is irreducible. Restricting to transverse sections preserves irreducibility and hence connectedness of these covers. Thus the monodromy action of a covering space on the intersection points is two-transitive. A general tangent hyperplane has just one contact of order two, by characteristic-zero projective biduality theorem; locally its intersection equation is together with the other simple points. A small loop in exchanges exactly those two points. Two-transitivity conjugates this transposition to every transposition, so the monodromy is the full symmetric group .
It follows that the finite étale cover of choosing a -subset is connected and irreducible. The condition of dependent evaluations is closed on that cover. It is proper because we have exhibited an independent subset; its image under the finite map is therefore a proper closed subset of . Outside that image every subset is independent. By the Abel theorem for divisors, distinct such subsets have distinct Abelian sums: means there is only one effective representative of its class. Thus the subset cover is birational to with its Gauss map of a theta divisor, and
More precisely it is exactly the restriction of over , since this étale cover is normal and has the same function field.
There can consequently be no branch component over . Its complement is the irreducible hypersurface of tangent hyperplanes, the projective dual variety of . Around a general point of , the transposition of two roots acts nontrivially on -subsets: choose a subset containing exactly one of those two roots. Hence the finite cover is ramified there. This proves the exact identity
The quoted biduality theorem now recovers as an embedded curve, and the canonical map identifies with . This proves reconstruction in the nonhyperelliptic case, including special curves rather than only a general curve in moduli.
Now suppose that is a hyperelliptic curve of genus . Let be its degree-two map. The Riemann-Hurwitz formula gives branch points . Choose an affine coordinate with no branch at infinity and write
The forms are holomorphic and independent; their orders at the two points at infinity show there are exactly of them. Thus the canonical map factors as followed by the degree- rational normal curve embedding
Write .
A general hyperplane section of consists of distinct points avoiding the branch points. The corresponding canonical divisor on consists of the pairs . Choosing one lift from each pair gives an effective divisor with : the evaluation matrix of the canonical forms has Vandermonde rank . If instead a degree- sub-divisor contains an entire pair, it has at least two sections, because that pair is a fibre of and moves in a pencil. Every choice not taking exactly one point from each pair contains an entire pair. Therefore the generic smooth theta points over the hyperplane are exactly the independent choices of one lift in each pair, and
The cover of these choices is again the restriction of to the open set of transverse sections of avoiding the : the identification on a dense open follows from the Abel theorem for divisors, and the choice cover is finite étale and normal. Its components would correspond to components of the function field of the irreducible theta divisor, so it is connected.
There are precisely two types of codimension-one branching. If a simple root crosses a branch point , its two lifts exchange, giving nontrivial monodromy on the choices. This yields the hyperplane
If two roots collide away from the , a local transverse parameter gives roots . Trivialize the two sheets of near . Choices selecting different sheets for the two roots are exchanged when the roots exchange; hence this also causes branching. Its locus is , the hypersurface of tangent hyperplanes to . Away from these two loci all roots and all their lifts vary locally without branching, so there are no other branch components. Consequently
For , is irreducible and nonlinear. Indeed, hyperplane sections of the degree- rational normal curve are binary forms of degree , and the repeated-root hypersurface is their irreducible discriminant of degree . Thus intrinsically separates into its single nonlinear component and its linear components. Biduality recovers , while dualizing recovers the point on .
These data recover the double cover uniquely. Identify with and move infinity away from the recovered branch set. Its function field is
Any other quadratic extension with the same simple branch points has a defining rational function whose ratio to this product has even valuation everywhere. On such a rational function is a constant times a square, and every nonzero complex constant is a square. The quadratic extensions are therefore isomorphic; their smooth projective models are isomorphic. This proves the hyperelliptic reconstruction. The branch locus also distinguishes the two cases: the nonhyperelliptic locus has one irreducible branch hypersurface, whereas the hyperelliptic locus has the additional linear components.
Finally handle the small genera. In genus two, the theta divisor is , the image of itself. The Abelian sum map on is injective because a principal divisor with would give a degree-one map to . Its derivative is nonzero since the canonical bundle has no base point, so this compact injective immersion identifies with . In genus one the same argument identifies with its Jacobian variety after choosing an origin. In genus zero the Riemann-Roch theorem supplies a function with one simple pole, giving . These cases complete the proof: the canonically polarized Jacobian determines the curve in every genus.
There is an error in the printed request: the cubic is not irreducible at , and three elements cannot be a field basis in that case. The integral-basis assertion is valid under the necessary hypothesis .
Under that hypothesis, is an algebraic integer and has finite index . Change of basis in the trace pairing gives the discriminant-index identity
The field discriminant is an integer. If the left side is square-free, cannot divide it unless . Thus is an integral basis when the cubic is irreducible and its discriminant is square-free.
For any root, gives
Since , and ,
On the other hand, the product of is : the product of the roots is one, and each of the three unordered root pairs contributes a negative squared difference. Therefore, without needing irreducibility for this polynomial calculation,
The three-root product is the field norm only when the cubic is the minimal polynomial.
At , the polynomial has no rational root (the only candidates are ), so it is irreducible, and its discriminant is . Hence is an integral basis.
At ,
If , then with integral basis . If , then and is an integral basis, since this quadratic order has square-free discriminant five; also . In neither case is the three-element list a basis. Its discriminant five is instead that of the degree-three algebra ; a basis there does not repair the printed assertion about the field .
Use the usual positive-definite convention for binary quadratic forms with negative discriminant. Reduction by integral changes of variables of determinant one produces a form with , and with on the boundary. For completeness, replacing by makes ; if then , interchange the variables, reversing one sign to keep determinant one. The new positive leading coefficient is smaller. Repeating terminates and gives a reduced form.
For discriminant , and the reduction inequalities give . Thus , is odd and , and . The boundary convention chooses . Hence there is exactly one positive-definite equivalence class, represented by .
For an odd prime represented by this form, the identity shows that is a square modulo . Here , because otherwise and would divide the represented number. Conversely, if is a square modulo , choose odd with . Then is an integral positive-definite form of discriminant . By the preceding reduction it is equivalent to the given form, so its value at is represented by that form. Quadratic reciprocity gives , whose value is one precisely for . The exceptional primes are represented: and . Therefore
Without the positive-definite convention, the first printed assertion is false: also has discriminant but cannot be equivalent to a positive-definite form, because an invertible change of variables preserves its sign.
Since , the ring of integers is with discriminant . Write and . Reduction of modulo two gives a double root, and modulo five gives the roots . Thus , , and . Directly , since . The element belongs to and and has field norm ten, so equality of ideal norms gives
The Minkowski bound is , so every ideal class contains an integral ideal of ideal norm at most five. Three is inert and gives no prime ideal of ideal norm three; ideals of ideal norm four are , principal. Thus the class group is generated by , and the displayed relations make every class either trivial or . The class is nontrivial: a generator would have field norm , but would give or , neither a quadratic residue. Hence the class number is two.
The unit has field norm . To show it is fundamental, suppose is a positive real embedding of a unit . Its conjugate is , so is a positive integer less than two. Thus , and forces . Only is greater than one, contrary to the strict inequality. Therefore no smaller unit exceeds one. Multiplying any positive unit by a suitable power of puts it in , proving that all units are .
An element of field norm has principal ideal or by prime ideal factorization. The first is generated by ; the second by , of field norm . Two generators of the same principal ideal differ by a unit. Therefore all, and only, the integer solutions are
The two displayed signs are independent; multiplication by a norm-minus-one unit exchanges the two field norm signs.
For a nonidentity element of represented by a determinant-one matrix, the classification is elliptic, parabolic or hyperbolic according as the squared trace is less than, equal to or greater than four. This follows either from the eigenvalues or from the discriminant of its fixed-point equation on the boundary of the upper half-plane.
The first matrix has determinant one and trace two, and is not scalar. It is therefore parabolic. Its fixed-point equation is
Thus its unique fixed point is on the real boundary, of multiplicity two; infinity is not fixed and there is no fixed point in the upper half-plane. The second matrix has determinant one and trace , so it is hyperbolic. Its two boundary fixed points are , consistently with this classification.
For and , a nonzero infected equilibrium requires . Using and population conservation gives
Eliminate before assessing stability, since perturbations must preserve the total population. The resulting Jacobian matrix at this endemic equilibrium is
Its trace is negative and its determinant positive, so both eigenvalues have negative real part: the endemic equilibrium is locally asymptotically stable on the fixed-population plane. The full three-variable system has a neutral direction corresponding to changing the conserved population.
Writing , as the trace is and the determinant is . Hence the discriminant . As , the trace is and the determinant is again , giving a positive discriminant . Both quantities are in both limits, but the eigenvalues are complex for sufficiently small and real for sufficiently large .
A rational point of order can be moved to , giving an integral model
Its two-isogeny descent partner and the two-isogeny formula are
The maps extend over the identity and , have these two points as their respective kernels, and direct substitution and the doubling formula give and .
Define the two-torsion square-class homomorphisms by
and similarly for , using . To see the homomorphism property, intersect a line with . Its three -coordinates have product , since the resulting monic cubic has constant term . The third point is the negative of the sum of the first two, and negation preserves , so the square classes multiply as required. If the line passes through , the other two coordinates have product , giving exactly the prescribed exceptional value. Vertical lines and repeated intersections give the identity and tangent cases.
We also verify the kernel needed for descent:
For a point on , the first coordinate of its image under is , a square. Conversely, if has square -coordinate, solving gives
Its discriminant is , so it has rational nonzero roots . Taking gives a point on ; choose the signs so its image has second coordinate . The exceptional point lies in the image exactly when is a square, as seen from the points with , . This proves the first kernel formula. Applying the same argument to and then rescaling the second partner by proves the second.
There are only finitely many candidate square classes. For a rational point write and , with . These denominator forms follow because a negative prime valuation of makes the leading cubic term the unique lowest-valuation term, forcing that valuation to be even. Then
If a prime divides but not , the second factor is a unit at that prime, so its valuation in is even. Thus every class in has a signed squarefree representative . This is the prime-support bound in two-isogeny descent; the analogous candidates for divide .
The class is realized by a nonexceptional rational point precisely when its quartic covering in a two-isogeny descent
has a rational point with . Indeed substitution , gives the original equation, and conversely the square-class condition gives these parameters. Clearing denominators permits integral with . The zero-coordinate cases account for the already known identity and two-torsion classes.
For clarity, the two-isogeny rank formula contains a factor , even when has full rational two-torsion. Put , which is or , and let be the rank. From the two kernel formulas and ,
Indeed the kernel of the induced map is . The intersection is , of size , so . The Mordell-Weil theorem gives . Therefore
This gives a practical rank procedure. Enumerate the finitely many signed squarefree divisors of , and test their quartics over the reals and local fields at and primes dividing the relevant coefficients or elliptic-curve discriminant. These necessary local conditions give finite Selmer group of an elliptic curve upper bounds. Search the remaining coverings for rational points, or search the curves directly, to realize square classes and obtain lower bounds. One can also exhibit independent points and check independence by the canonical height pairing. When the upper and lower bounds agree, the rank is determined. Local solubility alone does not guarantee a rational point, so a persistent gap can reflect a nontrivial Tate–Shafarevich group rather than a missing congruence test.
Now take , for an odd prime . The curve has four rational two-torsion points . Their images under are . The prime-support bound allows no other classes, so
Its partner is . Any nonidentity rational point with has , because . Also the exceptional point has square class . Thus its image is contained in the four positive candidate classes . The rank bound for a prime congruent-number elliptic curve follows:
For an exact rank-one example, take . The point on realizes the class . The five-adic obstructions for the congruent-number isogeny covers exclude the other two possible nontrivial classes. For , the quartic is
Modulo , if the right side is , a nonsquare. Hence and ; primitivity gives . Dividing by then gives , again impossible. For , the quartic is . At least one of is nonzero modulo , and their fourth powers are or , so their sum is or modulo . The right side has five-adic valuation , impossible for a square. Thus exactly, and
The dual isogeny maps to on , giving an explicit rational point associated with this rank-one computation.
Write . Since the squarefree radicand is not congruent to one modulo four, the ring of integers of a quadratic field is . The discriminant of the basis is , giving
Write . The ideals are , and , with residue fields of orders 2, 3 and 3 respectively. Thus their ideal norms are 2, 3 and 3.
The generators of are . They are multiples of 2, and the integer combination of 4 and 10 giving 2 proves . Similarly is generated by . The difference of the middle generators is 6; together with 9 this generates 3. Therefore the product is .
The ideals and are coprime, so their product equals their intersection. The element belongs to both, hence . Both ideals have norm 6, so equality holds. Likewise belongs to and , and its norm is . Thus
The discriminant is 40. The Minkowski bound for ideal classes in a real quadratic field is , so every ideal class has an integral representative of norm 1, 2 or 3. The only possible prime ideals with those norms are the ones just listed. Their classes satisfy and . Thus there are at most two classes. The ideal is not principal, since a generator would have norm , already ruled out. Therefore the class number is exactly two.
Away from the velocity-profile corners, and an unstable normal mode satisfies . On the upper half of the jet write
The exterior solution is selected by decay at infinity. From the given velocity components, the perturbation pressure is
Continuity of normal velocity and pressure therefore gives continuity of and pressure matching at a piecewise-linear shear interface:
At the slope of jumps from to zero. Using yields
The lower interface then follows by parity.
For the even sinuous mode of a planar jet, the derivative can jump at the center: being an even function does not require the one-sided derivative to vanish. The slope jump at is , whereas . Thus
Substitute into the upper-interface relation. With , division by and multiplication by give
Collecting powers of the phase velocity gives the required dispersion relation:
These equations can also be treated as a homogeneous two-by-two system in , avoiding division by any potentially zero amplitude or phase-velocity factor.
For the odd varicose mode of a planar jet, , so ; the central pressure matching at a piecewise-linear shear interface is automatic. The upper-interface equation becomes , and therefore
This phase velocity is real for every , so the discrete varicose mode of a planar jet is neutral and there is no exponentially unstable odd mode.
To determine the unstable even modes, factor the reduced discriminant of the quadratic:
The second factor vanishes at zero and has derivative for , so it is positive there. The first factor starts at , is strictly increasing, and tends to infinity. It therefore has one positive zero . The even phase velocities are
It follows that
For , one member of the complex conjugate pair grows at rate and the other decays. At the speeds merge and there is no exponential growth; for both speeds are real. This is inviscid instability of a triangular jet. The unstable sinuous mode of a planar jet also satisfies Howard's semicircle theorem with and : its speed lies in the upper half-disk centered at with radius .
Figure 1.
Triangular jet profile and unstable sinuous phase velocities inside the Howard semicircle
. Details.
inversion in a circle with center and radius is , with and infinity interchanged. It fixes the circle pointwise and reverses complex orientation. Each inversion is fractional-linear in , so composing two eliminates the conjugation and gives with nonzero determinant: a Möbius transformation.
Translate and rotate two distinct centers to , where , and let the radii be . The composition of inversion in the first circle followed by the second is
Its fixed points solve . The discriminant is
Disjoint circle boundaries require either or , so both factors have the same nonzero sign and . There are two distinct fixed points on the line of centers. If the circles are concentric, their radii differ and the composition is , with exactly the two fixed points . Thus the composition has two distinct fixed points, also for nested disjoint circles. Reversing the order inverts the Möbius transformation and leaves its fixed points unchanged.
Irreducibility modulo two implies irreducibility over by Gauss lemma for polynomials, so the Galois group embeds transitively in . The discriminant is , hence primes two and three are unramified. The factor-degree theorem says that square-free reduction with factor degrees provides an element of the Galois group with cycle lengths . The reductions already found therefore give a 4-cycle and a 3-cycle.
The group order divides 24 and is divisible by four and three, so is 12 or 24. An order-12 subgroup of has index two and must be : in its quotient every transposition has the same nontrivial image, since transpositions are conjugate and generate , making the quotient map the sign. But a 4-cycle is odd. Therefore the order is 24 and
The discriminant can also be obtained from and the resultant, giving .
Write and . Since their common square-free discriminant is congruent to zero or one modulo four, it must be one modulo four; thus are odd. It is not a rational square because a square-free integer square is only zero or one. Both polynomials are irreducible over , and their roots generate the same quadratic extension . Their integral roots are and , respectively. Each generated order is
the usual ring of integers for square-free . This identifies both quotients with the same ring. Explicitly, if denotes the image of modulo , send the image of modulo to . Equal discriminants give , and the opposite translation gives its inverse. The two quotient rings are isomorphic.
The cubic has no rational root: its only possible integer roots are , and neither works. It is therefore irreducible. Let be its image in ; it is an algebraic integer, so . The power-basis discriminant is
The index-discriminant formula from part (i) says that the square of divides 59. It must be one. Hence , with integral basis and field discriminant .
Write the inward radial velocity as , with . Steady spherical mass conservation and the radial Euler equations for an inviscid fluid give
For a barotropic fluid, define the sound speed and the barotropic enthalpy function by . Choosing , integration gives the Bernoulli equation
The continuity equation gives . Eliminating mass density from the momentum equation produces the barotropic spherical accretion equation
The sonic point is where the inward Mach number equals one. A smooth solution cannot cross the vanishing coefficient on the left with a nonzero right-hand side, so a regular sonic radius must obey
Conversely a continuous flow that is subsonic at infinity and supersonic inward must cross unit Mach number; smooth passage selects a simultaneous zero of the two sides, hence a regular sonic point. The associated critical accretion branch, rather than an arbitrary subsonic solution with smaller flux, is selected by these outer and inner requirements.
For the polytropic equation of state, set to distinguish its exponent from any separately specified perturbation exponent. Assume and first take the usual positive index . Then
For the attractive mixed potential, assume . The sonic condition and Bernoulli equation give
Hence polytropic accretion in a mixed inverse-power potential reduces to
For all nonzero coefficients in the left side are nonnegative for , with at least one strictly positive: there is no finite positive root. For the constant term is negative and the leading coefficient positive, giving exactly one positive root. Define and . Then
This is the finite-radius transonic range for the two positive attractive terms.
The critical point has real distinct crossing slopes. To check this rather than just solve the sonic equalities, put , and . Differentiating the barotropic spherical accretion equation at its simultaneous zero gives
Its discriminant is . The sonic Bernoulli equation says , so this discriminant is strictly positive for . The accretion crossing is the slope for which , namely the smaller root. The two slopes have opposite signs of this Mach-number derivative. This is the regular transonic saddle, not a singular vertical tangent.
The corresponding mass accretion rate follows without an additional arbitrary mass density normalization:
One may equivalently replace by . With a nonzero term, inner free fall has and continuity gives . Thus increases without bound precisely in the positive-index range , consistent with the required inner supersonic branch.
The assumptions on the parameters matter. If and , the finite sonic radius instead requires , the usual Bondi accretion range; is the zero-radius limiting critical case. If and , the condition remains . With arbitrary signs of , one must use the quadratic and retain only positive roots with and real admissible crossing slopes; attraction was not explicitly given as a sign condition in the source.
If positive polytropic index was not intended as an implicit physical convention, there is a further mathematical range. For , positive compressibility also permits , giving . Both potential terms attractive then give a positive sonic root as well:
The same mass accretion rate formula applies; its crossing-slope discriminant is positive. The range has and no physical acoustic sonic point, has constant pressure and zero sound speed, and does not define the printed equation of state. Thus the usual answer assumes and positive attractive coefficients; these conventions should not be confused with additional printed restrictions.