Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a wide range of galaxy properties, from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a well-ordered disk is present, we calculate the rotational profile <v_phi>(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM_enc/r) assuming spherical symmetry. In the simulated massive high-redshift galaxies, the gas rotation traces the circular velocity reasonably well at intermediate radii r~1-3 kpc, but the two quantities diverge significantly outside that range. At larger radii, gradients in the turbulent pressure can bias dynamical mass measurements low by ~10-40%. In the interior, the assumption of a spherically-symmetric gravitational potential becomes increasingly poor owing to a massive disk component, reducing the gas rotational velocities by >~10%. Finally, in the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

We use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses ($M_{\ast} = 10^9-10^{12} \, {\rm M_{\odot}}$), galaxy types, environments, and assembly histories. We find that the `two-phase' picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation -- namely, the stellar mass growth of galaxies below a few times $10^{11} \, {\rm M_{\odot}}$ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at $z=0$ contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for $M_{\ast} \approx 10^{12} \, {\rm M_{\odot}}$ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, $\sim$50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio $\mu > 1/4$), $\sim$20 per cent from minor mergers ($1/10 < \mu < 1/4$), $\sim$20 per cent from very minor mergers ($\mu < 1/10$), and $\sim$10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.

Galaxy populations at different cosmic epochs are often linked together by comoving cumulative number density in observational studies. Many theoretical works, however, have shown that the number densities of tracked galaxy populations evolve in bulk and spread out over time. We present a number density method for linking progenitor and descendant galaxy populations which takes both of these effects into account. We define probability distribution functions that capture the evolution and dispersion of galaxy populations in comoving number density space, and use these functions to assign galaxies at one redshift $z_f$ probabilities of being progenitors or descendants of a galaxy population at another redshift $z_0$. These probabilities are then used as weights for calculating distributions of physical properties such as stellar mass, star formation rate, or velocity dispersion within the progenitor/descendant population. We demonstrate that this probabilistic method provides more accurate predictions for the evolution of physical properties then either the assumption of a constant number density or the assumption of an evolving number density in a bin of fixed width by comparing the predictions against galaxy populations directly tracked through a cosmological simulation. We find that the constant number density method performs most poorly at recovering galaxy properties, the evolving number method density slightly better, and the probabilistic number density method best of all. The improvement is present for predictions of both stellar mass as well as inferred quantities such as star formation rate and velocity dispersion which were not included in the number density fits. We demonstrate that this method can also be applied robustly and easily to observational data, and provide a code package for doing so.

Galaxy comoving number density is commonly used to forge progenitor/descendant links between observed galaxy populations at different epochs. However, this method breaks down in the presence of galaxy mergers, or when galaxies experience stochastic growth rates. We present a simple analytic framework to treat the physical processes that drive the evolution and diffusion of galaxies within comoving number density space. The evolution in mass rank order of a galaxy population with time is influenced by the galaxy coagulation rate and galaxy "mass rank scatter" rate. We quantify the relative contribution of these two effects to the mass rank order evolution. We show that galaxy coagulation is dominant at lower redshifts and stellar masses, while scattered growth rates dominate the mass rank evolution at higher redshifts and stellar masses. For a galaxy population at $10^{10} M_\odot$, coagulation has been the dominant effect since $z=2.2$, but a galaxy population at $10^{11} M_\odot$ was dominated by mass rank scatter until $z=0.6$. We show that although the forward and backward median number density evolution tracks are asymmetric, the backward median number density evolution can be obtained by convolving the descendant distribution function with progenitor relative abundances. We tabulate fits for the median number density evolution and scatter which can be applied to improve the way galaxy populations are linked in multi-epoch observational datasets.

Massive quiescent galaxies have much smaller physical sizes at high redshift than today. The strong evolution of galaxy size may be caused by progenitor bias, major and minor mergers, adiabatic expansion, and/or renewed star formation, but it is difficult to test these theories observationally. Herein, we select a sample of 35 massive, compact galaxies ($M_* = 1-3 \times 10^{11}$ M$_\odot$, $M_*/R^{1.5} > 10^{10.5}$ M$_\odot$/kpc$^{1.5}$) at $z=2$ in the cosmological hydrodynamical simulation Illustris and trace them forward to $z=0$ to uncover their evolution and identify their descendants. By $z=0$, the original factor of 3 difference in stellar mass spreads to a factor of 20. The dark matter halo masses similarly spread from a factor of 5 to 40. The galaxies' evolutionary paths are diverse: about half acquire an ex-situ envelope and are the core of a more massive descendant, a third survive undisturbed and gain very little mass, 15% are consumed in a merger with a more massive galaxy, and a small remainder are thoroughly mixed by major mergers. The galaxies grow in size as well as mass, and only $\sim$10% remain compact by $z=0$. The majority of the size growth is driven by the acquisition of ex-situ mass. The most massive galaxies at $z=0$ are the most likely to have compact progenitors, but this trend possesses significant dispersion which precludes a direct linkage to compact galaxies at $z=2$. The compact galaxies' merger rates are influenced by their $z=2$ environments, so that isolated or satellite compact galaxies (which are protected from mergers) are the most likely to survive to the present day.

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

We present extensive radio observations of the nearby Type Ibc supernovae 2004cc, 2004dk, and 2004gq spanning 8-1900 days after explosion. Using a dynamical model developed for synchrotron emission from a slightly decelerated shockwave, we estimate the velocity and energy of the fastest ejecta and the density profile of the circumstellar medium. The shockwaves of all three supernovae are characterized by non-relativistic velocities of v ~ (0.1-25)c and associated energies of E ~ (2-10) * 1e47 erg, in line with the expectations for a typical homologous explosion. Smooth circumstellar density profiles are indicated by the early radio data and we estimate the progenitor mass loss rates to be ~ (0.6-13) * 1e-5 M_sun/yr (wind velocity 10^3 km/s). These estimates approach the saturation limit (~1e-4 M_sun/yr) for line-driven winds from Wolf-Rayet stars, the favored progenitors of SNe Ibc including those associated with long-duration GRBs. Intriguingly, at later epochs all three supernovae show evidence for abrupt radio variability that we attribute to large density modulations (factor of ~3-6) at circumstellar radii of r ~ (1-50) * 1e16 cm. If due to variable mass loss, these modulations are associated with progenitor activity on a timescale of ~ 10-100 years before explosion. We consider these results in the context of variable mass loss mechanisms including wind clumping, metallicity-independent continuum-driven ejections, and binary-induced modulations. It may also be possible that the SN shockwaves are dynamically interacting with wind termination shocks, however, this requires the environment to be highly pressurized and/or the progenitor to be rapidly rotating prior to explosion. The proximity of the density modulations to the explosion sites may suggest a synchronization between unusual progenitor mass loss and the SN explosion, reminiscent of Type IIn supernovae. [ABRIDGED]

The origin of highly-magnetized white dwarfs has remained a mystery since their initial discovery. Recent observations indicate that the formation of high-field magnetic white dwarfs is intimately related to strong binary interactions during post-main-sequence phases of stellar evolution. If a low-mass companion, such as a planet, brown dwarf, or low-mass star is engulfed by a post-main-sequence giant, the hydrodynamic drag in the envelope of the giant leads to a reduction of the companion's orbit. Sufficiently low-mass companions in-spiral until they are shredded by the strong gravitational tides near the white dwarf core. Subsequent formation of a super-Eddington accretion disk from the disrupted companion inside a common envelope can dramatically amplify magnetic fields via a dynamo. Here, we show that these disk-generated fields are sufficiently strong to explain the observed range of magnetic field strengths for isolated, high-field magnetic white dwarfs. A higher-mass binary analogue may also contribute to the origin of magnetar fields.