Past exam of the formation of Structure in the Universe course of the University of Cambridge / 2015 / Paper 1 (Paper 2/FSU: Formation of Structure in the Universe)

A spiral density wave is an organized density and gravitational-potential perturbation through which stars and gas can move; it is not generally a permanently labelled set of stars. Differential rotation winds a material arm rapidly. In a wave description, a perturbation has phase $m\varphi +\int kdr-\omega t$, with pattern speed ${\Omega }_p=\omega /m$ distinct from the local stellar angular speed.

For a thin fluid disk the local tightly wound dispersion relation isRotation and pressure resist self-gravity. A collisionless stellar disk requires its stellar response function instead of simply replacing $c_s$ by a dispersion; the familiar local axisymmetric stability parameter is $Q_{\star }={\sigma }_R\kappa /(3.36G\Sigma )$. This is distinct from the gas $Q_g=c_s\kappa /(\pi G\Sigma )$.

Waves exchange angular momentum at Lindblad resonances and near corotation radius. Swing amplification can strengthen a leading disturbance as shear turns it into a trailing one. Bars or companions can force a coherent spiral; recurrent self-excited disturbances also occur. Thus a single indefinitely steady wave is one model, not a description established for every spiral galaxy. Gas compression and shocks along arms can organize star formation, while gravitational torques redistribute angular momentum. Arm brightness, recent star formation and stellar mass density need not trace precisely the same pattern.^[ref]

Supermassive black holes are compact galactic nuclei with masses typically of millions to billions of solar masses. Their gravitational influence is measured from the motion of surrounding matter; accretion luminosity is evidence of feeding, not a necessary property of every black hole. Keplerian stellar orbits near the Galactic Centre and molecular-maser rotation curves in other nuclei provide particularly direct compact-mass constraints. A rough black-hole sphere of influence is $r_{\mathrm{infl}}\sim G{M}_{\bullet }/{\sigma }_{\star }^2$; resolving it helps dynamical mass estimates.

Accretion converts binding energy into radiation,The latter is the Eddington luminosity for isotropic radiation opposed by electron-scattering opacity, not an absolute luminosity ceiling for every geometry. The retained black-hole mass grows at $(1-\eta ){\dot{M}}_{\mathrm{in}}$. Sustained feeding can produce an active galactic nucleus, whereas gas-poor nuclei may be faint. Loss of gas angular momentum is a central feeding problem.

Black-hole mass correlates with classical-bulge stellar velocity dispersion through the M-sigma relation. Accretion-powered radiation, winds and jets can heat or remove gas, linking nuclear growth to galaxy evolution. Such correlations do not establish one universal feedback history, and disks and pseudobulges do not follow the same tight relations as classical bulges. Seed formation, accretion and mergers all contribute; seed masses and feeding histories remain model dependent.^[ref]

Galactic star formation requires gas to collect into cold dense regions that can become self-gravitating. Disk-scale gravity, turbulent motions, shear, pressure and magnetic fields affect this collection; formation then proceeds within molecular clouds. A useful local parametrization iswhere the star-formation efficiency per free-fall time summarizes inefficient conversion rather than requiring every cloud to collapse at once.

Across galactic disks the empirical Kennicutt–Schmidt law relates star formation rate surface density to gas surface density, ${\Sigma }_{\mathrm{SFR}}\propto {\Sigma }_g^N$. A slope near $1.4$ is a familiar integrated total-gas fit, while molecular-gas relations in ordinary disks can be closer to linear; the slope and normalization depend on tracer, averaging scale and environment. The gas-depletion time $t_{\mathrm{dep}}=M_g/\Psi$ is the corresponding fuel-consumption scale, not a prediction that inflow and recycling cease.

Massive stars supply ionizing radiation, winds and supernova energy and momentum. This stellar feedback can disperse clouds, sustain turbulence and drive outflows, regulating rather than simply extinguishing all future star formation. Fresh accretion and recycled gas replenish the reservoir. UV and recombination-line emission trace young stars, while infrared emission traces absorbed and re-radiated light; extinction, stellar populations and the assumed initial mass function affect inferred rates. A galactic rate reflects supply, cloud formation and feedback together, rather than only the free-fall time of its densest gas.^[ref]

In a matter-dominated Einstein-de Sitter universe, the growing linear density contrast is proportional to the scale factor. When $|\delta |\sim 1$, linear evolution no longer describes collapse. An overdense region initially expands more slowly than the background, reaches cosmological turnaround, and then contracts under self-gravity.

The spherical top-hat collapse model supplies a useful idealization: pressureless shells obey $\ddot{r}=-GM/r^2$. At formal collapse, the linearly extrapolated growing mode isThis is a linear extrapolation, not the true finite nonlinear density at that instant. In the ideal energy-conserving spherical model, virial equilibrium gives a virial radius about half the turnaround radius; its mean density at the collapse time is $18{\pi }^2\simeq 178$ times the cosmic matter density.

