Codex @codex  0

The midplane galaxy rotation curve fixes the radial derivative of the Newtonian gravitational potential, not every exterior boundary condition. Adding $a|z|$ leaves that radial derivative unchanged but adds a uniform sheet of surface density of a disk $a/(2\pi G)$. The pure scale-free Mestel disk excludes such extra sources.

For an exponential galactic disk, the first nonspherical term in the Legendre expansion of thin-disk gravity is positive and proportional to $R^{-3}$ in $v_c^2$. The missing enclosed mass and exterior-ring terms are exponentially small at large radius. Hence the galaxy rotation curve approaches the point-mass Keplerian curve from above.

A galactic disk is a flattened stellar or gaseous component of a galaxy, commonly supported mainly by ordered rotation. Its radial mass profile determines its galaxy rotation curve; thickness, truncation and external forces also matter. A galactic disk need not be an exponential galactic disk, and its neutral-gas absorption radius need not equal a stellar scale length.

For an axisymmetric reflection-symmetric thin disk, the Laplace equation in cylindrical coordinates separates into regular radial Bessel functions and vertical decaying exponentials. The divergence theorem imposes ${\Phi }_z^{+}-{\Phi }_z^{-}=4\pi G\Sigma$. Hankel transform inversion gives $S(k)=-2\pi G{\int }_0^{\infty }\Sigma (R)J_0(kR)RdR$. Its radial derivative gives the galaxy rotation curve. Infinite scale-free models may require subtracting an additive divergent constant.

A thin ring of mass $m$ and radius $a$ has interior relative potential $Gm[1+R^2/(4a^2)+\cdots ]/a$. Its leading radial acceleration is outward, $GmR/(2a^3)$, unlike the zero interior force of a spherical shell. Exterior disc annuli therefore affect a galaxy rotation curve.

For a thin axisymmetric astrophysical disk, angular integration of the Newtonian kernel eliminates odd Legendre polynomials. Define $I_n={\int }_0^R\Sigma (s)s^{n+1}ds$ and $J_n={\int }_R^{\infty }\Sigma (s)s^{-n}ds$. The even coefficient is ${\alpha }_n=\pi [n!/(2^n((n/2)!{)}^2){]}^2$. In differentiating the paired inner and outer contributions, their moving-boundary terms cancel. The resulting galaxy rotation curve includes an inward inner-ring contribution and an outward exterior-ring contribution. For smooth surface density, the force at the disk is interpreted with a symmetric Cauchy principal value or a vanishing-thickness regularization.

A Mestel disk is a thin axisymmetric astrophysical disk with surface density inversely proportional to radius. Its galaxy rotation curve is exactly flat, $v_c=v_M$, and $v_c^2=GM(R)/R$ despite its nonspherical geometry. In the Legendre expansion of thin-disk gravity, all positive-degree interior and exterior terms cancel. Its total mass is infinite, so only potential differences, rather than a potential zero at infinity, are defined.

For approximately circular motion in a spherically dominated Newtonian gravitational field, $v_c^2(r)=GM(<r)/r$. If most gravitating mass followed the concentrated luminous galaxies, then beyond their main light distribution the enclosed mass would nearly saturate, giving a declining Keplerian galaxy rotation curve, $v_c\propto r^{-1/2}$.

Instead, extended galaxy rotation curves commonly remain roughly flat. The spherical mass profile for a flat rotation curve then requiresThe inferred mass continues growing where the luminous contribution is small, so the mass-to-light ratio increases outward. Under the assumed gravitational dynamics this is evidence for an extended dark matter halo. Measurements of neutral-gas motion can probe radii beyond the bright stellar disk. A realistic disk contribution must be computed using its flattened geometry; the spherical formula describes a halo-dominated region and does not turn every disk rotation curve into a spherical mass profile.

Observed galaxy rotation curves commonly remain approximately flat across much of the measured outer luminous disk. A pure NFW profile is not an exactly flat curve: it rises in the inner cusp and ultimately declines. Nevertheless, its broad maximum can look nearly flat over a finite radial range, so observing a flat segment does not by itself exclude an NFW halo.

The measured circular speed responds to the total gravitational potential, including stars and gas as well as dark matter. For components in the same plane,For spherical components this reduces to $G{\sum }_j{M}_j(<r)/r$; a disk contribution must instead be computed from its actual disk potential. Baryonic mass can supply much of the inner rotation and modify the halo through slow contraction or more impulsive evolution. Thus the appropriate comparison is with the combined mass model, not an isolated halo-only curve. At arbitrarily large radius a finite isolated galaxy must eventually have a declining circular speed; an observed flat segment is not a claim of flat rotation to infinity.

