Codex @codex  0

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.

An approximately extended component surrounding a galactic disk, distinguished according to whether one means the stellar halo, gaseous halo or dark matter halo. These components have different constituents and scales.

A gas-rich galaxy merger involves a substantial reservoir of gas capable of radiative cooling. Dissipation and gravitational torques permit gas inflow and enhanced star formation, unlike the nearly collisionless stellar evolution in a dry galaxy merger. A surviving or replenished gas reservoir can subsequently form a galactic disk.

The $k_R=-1$ orbital resonance of an $m$-fold perturbation of a galactic disk. For near-circular stellar orbits, ${\Omega }_R$ is the radial epicyclic frequency. The harmonic changes radial action and azimuthal angular momentum in the ratio $\Delta J_R/\Delta J_{\phi }=-1/m$, so an orbit losing angular momentum gains random radial motion.

A major galaxy merger joins galaxies with comparable masses. Rapid changes in the gravitational potential can redistribute stellar orbital energies through violent relaxation, substantially disturbing preexisting galactic disks. The outcome depends on gas content, orbit and subsequent accretion; a major merger does not universally prohibit disk rebuilding.

An open cluster is a relatively low-mass group of stars formed together, typically in a galactic disk. Two-body relaxation, stellar evaporation, and external tidal disturbances disperse many open clusters on much shorter times than the lifetimes of their low-mass member stars.

The $k_R=+1$ orbital resonance of an $m$-fold perturbation of a galactic disk. For near-circular stellar orbits, ${\Omega }_R$ is the radial epicyclic frequency. The harmonic changes radial action and azimuthal angular momentum in the ratio $\Delta J_R/\Delta J_{\phi }=1/m$, so an orbit gaining angular momentum gains random radial motion.

Let the uniform thin galactic disk have outer radius $R_d$ and surface mass density $\Sigma =M/(\pi R_d^2)$, where $M$ is the mass calculated in a. If the redistribution retains the stated outer radius, set $R_d=16\mathrm{kpc}$; otherwise $M$ alone does not determine both $R_d$ and $\Sigma$.

A ring of radius $s$ and width $ds$ has mass $dM=2\pi \Sigma sds$. By axial symmetry its transverse gravitational force components cancel. At signed height $z$, its vertical force on a star of mass $m$ isIntegrating the rings gives, for $z\ne 0$,ThusFor a positive distance above the plane, this is $\mathbf{F}=-2\pi Gm\Sigma [1-z/\sqrt{z^2+R_d^2}]\hat{\mathbf{z}}$, directed toward the disk. With $R_d=16\mathrm{kpc}$, the inferred $\Sigma \simeq 178.3M_{\odot }\mathrm{p}{\mathrm{c}}^{-2}$ sets its normalization.

The far-field expansion gives $F_z\simeq -GMm\mathrm{sgn}(z)/z^2$, as required for a point-like total mass. Close to either face, $F_z\to \mp 2\pi Gm\Sigma$, the infinite-sheet limit. The ideal zero-thickness model has a jump in gravitational field at $z=0$; the formula is understood through its one-sided limits. Its symmetric mean is zero on the plane, and a disk of finite thickness resolves the discontinuity.

The restriction $k_{\phi }=\pm m$ describes an $m$-fold rotating perturbation of a galactic disk, for example a transient spiral density wave or, for $m=2$, a galactic bar. Its pattern speed is ${\Omega }_{\mathrm{p}}$, its amplitude peaks at $t_{\mathrm{peak}}$, and its Gaussian function envelope has duration of order $\tau$. The coefficients $u_{\mathbf{k}}(\mathbf{J})$ describe its strength along different stellar orbits. A real perturbing potential requires $u_{-\mathbf{k}}=u_{\mathbf{k}}^{\ast }$.

Write $g(t)=\exp [-(t-t_{\mathrm{peak}}{)}^2/(2{\tau }^2)]$. The product entering c isThe phase therefore measures detuning from the rotating pattern:Substitution givesThe strongest long-lived response occurs near an orbital resonance, ${\omega }_{\mathbf{k}}=0$, where successive force contributions do not rapidly cancel in phase.

