A dense 0.1-solar-mass star in a 51-minute-orbital-period eclipsing binary.

Of more than a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes1-9. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution before interacting with the white dwarf10-14, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs13-18. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute-orbital-period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater owing to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modelling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs.

Where are the most ancient stars in the Milky Way?

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disc galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars (z form ≳ 5) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at z = 0 than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed in situ but are accreted during hierarchical assembly. These ex situ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed in situ in the solar neighbourhood, and in some simulations, even in the galactic centre. Secondly, old stars formed in situ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at z ≳ 2, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines away from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before z = 5. We predict most of these targets to have higher metallicity (-3 < [Fe/H] < -2) than the most extreme metal-poor stars.

The origin of the diverse morphologies and kinematics of Milky Way-mass galaxies in the FIRE-2 simulations.

We use hydrodynamic cosmological zoom-in simulations from the Feedback in Realistic Environments project to explore the morphologies and kinematics of 15 Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90 per cent of their stellar mass within 2.5 kpc to well-ordered discs that reach ≳15 kpc. The gas in our galaxies always forms a thin, rotation-supported disc at z = 0, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disc-like orbits and with the radial extent of the stellar disc. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disc after z ~ 1 develop well-ordered stellar discs. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-z mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact discs in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in z = 0 MW-mass galaxies (even z = 0 bulge stars) form in a disc: ≳60-90 per cent of stars begin their lives rotationally supported.

Gas kinematics in FIRE simulated galaxies compared to spatially unresolved HI observations.

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotationally supported gas discs produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing V/σ and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of ~2000 low-redshift, HI-detected galaxies with Mstar = 107-11 M☉ and compare to the simulated galaxies. At Mstar ≳ 1010 M☉, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotationally supported discs. Both the observed and simulated line profiles become more Gaussian like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotational to dispersion support more strongly: at Mstar 108-10 M, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion supported. Most of the lower-mass-simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotational versus dispersion support.