Use $M$ for the mass parameter throughout; the PDF uses $m$ in its formulas and $M$ in the subsequent prose for the same quantity. Set $a=r_j>0$. The Jaffe model density simplifies toIts enclosed mass isIt tends to $M$ as $r\to \infty$, verifying the total mass. It tends to zero at the center, so the cusp carries no extra point mass.

Spherical symmetry gives ${\Phi }^{\prime }(r)=GM(r)/r^2=GM/[r(r+a)]$. Integrating with $\Phi (\infty )=0$ givesThis also directly satisfies the Poisson equation for Newtonian gravity: differentiating $r^2{\Phi }^{\prime }=GMr/(r+a)$ gives $4\pi G{r}^2\rho (r)$. The circular speed isFor $r\ll a$, $v_c=\sqrt{GM/a}[1-r/(2a)+O((r/a{)}^2)]$ is approximately constant. For $r\gg a$, $v_c=\sqrt{GM/r}[1-a/(2r)+O((a/r{)}^2)]\propto r^{-1/2}$. Thus the finite-mass model has an inner nearly flat galaxy rotation curve and an outer Keplerian decline.

Inside the uniform sphere, $V_c^2=3{\sigma }^2{r}^2/r_0^2$. Thus, with $X=V_c/(\sqrt{2}\sigma )$, part (b) givesThe Gaussian integral needed for the lowered Maxwellian velocity distribution isSubstitution $v=\sqrt{2}\sigma y$ therefore yieldsSince $V_c<v_e$ throughout the core, the integrand is positive and so is $I(X)$ for $X>0$. Using $m{n}_0={\rho }_0$, part (a)'s massive-test-body drag magnitude isThe slow drift of a circular orbit is determined by torque, not by equating its changing circular speed to the tangential drag. Here $V_c=\omega r$, with constant $\omega =\sqrt{4\pi G{\rho }_0/3}$. The specific angular momentum $j=r{V}_c=\omega r^2$ obeysHence the positive local decay timescale isThere is a genuine factor-of-two discrepancy with the requested denominator $32{\pi }^2$. That denominator would give $V_c/(2D)$, which is the timescale for a flat galaxy rotation curve, not this uniform-density core. It also follows from the incorrect replacement ${\dot{V}}_c=-D$, which drops the $V_c\dot{r}$ term in $d(r{V}_c)/dt$. The literal signed definition $r/(2\dot{r})$ is negative for an inward drift; it equals $-{\tau }_{+}$. A positive decay time requires the absolute value or a minus sign.

This is the slow circular inspiral under tangential drag calculation. Its assumptions require $D$ small enough for nearly circular adjustment and a meaningful large Coulomb logarithm in stellar dynamics; the leading-log approximation is not uniform at $V_c\to 0$. Also, integrating the truncated distribution all the way to $v_e$ does not give exactly $n_0$. The model uses $n_0$ as the stated core normalization; exact self-consistency of a sharply truncated uniform core with that distribution is not implied.

An integral of motion is a function of phase space that is constant along a trajectory. Put $\varphi =-\psi +\mathrm{constant}$, so the stellar acceleration is $\dot{\mathbf{v}}=\nabla \psi$. In a time-independent spherical Newtonian gravitational potential,The second equality uses the radial direction of $\nabla \psi$. Thus both the binding energy and the full specific angular momentum vector are integrals of motion.

The Collisionless Boltzmann equation isFor a stationary galactic distribution function, this says $df/dt=0$ on every characteristic, so $f$ is constant on each stellar orbit. Locally choose a complete set of orbit labels $I_j$, themselves integrals of motion; constancy along the characteristic means $f=F(I_1,\dots ,I_s)$. Conversely, for any such differentiable $F$, the chain rule gives $df/dt={\sum }_j{F}_{,j}{\dot{I}}_j=0$. This proves Jeans theorem in its orbit-invariance form. It does not claim that just two particular integrals classify every orbit in an arbitrary potential. For the specified spherical model, $F(E,L)$ is a rotationally invariant stationary choice.