In an axisymmetric cylindrical coordinate system, the steady Collisionless Boltzmann equation isLet $\nu =\int f{d}^{3}v$ and $\overline{A}={\nu }^{-1}\int Af{d}^{3}v$. Multiply by $v_R$ and integrate over velocities. The spatial terms give ${\partial }_R(\nu \overline{v_R^2})+{\partial }_z(\nu \overline{v_R{v}_z})$. Integration by parts in $v_R$ supplies $\nu {\Phi }_R-\nu \overline{v_{\varphi }^2}/R$, and the $v_{\varphi }$ derivative supplies $\nu \overline{v_R^2}/R$. The $v_z$ force term vanishes because $v_R$ does not depend on $v_z$. This assumes sufficient decay of the galactic distribution function to discard all velocity boundary terms. The radial Jeans equation is consequentlyReflection symmetry about $z=0$ makes $\nu$ even and $\overline{v_R{v}_z}$ odd. At the midplane the latter vanishes, and ${\nu }^{-1}{\partial }_z(\nu \overline{v_R{v}_z})={\partial }_z\overline{v_R{v}_z}$. Multiplying by $R/\nu$ provesAll quadratic velocities here are population moments, not individual stellar velocities. In particular, symmetry does not make the midplane tilt contribution to asymmetric drift disappear: an odd mixed moment can have a nonzero derivative at zero.

For application to the Milky Way, assume negligible mean radial and vertical flows, and define the circular speed by $v_c^2=R{\Phi }_R(R,0)$. Write $\overline{v_R^2}={\sigma }_R^2$ and $\overline{v_{\varphi }^2}={\overline{v}}_{\varphi }^2+{\sigma }_{\varphi }^2$. Rearranging gives the exact stress-support relationThus the stellar asymmetric drift $v_a=v_c-{\overline{v}}_{\varphi }$ is the right-hand side divided by $v_c+{\overline{v}}_{\varphi }$. For $v_a\ll v_c$, this denominator is approximately $2v_c$. If $\nu \propto e^{-R/h_{\nu }}$ and ${\sigma }_R^2\propto e^{-R/h_{\sigma }}$ locally, the bracket becomes $R/h_{\nu }+R/h_{\sigma }-1+{\sigma }_{\varphi }^2/{\sigma }_R^2$, making the density and dispersion-gradient contributions explicit.

A radially declining random-motion stress supplies some of the force needed to support a galactic disk, so its stars need less ordered azimuthal motion than a population on circular orbits. Hotter, often older, stellar populations therefore generally rotate more slowly than cold populations in the same gravitational potential; this is the observed stellar asymmetric drift. It is a stress-support effect, not a frictional slowing of each orbit. For a velocity ellipsoid aligned approximately with spherical coordinates, $\overline{v_R{v}_z}\simeq ({\sigma }_R^2-{\sigma }_z^2)z/R$, and the midplane tilt term contributes $-({\sigma }_R^2-{\sigma }_z^2)$ rather than zero.

Local measurements of mean rotation and all components of the velocity dispersion, combined with tracer-density and dispersion gradients, can therefore constrain $v_c$ and the radial Newtonian gravitational field. A hotter sample's mean rotation must not be identified directly with the circular speed. The tracer number density need not be the gravitating mass density. Nonaxisymmetric streaming, a nonstationary disk, poorly measured radial gradients, or a neglected tilt term can bias this inference; these are assumptions of the derivation rather than extra free corrections to the circular force.

First consider a genuinely one-dimensional, time-independent collisionless stellar system with Hamiltonian $E_z=v_z^2/2+\psi (z)$. Choose $\psi (0)=0$. For bound vertical oscillations, a phase-mixed population has equal occupations of the upward and downward branches of each orbit. By Jeans theorem, its stationary distribution is a function of the conserved vertical energy, $f_z(z,v_z)=F_z(E_z)$. Equivalently, the one-dimensional Collisionless Boltzmann equation $v_z{\partial }_z{f}_z-{\psi }^{\prime }(z){\partial }_{v_z}{f}_z=0$ makes $f_z$ constant along the energy curves. Branch-dependent incoming and outgoing populations or time-dependent waves would not satisfy this assumption.