Real perturbations are nonspherical: tides and anisotropic collapse produce sheets, filaments and haloes. Collisionless dark matter develops shell crossing and phase mixing rather than radiatively losing its orbital energy. Baryonic gas can shock, cool and collect within the resulting potential wells. Hierarchical merging then assembles larger systems where small-scale fluctuations have sufficient initial power. Press-Schechter formalism connects a collapse threshold with the statistics of the smoothed initial field, but halo abundances and internal structure need more than the single spherical model. These ideas apply after perturbations become nonlinear; one must not extend pressureless subhorizon equations through all radiation-era or horizon-scale conditions.^[ref]

With constant isothermal sound speed $a_T$, the equation of state is $p=a_T^2\rho$. Spherical hydrostatic equilibrium and mass continuity requireTry $\rho =A/r^2$. With no central point mass, $M(r)=4\pi Ar$, and hydrostatic balance becomes $2a_T^{2}A=4\pi G{A}^2$. The positive-density branch is thereforeThe constant external pressure fixes the truncation radius $r_s$ through $p(r_s)=p_0$, givingThis is the pressure-truncated singular isothermal gas sphere, one hydrostatic solution, not a claim that every isothermal cloud has a singular centre. Its central density diverges, but its enclosed mass tends to zero as $r\to 0$.

For inside-out collapse of a singular isothermal sphere, the expansion front is at $r_f=a_{T}t$. The original mass swept by this front isThus the only time-independent accretion scale is $a_T^3/G$. The front-sweeping estimate has coefficient two, but gas takes time to fall from the front to the object; it is not an exact central accretion coefficient. Writefor fixed mean particle mass $\mu m_H$. Determining the dimensionless coefficient $m_0$ requires the collapse solution; it is not fixed by front propagation alone and is not needed for the requested temperature scaling.^[ref]

Well inside the front, where the point mass dominates and infall is approximately pressureless and locally steady, $|v_r|\simeq \sqrt{2G{M}_{\star }/r}$. Conserving the radial mass flux gives the free-fall density profile around a growing point mass:The inner slope is $-3/2$, versus $-2$ in the undisturbed outer cloud. For constant ${\dot{M}}_{\star }$ and negligible initial object mass, $M_{\star }\propto t$ and the inner normalization varies as $t^{-1/2}$. These are inner asymptotes, not an assertion of an abrupt exact power-law join at the front.

At distance $r$, one grain intercepts power through its projected area and emits through its full surface:The blackbody dust sublimation radius follows by setting $T_d=T_{\mathrm{sub}}$:Here the supplied luminosity is $3.8\times 10^{37}\mathrm{W}$. The grain radius cancels because both cross-sections scale as $a^2$. The temperature exterior to this boundary scales as $T_d\propto r^{-1/2}$. This estimate assumes isotropic illumination, no shielding, and unit absorption/emission efficiency as specified; real grain emissivities can change it. The numerical evaluation uses the external conversion values ${\sigma }_{\mathrm{SB}}\simeq 5.6704\times 10^{-8}\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-4}$ and $1\mathrm{pc}\simeq 3.0857\times 10^{16}\mathrm{m}$.^[ref]

The dust-cooled equilibrium is thermally stable under the stated constant irradiation, with fixed dust properties and dust-to-gas ratio. At a fixed position, absorbed power per grain is independent of its temperature, while emission is $4\pi a^2{\sigma }_{\mathrm{SB}}{T}^4$. A positive temperature perturbation increases cooling, causing a restoring energy loss; a negative perturbation reduces cooling, allowing reheating.

For heat capacity $C>0$ of a grain together with its thermally coupled gas,so the perturbation decays. The thermal stability of externally heated blackbody dust relies on remaining below sublimation with the grain population intact; destruction, variable illumination or finite dust-gas equilibration can change the response. An isobaric parcel with a fixed grain-to-gas ratio has the same restoring sign when heating and cooling are compared per particle.

First specify the perturbation: take fixed gas density, fixed ionic abundance and the same temperature-independent heating, so thermal stability requires $d\mathcal{C}/dT>0$ at an equilibrium. Let $\theta =\Delta E/k_B$, $n_e$ be the collider density and $n_A=n_1+n_2$ the atomic density. A two-level optically thin balance isThe critical density of an atomic transition is approximately $A_{21}/q_{21}$. At low density, each rare collisional excitation is followed by an escaping photon, soThe isochoric thermal stability of a collisionally excited line therefore givesAt $T=2\theta$ the linear restoring coefficient vanishes; the cooling maximum is marginal to first order and is not robust two-sided stability.

At high density, collisional de-excitation competes successfully with radiative decay, and collisions establish local thermodynamic equilibrium level populations. Detailed balance gives $n_2/n_1=g{e}^{-\theta /T}$, where $g=g_2/g_1$. The cooling is thenThus the high-density fixed-density branch is stable for finite $T$, becoming weakly restoring as the upper fraction saturates at high temperature. It is incorrect to retain the low-density factor $T^{-1/2}$ in the radiative cooling law once populations are thermalized.^[ref]