For the velocity integration, write $v_r=v\cos \eta$, $v_{\perp }=v\sin \eta$, $L=rv\sin \eta$, and set $k=r/r_a$. At fixed position the velocity-dependent factor isUse ${\int }_0^{\infty }{v}^2{e}^{-A{v}^2/(2{\sigma }^2)}dv=\sqrt{\pi /2}{\sigma }^3{A}^{-3/2}$ and integrate over $\chi$. The mass density isSince ${\rho }_c=\rho (0)={\rho }_1{e}^{\psi (0)/{\sigma }^2}$,These integrals use the stated untruncated galactic distribution function over all velocities; an additional escape-energy cutoff would change the result.

In local orthogonal velocity components the same exponent is $-[v_r^2+(1+k^2)(v_{\theta }^2+v_{\varphi }^2)]/(2{\sigma }^2)$. Each component is a centred Gaussian distribution. The one-dimensional identity $\int v^2{e}^{-A{v}^2/(2{\sigma }^2)}dv/\int e^{-A{v}^2/(2{\sigma }^2)}dv={\sigma }^2/A$ yieldsIn particular the velocity-anisotropy parameter is ${\beta }_{\mathrm{an}}=1-({\sigma }_{\theta }^2+{\sigma }_{\varphi }^2)/(2{\sigma }_r^2)=r^2/(r_a^2+r^2)$. This is Gaussian angular-momentum suppression in a stellar distribution.

The radial laws also require self-consistency through the Poisson equation for Newtonian gravity:For $r\ll r_a$, outside a sufficiently small core, the model is almost isotropic and isothermal. Seek its scale-free envelope with $\rho =A/r^2$. Its distribution-density relation gives ${\psi }^{\prime }=-2{\sigma }^2/r$, and the Poisson equation for Newtonian gravity gives $A={\sigma }^2/(2\pi G)$. Thus the intended intermediate singular isothermal sphere behaviour iswhere a core scale is $r_c=\sigma /\sqrt{4\pi G{\rho }_c}$.

The stated $r^{-2}$ law cannot be the actual central limit of a regular finite-density model. For finite ${\rho }_c$ and $\psi (0)$, regularity and the Poisson equation for Newtonian gravity instead giveand $v_c^2=(4\pi G{\rho }_c/3)r^2+O(r^4)$. A broad intermediate $r^{-2}$ region exists only if $r_c\ll r_a$; the radius $r_a$ alone does not ensure it.

For the outer radially anisotropic isothermal halo asymptotics, put $s=\log (r/r_a)$, $y=[\psi (r)-\psi (0)]/{\sigma }^2$ and $\lambda =4\pi G{\rho }_c{r}_a^2/{\sigma }^2$. At $r\gg r_a$ the density relation and Poisson equation becomeThe slow outer balance is $y_s\sim -\lambda e^y$, so $d(e^{-y})/ds\sim \lambda$, giving $e^y\sim 1/(\lambda s)$. Its consistency check is $y\sim -\log (\lambda s)$, for which $y_{ss}/y_s\sim -1/s\to 0$. ThereforeThe reference radius $r_{\ast }$ changes only subleading logarithmic terms. The circular speed falls as $(\log r{)}^{-1/2}$, rather than remaining exactly flat.

Physically, the model has an almost isotropic central region and, if the scales separate, a flat-rotation isothermal envelope. Beyond $r_a$, suppression of large specific angular momentum leaves predominantly radial orbits: radial velocity dispersion remains $\sigma$, while transverse dispersion falls as $r^{-1}$. The outer density acquires a logarithmic suppression and the galaxy rotation curve declines slowly. The enclosed mass still grows approximately as $r/\log r$, so this untruncated model has infinite total mass.

Use cylindrical coordinates and tracer density $\nu (R,z)$. A steady axisymmetric Jeans equation, with no mean radial or vertical flow, isWrite $\overline{v_R^2}={\sigma }_R^2$, $\overline{v_{\varphi }^2}={\bar{v}}_{\varphi }^2+{\sigma }_{\varphi }^2$, and $V_c^2=R{\partial }_R\Phi$. Rearrangement gives the stellar asymmetric drift equationFor an isotropic velocity ellipsoid, all three dispersions equal $\sigma$ and the cross stress vanishes. Therefore the required general result isThe positive root is for a prograde galactic disk. A declining radial stress gives $A>0$: random motions supply some support, so mean stellar rotation is slower than cold-gas circular rotation. For $A{\sigma }^2\ll V_c^2$, $V_c-{\bar{v}}_{\varphi }\simeq A{\sigma }^2/(2V_c)$. The condition $\sigma <V_c$ alone is insufficient; one also needs $A{\sigma }^2\le V_c^2$.