Integrate both signs of $v_z$ and change variable to $E=v_z^2/2+\psi$:The density is thus an Abel transform of the vertical-energy distribution. To invert it without simply quoting an inversion rule, apply a second integral and exchange its order:The inner integral is $\pi$, as follows by setting $\psi =E+(t-E){\sin }^2\theta$. Writing $J(E)={\int }_0^{\infty }\nu (E+s)s^{-1/2}ds$ before differentiating avoids an artificial infinite endpoint term. Under sufficient differentiability and decay,The vertical-energy Abel inversion of a stellar distribution is thereforeThis proves the normalization, including the factor from the two velocity signs. The additive origin of $\psi$ merely changes the origin of energy.

For this expression to define a physical phase-space distribution function, $\nu$ must be expressible as a single-valued function of $\psi$, the integrals and differentiated integral must converge, and the resulting $F_z$ must be nonnegative. Reflection symmetry and a potential increasing with $|z|$ provide the usual single-valued relation on either side of the disk. A sufficiently smooth decreasing $\nu (\psi )$ is a sufficient condition for nonnegativity; it is not asserted to be necessary for every positive energy distribution. If there is a finite escape energy, replace the infinite limits by that energy and impose the appropriate cutoff and boundary conditions. The displayed infinite-limit formula presumes either a confining potential or a distribution extended by zero above its allowed energy, with no omitted boundary population.

A three-dimensional galactic disk does not automatically meet these hypotheses. Exact separation $\Phi (R,z)={\Phi }_{\parallel }(R)+\psi (z)$ makes $E_z$ an integral of motion. If the full galactic distribution function factorizes into an in-plane factor and a vertical factor $F_z(E_z)$, integrating the in-plane velocities gives the required one-dimensional distribution after its normalization is absorbed into $\nu$. More generally the vertical marginal must depend only on $E_z$ in the sample being analysed. Near the solar radius, approximate separation can be useful when the sampled heights and radial excursions are small relative to the scales over which the vertical force changes with $R$. If only a slowly varying vertical action is conserved, replacing it by one fixed $E_z$ requires an additional controlled local approximation.

The population must be sufficiently stationary and phase mixed, with negligible mean vertical flow, and the data must measure a well-defined tracer population with a known selection function. Radial coupling and the tilted velocity ellipsoid must be negligible or modelled. In the full vertical Jeans equation, that coupling is $R^{-1}{\partial }_R(R\nu \overline{v_R{v}_z})$; it cannot generally be discarded merely because a mixed moment vanishes at the midplane. The tracer number density $\nu$ is not the mass density producing the potential.

To infer the local mass, measure the tracer density as a function of height and its vertical velocity distribution near the midplane. With $\psi (0)=0$, the latter directly samples $F_z(v_z^2/2)$. Insert that measured distribution into the forward integral for $\nu (\psi )$, and match it to the observed $\nu (z)$ to infer $\psi (z)$. Equivalently, for trial potentials compute $d\nu /d\psi =(d\nu /dz)/(d\psi /dz)$ away from the plane, apply the inversion, and require agreement with the measured velocities. Density data alone do not determine both an unknown potential and an unknown energy distribution.

Once the vertical force is known, use the full axisymmetric Poisson equation for Newtonian gravity, rather than identifying every vertical derivative with density in isolation:At the midplane, $v_c^2(R)=R{\partial }_R\Phi (R,0)$, soIn terms of the circular-orbit Oort constants, $A=\frac{1}{2}(v_c/R-d{v}_c/dR)$ and $B=-\frac{1}{2}(v_c/R+d{v}_c/dR)$, the radial term equals $2(B^2-A^2)$. It vanishes for a locally flat rotation curve, not for an arbitrary rotating disk. An integrated version gives the dynamical surface density from vertical stellar motions:Here the actual vertical acceleration is $-{\psi }^{\prime }(z)$. The sign of the radial correction follows directly from the Poisson equation for Newtonian gravity.

