The probability that an isotropically oriented circular orbit produces an exoplanet transit, assuming . Uniformity of gives the displayed result. Orbital eccentricity, a restriction to full rather than grazing crossings, or a non-isotropic orientation distribution changes the expression.
Gravitational field of a finite uniform line segment Created 2026-10-18 Updated 2026-10-20
Spread mass uniformly along on an axis, with linear mass density . In its perpendicular bisector plane, axial field components cancel, and the radial component is
This is an exact line mass model and a thin-rod approximation for a physical cylinder. A centred circular orbit in this symmetry plane has . Thus for , and for . For a finite-width rod, the first limit is useful only while its transverse size remains much smaller than .
A circular infinite cylinder of radius and uniform mass density has linear mass density . A coaxial Gaussian surface of length has inward radial gravitational field and no end-cap flux. Gauss's law for gravity gives outside, hence the displayed result. Inside, enclosed linear mass density is , so . An exterior circular orbit has speed independent of radius. The exterior formula presupposes circular axial symmetry; area alone would not specify the near field of a noncircular rod.
An orbit around an isolated source constrains its standard gravitational parameter rather than independently measuring the source's mass and the gravitational constant. In the test particle approximation, a circular orbit of radius and orbital period satisfies
The left equation equates the required centripetal acceleration to the gravitational acceleration. More generally, a Kepler orbit depends on the same . In the exact relative two-body problem, replace by .
For any positive dimensionless number , replacing by and by leaves and every predicted test particle trajectory unchanged. This is the gravitational mass degeneracy: more observations of the same orbit can refine without separating its factors. A laboratory source whose mass is calibrated independently supplies the missing information.
Consequently an astronomical mass expressed in kilograms as inherits measurement uncertainty from , even if is known much more precisely. To first order,
When the relative standard uncertainty of is negligible, that of is approximately . Conversely, the ratio of two source masses inferred using the same is , so this common factor cancels. That distinction explains how excellent orbit determination can coexist with less precise absolute astronomical masses.
For a binary star viewed far from its centre of mass, let be the separation and the reduced mass, the total mass, and . The Taylor expansion of the point-mass Newtonian gravitational potential has zero dipole because . Its second moment is , so
This holds for an instantaneous separation in any orbit, not only a circular orbit. The mass quadrupole moment convention gives .
Facing a companion and facing the source of light are different relations. The Moon keeps its near side of the Moon approximately toward Earth, while the direction to the distant Sun is nearly fixed during the idealized orbit shown below. A fixed surface patch therefore turns toward, across and away from sunlight. The far side of the Moon has day and night too; its name does not mean permanent darkness. [ref].
Figure 1.
A face fixed toward Earth still alternates between sunlight and darkness
. Four snapshots of one circular-orbit cartoon viewed from above the orbital plane. The orange dot marks the same Earth-facing surface patch in every snapshot; the Sun is far to the left and all rays travel rightward. White and grey show directly lit and unlit hemispheres. The diagram neglects Earth's motion around the Sun, lunar libration, eclipses, terrain and reflected light. Sizes and distances are arbitrary; this is a geometry cartoon, not an ephemeris or photograph. Details.
If the Sun direction has fixed angle and the marked surface patch has outward normal vector at angle , its direct illumination is proportional to . With synchronous rotation, , so the patch goes through day and night even though it always faces Earth. In reality the Earth-Moon system advances around the Sun: the lunar spin and orbital cycle relative to distant stars is about Earth days, while the lunar phase cycle has a synodic period of about Earth days. These numbers describe the actual system, not the fixed-sunlight cartoon. [ref].
For a synchronously rotating exoplanet whose orbital companion is itself the illuminating star, the direction to the star turns with the orbit. On an ideal circular orbit, with perpendicular spin axis and no additional oscillations, this direction is fixed relative to the surface, giving persistent day and night hemispheres. Day-night heat redistribution then helps determine the nightside temperature; permanent lack of direct starlight does not imply zero temperature. Orbital eccentricity and departures from this ideal spin geometry modify that simple picture.
Kepler's first law 2026-10-13
The bound relative Kepler orbit is an ellipse with the attracting centre at one focus of a conic. In polar coordinates, its separation is the displayed function of true anomaly , the angle measured from periapsis. For and , it becomes , giving at and at . Setting orbital eccentricity to zero gives a circular orbit. This specifies the path; Kepler's second law specifies how it is traversed.
Kepler's third law 2026-10-13
For an isolated two-body problem in Newtonian gravity, the relative elliptic orbit's orbital period obeys
Here is the relative semi-major axis, and is the sum of the two masses. The separation varies during an elliptic orbit, so replacing by works only for a circular orbit, where . Likewise, either body's individual semi-major axis about the centre of mass must first be converted to the relative one. When is negligible compared with , the denominator can be approximated by .
For fixed gravitational parameter, the orbital period depends on the semi-major axis and is independent of orbital eccentricity. Our , example gives both for and for a circular orbit with . The eccentric body's rapidly changing speed compensates for its different path. A compact way to see the cancellation uses Kepler's second law: divide the ellipse's area by its areal velocity . With and specific angular momentum , the result is .
Compare two separate Kepler orbits with the same gravitational parameter, and :
With , their orbital periods are and ; measured in units of , these are and . Quadrupling the size therefore takes eight times as long to complete a circuit. This power law makes comparisons possible without first choosing physical units.
Figure 1.
Quadrupling the relative semi-major axis multiplies the orbital period by eight at fixed total mass
. Details.
The same relation can be inverted to weigh a system: . An observed relative semi-major axis and orbital period determine the total mass, while additional information is needed to separate the two masses. This is useful for binary stars, as discussed in NASA, Orbits and Kepler's Laws: measuring masses. Geometry and a clock together become a weighing instrument.
Last-orbit frequency of a compact binary Created 2026-09-24 Updated 2026-10-14
As a heuristic approximation, inserting into the circular orbit form of Kepler's third law, with here denoting the total mass, gives
The dominant quadrupole gravitational-wave frequency is twice the orbital frequency, .
For a spherical body of radius , mass and mean density , a low circular orbit has period , by equating the gravitational acceleration from Newton's law of universal gravitation to the centripetal acceleration. Substitution proves the displayed relation. Thus similar low-orbit periods indicate similar mean densities, even for bodies with very different masses and radii. At altitude , the right side is multiplied by ; the simple density inference requires and neglects rotation, nonsphericity and other gravitating bodies.
Lunar phase 2026-10-14
A lunar phase is the appearance of the sunlit portion of the Moon as seen from Earth. Outside a lunar eclipse, the distant Sun illuminates approximately one hemisphere of the Moon, but an observer sees a changing overlap of that hemisphere with the hemisphere facing Earth. The lunar phase therefore changes even though synchronous rotation keeps approximately the same near side of the Moon facing Earth. [ref].
The cycle runs through new Moon, waxing crescent, first quarter, waxing gibbous, full Moon, waning gibbous, last quarter and waning crescent. “Waxing” means that the lunar illuminated fraction is increasing; “waning” means it is decreasing. A crescent has less than half its apparent disk sunlit, and a gibbous lunar phase has more than half. First and last quarter each show half a disk: “quarter” refers to progress around the cycle, rather than a quarter of the disk being lit. [ref].
An instructive circular orbit model puts Earth at the origin, the Sun far away along , and the Moon at . The sunward position , , gives new Moon; the opposite position , , gives full Moon. Counterclockwise motion gives first quarter at and last quarter at . The Sun's direction is fixed during this snapshot of the cycle. Every small globe in the orbital view is illuminated on its half; the adjacent apparent disks are separate views from Earth, rather than globes viewed from above.
Figure 1.
Eight lunar phases in a circular orbit and the corresponding apparent disks from Earth
. Details.
The disk views choose north up, with the upward direction perpendicular to the orbital plane. In this convention waxing lunar phases are bright on the right and waning lunar phases on the left. Local horizon orientation and viewing from different locations on Earth can rotate the apparent disk; its lit edge always points toward the projected direction of the Sun. The lunar illuminated fraction is unaffected by this rotation.
The useful angle is the phase angle , with the Sun, the Moon and the observer at Earth's center: its vertex is at the Moon. Full Moon has , quarter has , and new Moon has in the ideal model. By contrast, the elongation has its vertex at Earth. For a distant Sun, its rays are parallel, so
Here is the lunar illuminated fraction: a fraction of projected disk area, not a fraction of the Moon's total surface area or its received radiant flux relative to full Moon. The phase angle and unsigned elongation each lie between and , so either alone gives the illuminated amount without distinguishing waxing from waning. The orbital direction resolves that ambiguity.
Figure 2.
Projected illuminated fraction versus the Sun–Moon–observer phase angle at the Moon
. Details.
The lunar illuminated fraction has an elegant geometric derivation. Model the Moon as a sphere of unit radius and view it along the direction. Take the Sun's unit direction to be . A visible surface point above apparent disk coordinates has . Direct illumination requires the outward normal's dot product with the sunward direction to be positive. The day-night terminator therefore obeys . At a fixed , set ; the visible day-night terminator is , and the illuminated strip extends from to the right edge . Integrating its width gives
The integral is the area of half a unit circle. Thus a three-dimensional visibility problem reduces to a circle's area and a cosine. The derivation assumes parallel sunlight, a distant observer and a sphere, and ignores a lunar eclipse. Real observed brightness additionally depends on surface scattering of light and observing geometry, so does not promise half the full-Moon radiant flux.
For a worked example, elongation gives phase angle and : a crescent with of its apparent disk sunlit. At elongation , the phase angle is and : a gibbous lunar phase. At , . In each case the same occurs once while waxing and once while waning in the ideal cycle.
Ordinary lunar phases are not caused by Earth's shadow. A lunar eclipse occurs only when the Moon actually enters that shadow near full Moon. The orbital inclination of the Moon is about relative to Earth's orbital plane around the Sun; most full Moons pass above or below the shadow. This is why full Moon does not imply a lunar eclipse every month. [ref].
The phase cycle is a synodic month, about days, whereas the Moon's orbital period relative to distant stars is a sidereal month, about days. During one sidereal month, Earth advances around the Sun, so the Moon must move farther to repeat the sunward alignment. For coplanar circular orbits traversed in the same direction, subtract their mean motions: , with . Hence
Using and gives , consistent with the mean synodic month of days. After one sidereal month, the sunward direction has advanced about ; catching up takes about more days. These are mean-period estimates: orbital eccentricity and changing geometry make individual synodic months vary. [ref]. [ref].
The lunar phase is also a practical clock for observing. Full Moon is opposite the Sun and rises roughly at sunset; first quarter is about a right angle from the Sun and is convenient in the evening, while last quarter is convenient before dawn. These are rough guides, not exact rise-time predictions for every place and season. [ref]. Historically, repeated lunar phases provided the month in lunar calendars, while half-lit disks gave Aristarchus of Samos a way to infer astronomical distances. [ref].
Newtonian potential of a point mass Created 2026-09-29 Updated 2026-10-14
The Newtonian potential of a point mass at the origin is
It is the free-space Green function of the Laplacian with the normalization of the Poisson equation for Newtonian gravity. Differentiating gives inward gravitational acceleration . With , the normalized Newtonian gravitational potential is ; at it equals .
For a radially outward launch from , conservation of mechanical energy in the fixed-source test particle approximation gives the minimum escape velocity:
At this threshold the speed tends to zero only at infinite separation. The model neglects air resistance, external tidal forces, and general relativity. [ref]
A circular orbit at instead requires tangential speed . Choosing units with gives , , and . A radial launch at speed reaches and returns; the same speed directed tangentially gives a circular orbit. [ref]
Figure 1.
Energy balance and outward radial speeds below, at, and above escape
. Left: at the escape velocity, positive kinetic energy exactly balances negative Newtonian gravitational potential energy. Right: radial launches at normalized speeds , , and respectively turn at , approach zero speed at infinity, and retain a positive limiting speed. The plotted trajectories use a fixed source and the test particle approximation. Details.
Orbital speed 2026-10-18
Orbital speed is the magnitude of the velocity of an orbiting body in a specified reference frame. For a Newtonian circular orbit around a point source of mass , equating to gives . In a curved spacetime, a speed inferred from coordinate radius and coordinate time generally differs from a speed measured by a local observer. For a circular orbit in the Schwarzschild metric, a static local observer measures ; uses the clock normalization at infinity instead.
Paczyński-Wiita circular orbit Created 2026-10-05 Updated 2026-10-07
Let . In the Paczyński-Wiita potential, circular orbit balance and the radial epicyclic frequency formula give
Thus the innermost stable circular orbit lies at . Its specific orbital energy is , giving a zero-torque radiative efficiency of black-hole accretion of . The model reproduces the Schwarzschild inner orbit exactly but approximates its relativistic energy and frequency. Abramowicz's derivation explains the potential's relation to the Schwarzschild effective potential.
Expand a unit-mass particle's Lagrangian about an axisymmetric circular orbit of radius and frequency , using and . Circular-orbit balance removes the linear radial term; total-time-derivative invariance of a Lagrangian removes . To second order,
Here is the orbital shear parameter and is the vertical epicyclic frequency. The Euler-Lagrange equations retain the local Coriolis acceleration and tidal gravity.
For a Newtonian circular orbit, , so the specific orbital energy is . Half of the released potential energy remains in circular orbital motion. The binding-energy estimate between the specified edges is therefore
Using and one Julian year gives . Thus the radiative luminosity scale is , approximately or . The black-hole mass cancels because both edge radii scale with . The specified inner edge is ; it must not be silently replaced with a different relativistic orbit radius.
For a finite viscous disc, the precise radiated power also depends on boundary work: torque transport of energy in a steady accretion disk means that orbital binding-energy release need not all be radiated inside the two stated edges. If the intended model is a radiatively efficient steady Keplerian accretion disk with a zero-torque inner boundary condition, integrate the one-face flux over both faces. It gives
The difference from is the outward torque power at the outer edge. The question supplies neither boundary torque nor the external material-supply energy budget, so the two calculations state their assumptions explicitly: about W is the binding-energy estimate, while W is the finite-disc zero-inner-torque dissipation result. Both neglect radiation advection, winds and relativistic corrections as required for this model.[ref]
A stable sinusoidal radial-velocity semi-amplitude can be stellar reflex motion caused by a companion on a nearly circular orbit. The amplitude is compatible with a giant planet; a low-mass star or brown dwarf viewed nearly face-on can mimic the same projected motion, but requires an increasingly unlikely orientation. Starspots and magnetically suppressed convection can produce a rotational stellar activity signal on a ten-day timescale; contemporaneous photometry, activity indicators and stellar spectral-line bisector and other profile changes distinguish these from a clean common Doppler shift. Solar-like stellar oscillations ordinarily have much shorter periods, so a coherent ten-day oscillation of an otherwise ordinary solar-type dwarf is less plausible. Instrumental drift, sampling aliases and blended spectra also need checking. A sinusoid alone does not establish a planetary origin, and the data do not assign numerical probabilities to these alternatives.
For the minimum companion mass from radial velocity, take a circular two-body orbit and initially . The exact relation and its low-mass limit are
With ,
Here is the relative orbital semi-major axis. Nonzero orbital eccentricity multiplies the inferred by at fixed measured , but the sinusoidal model motivates here.
For the radial-velocity sensitivity limit at fixed planet mass, use and interpret the stated detection scale as a semi-amplitude threshold . The largest distance occurs for :
For an inclination , this boundary is smaller by . The result assumes adequate time coverage and neglects the stipulated stellar noise. If the detection scale instead denotes a full peak-to-peak excursion, set to half that scale; then increases fourfold.
The stellar-mass scaling of radial-velocity planet estimates changes these inferences because Kepler's third law and the stellar reflex response depend on the stellar mass. At fixed , and , whereas the fixed-mass detectability boundary scales as . Therefore
The observed period and amplitude are unchanged. The stellar radius does not enter these leading dynamical estimates, but is crucial for the exoplanet transit inference.
For the transiting companion, assume an opaque spherical planet, negligible transit dilution and negligible limb darkening. Its radius ratio is , hence
Using the revised stellar mass and a circular orbit gives and . Interpret the duration as the total first-to-fourth-contact duration . The impact parameter from total transit duration follows from the total duration of a circular planetary transit:
This is noncentral but not grazing: . The small-planet straight-chord approximation gives ; retaining finite and orbital curvature gives the stated value. A duration measured between other contacts needs the corresponding contact formula. Unknown orbital eccentricity and argument of periapsis change the transverse speed and prevent a unique unconstrained ; the numerical estimate is conditional on the stated circular model.[ref]
Figure 1.
The measured transit duration selects a short off-centre chord when finite planet size is retained
. Source. Details.
Because the transit constrains , the RV nondetection now bounds the actual planetary mass, rather than only its projected mass. The bulk-density upper limit from a transit and an RV nondetection is
These bounds assume the instrumental threshold can be applied to a fitted signal at the transit period; no statistical confidence level is specified. Such a large radius at this low mass is inconsistent with an ordinary compact rock-and-iron body. A low-density volatile-rich planet, possibly with an extended H/He envelope, is indicated under the planetary interpretation. The observations do not uniquely determine interior composition, atmospheric abundance or formation history; blended eclipses and activity still require validation. There is no luminosity or effective temperature supplied for the host, so a numerical planetary equilibrium temperature cannot be inferred from these data alone.
The numerical conversions above use externally documented , , , , , and . They are conversion constants, not additional measured system parameters.[ref][ref]
Adopt circular orbits, negligible planetary masses, spherical bodies, uniform stellar brightness and no transit dilution. Interpret the measured duration as the total duration of a circular planetary transit, from first to fourth contact, . These conventions matter especially for the off-centre crossing: finite planet size cannot automatically be ignored there. Unknown orbital eccentricity or limb darkening would require further information to infer unique stellar radii.
First convert the apparent magnitude difference into a fractional exoplanet transit depth, rather than treating a milli-magnitude as a part in a thousand:
Both crossings are non-grazing because their transit impact parameters satisfy , so the full projected planetary disc can block the star. Equal depths determine equal radius ratios, not equal absolute planet radii.
For orbital inclination measured from the sky plane, . Choosing mid-transit as orbital phase zero, the projected separation obeys
At first and fourth contact, . The full crossing duration therefore is
For the numerical scale use solar radius and nominal solar mass parameter , with one day equal to . These are stated conversion choices, not additional measurements of the two systems.[ref]
Kepler's third law gives . For a solar-mass, solar-radius host, this gives at five days and at two days. Substitution in the contact formula yields
The transit duration of a small planet instead takes both and , giving . It predicts and , respectively. The difference is appreciable for B because the fractional change in its crossing chord is larger. Durations defined by mid-ingress to mid-egress or by second to third contact should not be silently equated with .
Figure 1.
Finite planet size lengthens the contact-to-contact chord, particularly off centre
. Details.
The original diagram uses a straight projected track to display the contact geometry; the duration formula above retains orbital curvature.
To infer the actual stellar radii, invert the contact formula before applying the supplied mass-radius scaling. Let and ; then
The measured durations imply and . Normalize the stipulated constant mass-radius ratio to the Sun, so . Kepler's third law then gives
Thus the stellar radius from transit duration and a linear mass-radius relation is
The same assumed scaling gives and . A's longer crossing requires a lower mean stellar mass density, while B's shorter full-contact crossing requires a higher one. They are not both solar replicas under this model. Their planets have and .
For a consistently applied point-planet approximation, the duration scales as and the inferred radii would instead be , . The latter is close to solar. Consequently the assertion of a substantial non-solar radius for B depends on which duration convention and finite-size approximation are intended. The contact-based values above use the measured depth rather than discarding it. If eccentric orbits are allowed, the transit speed introduces an additional factor depending on eccentricity and orientation, so the two durations alone cannot establish those radii uniquely.
For an outer coplanar planet around B, the common orbital inclination gives . A point-planet centre must cross the stellar disc, requiring . The usual coplanar transit limit is therefore
This is a limiting value approached from below; a tangent path has no duration. For an outer planet of finite, unspecified radius ratio , any overlap instead requires , giving . If it had the same ratio , that geometric limit would be days; requiring full overlap would give days. Because neither its radius nor a photometric detection threshold is supplied, no unique finite-size practical maximum period can be inferred beyond this conditional expression. The stellar radius cancels from the point-planet period-ratio bound.
For B's inner planet, a reflected-light planetary phase curve has
where is geometric albedo and is phase angle. Near exoplanet secondary eclipse the visible hemisphere is almost fully illuminated. Using the standard full-phase approximation and the contact-based ,
The exact lost fraction of combined stellar and planetary light is , which changes the result negligibly at this contrast. The corresponding magnitude depth is about milli-magnitudes. The circular orbit has the same impact distance at eclipse as at transit, and , so the star can hide the entire planetary disc.
The inclination inferred for B is , putting closest approach to full phase at . A fully precise prediction would therefore also require the reflected-light phase function, which the question does not specify; the boxed value uses . The consistent point-planet duration approximation yields about instead, through its different inferred . It is important to keep the geometry used for the radius inference and the reflected-light calculation consistent.
Finally, the stellar radial-velocity semi-amplitude for a low-mass companion on a circular orbit is
The two stars' equal apparent magnitudes make a comparison at similar photon statistics plausible, though spectral-line widths, stellar variability and observing cadence can still differ. With the contact-based stellar estimates,
If the planets have comparable mean densities, A is easier: their equal depth then gives , so . The larger planetary mass outweighs A's longer orbital period and larger host mass. For equal planetary masses, B is easier by a factor in amplitude. The supplied transit data give planet radii, not masses, so the ranking is conditional on a planetary mass-radius assumption; a unique unconditional ranking would invent information.
Sub-Keplerian rotation means that the gas rotates more slowly than a pressure-free particle on a circular orbit at the same cylindrical radius. Let . Neglecting radial acceleration, radial force balance in the gas is
An outwardly decreasing pressure, , supplies part of the inward-gravity balance through an outward pressure force. Less centrifugal support is needed, so . Using , define the pressure-gradient parameter of a protoplanetary disk by
Then
For a thin disk and a smooth radial pressure profile, is small, of order . The factor of two is part of this explicit convention; changing the convention for changes the drift formula correspondingly.
A solid particle feels essentially no continuum gas-pressure support. Left alone, it would orbit at , so in the gas it feels a headwind. Drag removes angular momentum, causing radial drift of solids in a gas disk. Define Stokes number using the aerodynamic stopping time . In a dilute test-particle limit with linear drag, negligible radial gas flow and small departures from a circular Keplerian orbit, let . The steady local radial and azimuthal balances are
Solving gives
The sign is inward for the usual negative pressure gradient. Very small particles, , closely follow the gas and have a small relative drift . Very weakly coupled large bodies, , also drift slowly, with . The fastest pressure-driven inward drift occurs at , where . The corresponding physical grain or aggregate size depends on gas density, orbital radius and the drag regime; “metre-size barrier” is a conventional description of some disk models, not a universal size.
Figure 1.
Pressure support creates a headwind, with maximum inward solid drift at unit Stokes number
. Details.
The drift timescale in this model is
It can be much shorter than the disk gas lifetime, threatening loss of the material needed to make planetesimals. Growth through the vulnerable coupling range must compete with this loss, while differential drift generates collisions that may stick, bounce or fragment. Radial redistribution also changes solid-to-gas ratios and moves material across condensation fronts. These are consequences for planet formation, rather than merely a small correction to an orbital frequency.
The same mechanism can help concentrate solids. At a local radial pressure maximum, , so this component of drift vanishes. On the inside, a positive pressure gradient gives outward drift; on the outside, a negative gradient gives inward drift. Particles therefore converge toward a dust trap. Turbulent diffusion and inward gas advection can limit retention. With a nonzero radial gas speed , the test-particle result generalizes to , so gas advection need not vanish even at a pressure maximum.
Finally, if local solid loading becomes large, backreaction on the gas can no longer be neglected. Coupled gas-particle drift can support streaming instability and strong solid concentrations, potentially enabling gravitational collapse into planetesimals. The simple test-particle drift formula is a useful starting point, but not a description of that nonlinear high-loading regime. Pressure support creates both a drift obstacle to growth and opportunities for concentration where the pressure structure changes.
Radial epicyclic motion. At fixed specific angular momentum, use and . A circular orbit at obeys . Expanding gives , where
For the radial epicyclic frequency describes simple harmonic motion. The angular correction follows from conservation of angular momentum, , and describes an epicyclic motion around the guiding circular orbit.
Stabilizing terms. In the given local fluid dispersion relation, is rotational support, important at long wavelength; is random-motion or pressure support, important at short wavelength. The term is destabilizing self-gravity. Although the wording refers to stellar dispersion, this algebra uses the supplied fluid approximation; a collisionless stellar treatment has a reduction factor and a different numerical normalization.[ref]
Axisymmetric threshold. For , set . Completing the square gives
Thus gives local axisymmetric oscillatory stability; gives a growing band, with
At , the minimizing mode has zero frequency and is marginal, not a strictly restoring oscillation. Local axisymmetric stability does not establish stability to every global nonaxisymmetric mode.
long-wave and short-wave branches at marginal Toomre stability. For and ,
The long-wavelength branch, using the minus sign, requires
These endpoints are the inner and outer Lindblad resonance conditions. The stated dispersion relation does not restrict every wave solution to that interval. Its plus-sign short-wavelength branch has even outside it: for example, gives and satisfies the equation. Additional branch selection or resonant absorption physics would be needed for an absolute confinement claim. No such condition is supplied.
Toomre stability radius of a flat-rotation disk. With , . Using ,
The supplied local model therefore predicts stable radii below about , marginality at that radius, and axisymmetric instability farther out while its assumed constant surface density and thin-disk approximation remain applicable. An infinite constant-density, flat-rotation disk is an idealization, not a global self-consistent galaxy model.
Figure 1.
Original marginal-disk dispersion branches and radial Toomre parameter, showing the long-wave restriction and the surviving short-wave branch
. Details.