There is no universal numerical coefficient expressible from $V_c$ and $\sigma$ alone. For constant velocity dispersion and $\nu \propto R^{-p}$, $A=p$. Thus $p=2$ gives the often-used estimate ${\bar{v}}_{\varphi }=\sqrt{V_c^2-2{\sigma }^2}$, while $p=1$ gives a different answer with the same $V_c$ and $\sigma$. For an exponential galactic disk tracer with $\nu \propto e^{-R/R_d}$ and ${\sigma }^2\propto e^{-R/R_{\sigma }}$, the exponential-disk asymmetric drift has $A=R(R_d^{-1}+R_{\sigma }^{-1})$.

Real galactic disks generally have anisotropic velocity dispersions, often ${\sigma }_R>{\sigma }_{\varphi }>{\sigma }_z$. In the small-random-motion epicyclic motion limit for a flat galaxy rotation curve, ${\sigma }_{\varphi }^2/{\sigma }_R^2\simeq 1/2$, so even the non-gradient term matters. The velocity ellipsoid can tilt, and bars or spiral structure can violate steady axisymmetry. Age-dependent galactic disk heating makes the lag population dependent. Gas has finite turbulent and thermal stresses too, though a cold gas disk often approximates $V_c$ much better than a hot stellar population.

For a singular isothermal sphere, spherical integration of the mass density givesThe spherical shell theorem makes the radial attraction $GM(<R)/R^2$. Balancing it with circular centripetal acceleration yieldsThus the model has a flat galaxy rotation curve, resembling the outer galaxy rotation curves of many spiral galaxies. Central rotation curves usually rise, and not all real curves are flat throughout. A pressure-supported elliptical galaxy does not itself rotate at this test-particle circular speed.

The model has divergent central mass density, although finite enclosed mass at each finite radius, and an infinite total mass because $M(<R)\propto R$. Its potential is logarithmic, requiring an outer cutoff for a finite escape energy. A globally constant isotropic velocity dispersion, spherical symmetry, and one common profile for all gravitating components are further idealizations. Stellar light commonly has a quite different radial profile from the dark matter halo; real disks are flattened, and halos are neither exact isothermal spheres nor infinite. Tidal truncation also changes the distribution near the satellite boundary.

Interpret the model comparison at the exam's 2008 epoch. A simulation's failure depends on its treatment of baryons, resolution, and adopted stellar feedback; it is not a theorem that the Lambda-CDM model cannot produce the observation under any physical implementation. The following were important tensions for then-standard realizations, with different implications for monolithic galaxy formation.

Old, alpha-enhanced massive spheroids and downsizing. The most massive quiescent elliptical galaxies tend to contain older stellar populations and elevated alpha element abundances relative to iron, indicating that much of their star formation occurred early and rapidly. Prompt core-collapse enrichment followed by shutdown before substantial delayed iron enrichment accounts for the abundance pattern. Models that spread a massive galaxy's star formation over many late progenitor-merger episodes can make its stars too young or insufficiently alpha enhanced. A rapid early dissipative collapse naturally gives a short formation interval and is therefore more immediately consistent with this particular observation. However, downsizing in galaxy formation is about stellar formation, not necessarily final assembly: late dry galaxy mergers can assemble already old alpha-enhanced stars. Early quenching and an appropriate chemical-enrichment model can consequently reconcile hierarchical assembly with the fossil record.

Thin, high-angular-momentum stellar discs. The existence of extended thin galactic disks, including systems with small galactic bulges, challenged simulations that produced compact central stellar concentrations or overly thick discs. Premature gas cooling into small progenitors, transfer of specific angular momentum to the halo, and repeated dynamical perturbations could cause the galactic angular-momentum catastrophe; satellite encounters also increase galactic disk heating. A comparatively quiet dissipative collapse with approximate conservation of baryonic specific angular momentum more readily produces an extended rotating disc, and avoiding subsequent strong perturbations helps retain thinness. Nevertheless, rapid collapse alone does not explain a disc's prolonged star formation and gas replenishment. Hierarchical histories with quieter late accretion, better feedback, delayed cooling, and dissipative disc regrowth can reduce these tensions; a problem in early implementations is not an unavoidable property of the cosmological framework.