For an illustrative isothermal tracer with constant vertical velocity dispersion $s$, takeThe velocity integral gives $\nu (z)={\nu }_0{e}^{-\psi (z)/s^2}$, hence $\psi (z)=-s^2\ln [\nu (z)/{\nu }_0]$. For a locally flat rotation curve this yieldsThis example shows explicitly how a measured tracer scale and velocity scale constrain the gravitating density; using the tracer's own density as the mass source would instead impose an unjustified self-gravitating model.

The inferred density includes luminous and faint stars, gas and dark matter. Comparing it with an independently estimated baryonic census tests for extra gravitating mass. An excess need not be a dark disk: a dark matter halo also contributes to the local force. Conversely, agreement with the baryonic density locally does not exclude a halo at other radii or heights. Vertical disequilibrium, radial mixing, errors in selection or tracer populations, and neglected velocity tilt can mimic a mass discrepancy, so the one-dimensional equilibrium assumptions must be tested when interpreting such an inference.

The qualitative prediction of galactic chemical evolution is useful: later stellar populations tend to inherit more enriched gas, and old stellar populations can preserve low stellar metallicity. The closed-box metallicity distribution, however, predicts many more metal-poor long-lived dwarfs in the solar neighbourhood than are observed. This is the G-dwarf problem, also apparent with other long-lived tracers. A single well-mixed box cannot describe the different age and abundance distributions of the galactic disk, galactic bulge, and stellar halo of the Milky Way. Nor need its exponentially declining star formation rate match the history of ongoing disk star formation.

Prolonged relatively metal-poor galactic gas inflow supplies fuel while diluting the enrichment, allowing many stars to form after the earliest enrichment stage; pre-enrichment of the initial reservoir also removes a pristine tail. Galactic outflows lower the effective yield, but a constant well-mixed outflow alone still gives an exponential metallicity distribution function, and is not automatically a cure for the G-dwarf problem. Spatially varying star formation, radial transport and migration, and selection of local survivors complicate a comparison with one global box. A varying initial mass function or stellar yield changes the enrichment law. Finally, finite stellar lifetimes and Type Ia supernovae make iron enrichment delayed, so an observed iron-abundance distribution is not interchangeable with the instantaneous total-metal distribution. Gas exchange, spatially distinct populations, and delayed enrichment are essential qualifications to the closed-box prediction.

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.

Measure the logarithmic surface brightness slope in resolved radial annuli after accounting for instrumental blurring and foreground emission. Between a halo core and its outer boundary, dark-matter decay radiation predicts $d\ln I/d\ln R=-1$, whereas dark-matter annihilation radiation predicts $-3$; doubling radius therefore reduces the respective surface brightness by factors $2$ and $8$. These are quantitative templates to compare with independently mapped stellar and gaseous components.

An exponential galactic disk has $I_{\ast }(R)\propto e^{-R/R_d}$, so $d\ln I_{\ast }/d\ln R=-R/R_d$ and $I_{\ast }(2R)/I_{\ast }(R)=e^{-R/R_d}$. At $R=5R_d$ the latter is about $0.0067$, much faster decline than either halo template, and the projected emission should also be flattened along the galactic disk. A galactic bulge described by a Sérsic profile has $d\ln I_{\ast }/d\ln R=-(b_n/n)(R/R_e{)}^{1/n}$; for $n=4$ and $b_4\simeq 7.67$, this is about $-3.4$ at $10R_e$ and continues to steepen, whereas a pure halo power law keeps a fixed slope.

For a stellar halo with outer luminosity density proportional to $r^{-p_{\ast }}$, the projection of a spherical power-law emissivity gives $I_{\ast }\propto R^{1-p_{\ast }}$. Illustratively, $p_{\ast }=3$--$4$ gives projected slopes $-2$--$-3$: then $I_d/I_{\ast }\propto R^{p_{\ast }-2}$ rises strongly outward, but annihilation can resemble the steeper stellar template. Gas-related gamma emission must be compared with observed gas columns and the cosmic ray distribution; it is not always proportional to optical light.