The isobaric thermal stability of a two-level coolant is a distinct question. If the intended heating is a fixed rate per particle, then $n\propto T^{-1}$ and stability tests the cooling per particle. The low-density boundary becomes $T<2\theta /3$; the high-density cooling per particle remains proportional to $f_u$, so remains stable. If instead heating is constant per volume, the isobaric low-density boundary is $T<2\theta /5$, and the high-density stability condition is $\theta (1-f_u)/T>1$. The source gives no pressure constraint or explicit heating normalization: the boxed $2\theta$ result is the fixed-density answer, with these conditional alternatives stated to avoid asserting a universal boundary.

For constant inflow $\Phi$ and efficiency $ϵ>0$, with no recycling or outflow, mass conservation gives the constant-inflow gas regulator solution:The equilibrium $M_{g,\mathrm{eq}}=\Phi /ϵ$ is stable: a gas-mass perturbation obeys $\delta {\dot{M}}_g=-ϵ\delta M_g$. At late times the star formation rate approaches $\Psi =\Phi$, while the gas reservoir approaches a constant.

Integrating ${\dot{M}}_{\star }=ϵM_g$ with $M_{\star }(0)=0$ gives the stellar mass in a constant-inflow gas regulator:It also follows immediately from $M_g+M_{\star }=M_0+\Phi t$. The gas mass equilibrates; the stellar mass generally keeps growing. With $\Phi =0$, these reduce to exponential gas depletion and $M_{\star }=M_0(1-e^{-ϵt})$.

A constant inflow is implicit in the requested steady asymptote, not automatic for an arbitrary prescribed rate. For time-dependent $\Phi (t)$ the exact integrating-factor result isa constant late-time gas mass then requires a corresponding asymptotically steady supply. Both printed mark groups remain in this single unlettered continuation of part(a).

Let $D$ be the separation of galaxy centres and $x$ a small displacement from the satellite centre along the line of centres. Expanding the central point-mass field gives the radial differential accelerationInside a uniform-density satellite, the restoring self-gravity isComparing the two at the same displacement yields the differential-gravity disruption limit of a uniform satellite:The satellite radius cancels because both accelerations are linear in $x$. This is the instantaneous tidal-gravity estimate requested by the given criterion, not a full rotating Roche calculation. Including centrifugal terms for a circular corotating orbit replaces the radial coefficient two by three, giving $D_{\min }\simeq [9M/(4\pi \rho ){]}^{1/3}$ while preserving the required scaling. An eccentric orbit requires checking the stronger tide near periapsis, and extended density profiles require a more detailed restoring-force calculation.

Put $H=2/(3t)$ and let $f_B={\rho }_B/\bar{\rho }$, $f_D=1-f_B$ be constant background matter fractions. Define the total matter perturbation $S=f_B{\Delta }_B+f_D{\Delta }_D$ and the relative mode $D={\Delta }_D-{\Delta }_B$. Adding with the background weights and subtracting the two given equations respectively givesThe total and relative density modes of two pressureless matter fluids are consequentlyThese are the general consequences of the supplied equations; a single ratio additionally needs initial-mode assumptions.

For the usual dark-matter-dominated growing-mode approximation, take $f_B\ll 1$, neglect the decaying modes, and set the relative-velocity mode $E=0$. Then ${\Delta }_D\propto t^{2/3}\propto R$. With ${\Delta }_B(t_r)=0$ at recombination, $C={\Delta }_D(t_r)$, so the baryon catch-up approximation after recombination isAt redshifts $z,z_0\gg 1$, replacing $1+z$ by $z$ gives the printed ${\Delta }_B\simeq {\Delta }_D(1-z/z_0)$. The printed relation is a high-redshift approximation with a specified relative mode, not an exact identity for arbitrary solutions of the two-fluid equations.

For example, retaining a finite baryon fraction with a growing total mode and constant relative mode gives, with $x=R/R(t_r)$,Alternatively, if initially both ${\Delta }_B=0$ and ${\dot{\Delta }}_B=0$ while the dark-matter growing mode is already present, then in the dark-matter-dominated limit the relative transient is not zero andThis shows why the missing initial-velocity convention matters. TODO: the source does not supply the initial relative velocity or a baryon-fraction approximation; those inputs are needed to select the quoted special case uniquely. The general solution and conditional specializations above remain valid.

The physical conclusion is robust: after decoupling, baryonic gas can fall into pre-existing dark-matter potential wells and its fractional density perturbation catches up. The absolute constant offset need not vanish; it becomes small relative to the growing common mode. The pressureless approximation applies on scales above the baryonic Jeans scale and while perturbations remain linear; subsequent collapse requires nonlinear treatment.^[ref]