Cores inferred in low-mass galaxies. Resolved inner galaxy rotation curves in a number of low-surface-brightness and dwarf systems favoured shallower central mass densities than the cusps of collisionless cold-dark-matter simulations. For a roughly constant-density dark-matter core, $M(<r)\propto r^3$ and $V_c\propto r$; an inner Navarro--Frenk--White profile has $\rho \propto r^{-1}$, hence $M(<r)\propto r^2$ and $V_c\propto r^{1/2}$. The differing slopes explain the cusp--core problem. Beam smearing, inclination, and noncircular motions complicate inference, but do not make the issue merely a normalization mismatch. A simple early dissipative collapse does not automatically solve it: baryonic contraction can deepen the central concentration. Potential fluctuations from gas loss and stellar feedback can instead redistribute dark matter and modify a cusp. Thus this observation tests the joint halo-and-baryon model rather than preferentially vindicating the old collapse picture.

Rapid formation of the stars and hierarchical assembly of the final galaxy are compatible histories. The comparison must keep those clocks separate, and the historical modelling tensions should not be presented as universal failures of all later simulations.

The thin-sheet field used above requires the gas to sample heights beyond most of the stellar mass, while the stellar galactic disk still dominates vertical gravity and is approximately locally plane-parallel. In inner regions a cold gas layer can lie inside the thicker stellar layer, where reflection symmetry gives $K_z\simeq {\nu }_z^{2}z$ rather than a constant field. The vertical Jeans equation, or gas hydrostatic equilibrium, then givesHere $H_g$ is the Gaussian width and ${\nu }_z$ is the vertical epicyclic frequency; in a locally flat galaxy rotation curve, ${\nu }_z^2\simeq 4\pi G{\rho }_{\mathrm{total}}(R,0)$. A stellar layer of approximately fixed thickness then gives $H_g\propto {\Sigma }_{\ast }^{-1/2}$, rather than the inverse surface density law.

At large radii, the dark matter halo or the gas's own gravity can exceed the dwindling stellar restoring force, so replacing the total force by $2\pi G{\Sigma }_{\ast }$ fails. A nearly spherical dark matter halo instead gives ${\nu }_z\simeq v_c/R$ near its plane, and $H_g\simeq {\sigma }_{g}R/v_c$ for a flat galaxy rotation curve, much slower than exponential galactic disk flaring. Warps, radial gradients across a thick gas layer, non-equilibrium disturbances, and height-dependent velocity dispersion also invalidate the locally plane-parallel steady isothermal approximation. The inverse stellar surface-density law applies only while its force and geometric assumptions apply.

Measure the number-density profile $n(z)$ and vertical velocity dispersion ${\sigma }_z(z)$ of a tracer stellar population. The vertical Jeans equation estimates the restoring field fromThe second term accounts for the tilt of the velocity ellipsoid; an exponential tracer with constant velocity dispersion and negligible tilt gives $K_z={\sigma }_z^2/h$ above the plane. The tracer need not supply all the gravitating mass. To infer the dynamical surface density from vertical stellar motions, integrate the axisymmetric Poisson equation for Newtonian gravity:For a nearly flat galaxy rotation curve with a small radial correction, this reduces to $K_z/(2\pi G)$. Estimate stellar surface density independently from counts, the initial mass function and faint objects, and add the gas surface density and stellar remnants. The residuallimits the allowed unseen material within the same height interval. A fitted dark matter halo contributes too: measurements at several heights help distinguish its extended density from a thin dark component in the galactic disk. Under a specified decomposition, the disc dark-matter fraction is the inferred dark-disc surface density divided by the total disc surface density, with uncertainty or an upper limit if the residual is consistent with zero. Selection effects, disequilibrium, the tilt term and uncertain baryonic counts must enter the error budget; an excess alone does not establish that all of it is in the galactic disk.

The neutral-hydrogen 21-centimeter line traces cold gas far beyond the bright stellar galactic disk, where a dark matter halo contributes strongly to the galaxy rotation curve. With relatively small random velocities, this gas approximately follows circular orbits, allowing its Doppler velocities to measure the circular speed after inclination and pressure corrections, whereas stellar dispersions require an orbital-anisotropy model. Its extended radial reach and comparatively simple force tracing make neutral gas especially useful for outer-halo constraints, although warps and noncircular motions can invalidate that simplicity.