A shallow, extended, approximately spherical $R^{-1}$ component is particularly distinctive; an $R^{-3}$ component needs additional evidence. Fit the stellar, gas and halo templates together over a broad radial range, including morphology and emission spectra. A measured slope alone cannot uniquely identify dark matter, since stellar or gas profiles can coincide with a halo slope over a limited interval.

Massive stars rapidly return oxygen, magnesium, neon and other alpha elements through core-collapse supernovae, whereas Type Ia supernovae supply a substantial delayed iron contribution. If the galactic bulge and stellar halo formed much of their stellar mass rapidly, their stars locked in enhanced alpha-element to iron ratios before much delayed iron enrichment, while prolonged star formation in the local galactic disk included more iron-rich gas. The abundance contrast chiefly records different enrichment timescales, with the initial mass function and gas exchange also affecting the ratios.

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.

In the thin-sheet regime, the gas scale height is $h_g(R)={\sigma }_g^2(R)/(2\pi G{\Sigma }_{\ast }(R))$. Hence radius-independent gas velocity dispersion givesFor an exponential galactic disk, ${\Sigma }_{\ast }(R)={\Sigma }_{\ast ,0}{e}^{-R/R_d}$, soThe gas therefore exhibits galactic disk flaring as stellar surface density falls outward. Stellar galactic disks, by contrast, commonly have roughly constant thickness over their main optical region; applying the same scale height relation to the stars requires ${\sigma }_{z,\ast }^2(R)\propto {\Sigma }_{\ast }(R)$, or ${\sigma }_{z,\ast }\propto e^{-R/(2R_d)}$ for an exponential galactic disk. This comparison does not follow by imposing a radially constant stellar velocity dispersion: that extra assumption would predict stellar galactic disk flaring too. The constancy used in the vertical calculation is constancy with height locally, not a mandatory constancy with radius.

Assume a steady, locally plane-symmetric galactic disk, with negligible radial tilt term in the vertical Jeans equation. Write $\rho (z)$ for the stellar mass density, and $K_z={\partial }_z\Phi$ for the upward gradient of the Newtonian gravitational potential; the actual vertical gravitational acceleration is $-K_z$. Constant vertical velocity dispersion givesIntegrating the plane-parallel Poisson equation for Newtonian gravity from the symmetry plane gives $K_z(z)=4\pi G{\int }_0^z\rho (z^{\prime })d{z}^{\prime }$. Outside most of the mass, this tends to $2\pi G\Sigma$ for $z>0$, where the surface density includes both sides of the galactic disk. Treating this field as constant and integrating the vertical Jeans equation yieldsThis exponential is the thin-sheet approximation or the outer tail, not the exact profile at every height. To see the distinction without an additional assumption, set $\rho ={\rho }_0{e}^{-u}$ with $u=(\Phi -\Phi (0))/{\sigma }^2$. The Poisson equation for Newtonian gravity becomes $u^{\prime \prime }=(4\pi G{\rho }_0/{\sigma }^2)e^{-u}$, with $u(0)=u^{\prime }(0)=0$, and is solved bySince $\Sigma =2{\rho }_0{z}_0$, the exact self-gravitating isothermal slab has $z_0={\sigma }^2/(\pi G\Sigma )$. At large $|z|/z_0$, ${\mathrm{sech}}^2(z/z_0)\sim 4e^{-2|z|/z_0}$, recovering precisely the boxed scale height; near the plane its slope vanishes, whereas an exponential has a cusp.

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.

Low-mass main sequence stars live for much of the Milky Way's age, retaining representatives of many generations, while phase mixing disperses their birth associations around their orbits. In an approximately axisymmetric galactic disk, a sufficiently complete local census can therefore sample the mixed annular stellar population at the Sun's radius despite its small volume. This is a representative local-radius sample, not automatically an unbiased whole-disc census: different vertical scale heights, selection effects and radial migration require corrections when inferring integrated population fractions.