ABSTRACT
Planck data provide precise constraints on cosmological parameters when assuming the base ΛCDM model, including a 0.17 per cent measurement of the age of the Universe, $t_0=13.797 \pm ...0.023\, {\rm Gyr}$. However, the persistence of the ‘Hubble tension’ calls the base ΛCDM model’s completeness into question and has spurred interest in models such as early dark energy (EDE) that modify the assumed expansion history of the Universe. We investigate the effect of EDE on the redshift–time relation z↔t and find that it differs from the base ΛCDM model by at least ${\approx } 4{{\ \rm per\ cent}}$ at all t and z. As long as EDE remains observationally viable, any inferred t ← z or z ← t quoted to a higher level of precision do not reflect the current status of our understanding of cosmology. This uncertainty has important astrophysical implications: the reionization epoch – 10 > z > 6 – corresponds to disjoint lookback time periods in the base ΛCDM and EDE models, and the EDE value of t0 = 13.25 ± 0.17 Gyr is in tension with published ages of some stars, star clusters, and ultrafaint dwarf galaxies. However, most published stellar ages do not include an uncertainty in accuracy (due to, e.g. uncertain distances and stellar physics) that is estimated to be $\sim 7\!-\!10{{\ \rm per\ cent}}$, potentially reconciling stellar ages with $t_{0,\rm EDE}$. We discuss how the big data era for stars is providing extremely precise ages ($\lt 1{{\ \rm per\ cent}}$) and how improved distances and treatment of stellar physics such as convection could result in ages accurate to $4\!-\!5{{\ \rm per\ cent}}$, comparable to the current accuracy of t↔z. Such precise and accurate stellar ages can provide detailed insight into the high-redshift Universe independent of a cosmological model.
Abstract
Motivated by the stellar fossil record of Local Group (LG) dwarf galaxies, we show that the star-forming ancestors of the faintest ultra-faint dwarf galaxies (UFDs; MV
∼ −2 or M
⋆ ∼ 102 at ...z = 0) had ultraviolet (UV) luminosities of M
UV ∼ −3 to −6 during reionization (z ∼ 6–10). The existence of such faint galaxies has substantial implications for early epochs of galaxy formation and reionization. If the faint-end slopes of the UV luminosity functions (UVLFs) during reionization are steep (α ≲ −2) to M
UV ∼ −3, then (i) the ancestors of UFDs produced >50 per cent of UV flux from galaxies; (ii) galaxies can maintain reionization with escape fractions that are more than two times lower than currently adopted values; (iii) direct Hubble Space Telescope and James Webb Space Telescope observations may detect only ∼10–50 per cent of the UV light from galaxies; and (iv) the cosmic star formation history increases by ≳ 4–6 at z ≳ 6. Significant flux from UFDs, and resultant tensions with LG dwarf galaxy counts, is reduced if the high-redshift UVLF turns over. Independent of the UVLF shape, the existence of a large population of UFDs requires a non-zero luminosity function to M
UV ∼ −3 during reionization.
With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest-mass, ultrafaint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these ...faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their six-dimensional phase-space information, we calculate the orbital histories of 13 ultrafaint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: (i) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); (ii) Magellanic satellites that were recently captured by the MCs < 1 Gyr ago (Reticulum 2, Phoenix 2); (iii) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and (iv) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites.
Abstract
We study the z = 0 gas kinematics, morphology and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning Mstar = 106–11 ...M⊙. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest mass galaxies (Mstar < 108 M⊙), gas fails to form a morphological disc and is primarily dispersion and pressure supported. At intermediate masses (Mstar = 108–10 M⊙), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating discs to irregular spheroids with negligible net rotation. All the high-mass (Mstar = 1010–11 M⊙) galaxies form rotationally supported gas discs. Many of the haloes whose galaxies fail to form discs harbour high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark matter halo increases significantly with galaxy mass, from 〈jgas〉/〈jDM〉 ∼ 0.1 at $M_{\rm star}=10^{6\text{--}7}\, \rm M_{{\odot }}$ to 〈jgas〉/〈jDM〉 ∼ 2 at Mstar = 1010–11 M⊙. The reduced rotational support in the lowest mass galaxies owes to (a) stellar feedback and the UV background suppressing the accretion of high angular momentum gas at late times, and (b) stellar feedback driving large non-circular gas motions. We broadly reproduce the observed scaling relations between galaxy mass, gas rotation velocity, size and angular momentum, but may somewhat underpredict the incidence of disky, high angular momentum galaxies at the lowest observed masses (Mstar = (106–2 × 107) M⊙). Stars form preferentially from low angular momentum gas near the galactic centre and are less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates the simulated galaxies’ stellar angular momentum, particularly at low masses.
Abstract
We measure the ionizing photon production efficiency (
ξ
ion
) of low-mass galaxies (10
7.8
–10
9.8
M
⊙
) at 1.4 <
z
< 2.7 to better understand the contribution of dwarf galaxies to the ...ionizing background and reionization. We target galaxies that are magnified by strong-lensing galaxy clusters and use Keck/MOSFIRE to measure nebular emission-line fluxes and Hubble Space Telescope to measure the rest-UV and rest-optical photometry. We present two methods of stacking. First, we take the average of the log of H
α
-to-UV luminosity ratios (
L
H
α
/
L
UV
) of galaxies to determine the standard log(
ξ
ion
). Second, we take the logarithm of the total
L
H
α
over the total
L
UV
. We prefer the latter, as it provides the total ionizing UV luminosity density of galaxies when multiplied by the nonionizing UV luminosity density. log(
ξ
ion
) calculated from the second method is ∼0.2 dex higher than the first method. We do not find any strong dependence between log(
ξ
ion
) and stellar mass, far-UV magnitude (
M
UV
), or UV spectral slope (
β
). We report a value of log(
ξ
ion
) ∼ 25.47 ± 0.09 for our UV-complete sample (
) and ∼25.37 ± 0.11 for our mass-complete sample (7.8 < log(
M
*
) < 9.8). These values are consistent with measurements of more massive, more luminous galaxies in other high-redshift studies that use the same stacking technique. Our log(
ξ
ion
) is 0.2–0.3 dex higher than low-redshift galaxies of similar mass, indicating an evolution in the stellar properties, possibly due to metallicity or age. We also find a correlation between log(
ξ
ion
) and the equivalent widths of H
α
and O
iii
λ
5007 fluxes, confirming that these equivalent widths can be used to estimate
ξ
ion
.
We use data from Gaia's second data release (DR2) to constrain the initial-final mass relation (IFMR) for field stars with initial masses 0.9 min/M 8. Precise parallaxes have revealed unprecedented ...substructure in the white dwarf (WD) cooling sequence on the color-magnitude diagram (CMD). Some of this substructure stems from the diversity of WD atmospheric compositions, but the CMD remains bimodal even when only spectroscopically confirmed DA WDs are considered. We develop a generative model to predict the CMD for DA WDs as a function of the initial mass function, stellar age distribution, and a flexibly parameterized IFMR. We then fit the CMD of 1100 bright DA WDs within 100 pc, for which atmospheric composition and completeness are well understood. The resulting best-fit IFMR flattens at 3.5 min/M 5.5, producing a secondary peak in the WD mass distribution at mWD ∼ 0.8 M . Our IFMR is broadly consistent with weaker constraints obtained from binaries and star clusters in previous work but represents the clearest observational evidence obtained to date of theoretically predicted nonlinearity in the IFMR. A visibly bimodal CMD is only predicted for mixed-age stellar populations: in single-age clusters, more massive WDs reach the bottom of the cooling sequence before the first lower-mass WDs appear. This may explain why bimodal cooling sequences have thus far evaded detection in cluster CMDs.
Abstract
We develop a data-driven spectral model for identifying and characterizing spatially unresolved multiple-star systems and apply it to APOGEE DR13 spectra of main-sequence stars. Binaries and ...triples are identified as targets whose spectra can be significantly better fit by a superposition of two or three model spectra, drawn from the same isochrone, than any single-star model. From an initial sample of ∼20 000 main-sequence targets, we identify ∼2500 binaries in which both the primary and secondary stars contribute detectably to the spectrum, simultaneously fitting for the velocities and stellar parameters of both components. We additionally identify and fit ∼200 triple systems, as well as ∼700 velocity-variable systems in which the secondary does not contribute detectably to the spectrum. Our model simplifies the process of simultaneously fitting single- or multi-epoch spectra with composite models and does not depend on a velocity offset between the two components of a binary, making it sensitive to traditionally undetectable systems with periods of hundreds or thousands of years. In agreement with conventional expectations, almost all the spectrally identified binaries with measured parallaxes fall above the main sequence in the colour–magnitude diagram. We find excellent agreement between spectrally and dynamically inferred mass ratios for the ∼600 binaries in which a dynamical mass ratio can be measured from multi-epoch radial velocities. We obtain full orbital solutions for 64 systems, including 14 close binaries within hierarchical triples. We make available catalogues of stellar parameters, abundances, mass ratios, and orbital parameters.
Abstract
We use Monte Carlo simulations to explore the statistical challenges of constraining the characteristic mass (m
c
) and width (σ) of a lognormal sub-solar initial mass function (IMF) in ...Local Group dwarf galaxies using direct star counts. For a typical Milky Way (MW) satellite (M
V
= −8), jointly constraining m
c
and σ to a precision of ≲ 20 per cent requires that observations be complete to ≲ 0.2 M⊙, if the IMF is similar to the MW IMF. A similar statistical precision can be obtained if observations are only complete down to 0.4 M⊙, but this requires measurement of nearly 100× more stars, and thus, a significantly more massive satellite (M
V
∼ −12). In the absence of sufficiently deep data to constrain the low-mass turnover, it is common practice to fit a single-sloped power law to the low-mass IMF, or to fit m
c
for a lognormal while holding σ fixed. We show that the former approximation leads to best-fitting power-law slopes that vary with the mass range observed and can largely explain existing claims of low-mass IMF variations in MW satellites, even if satellite galaxies have the same IMF as the MW. In addition, fixing σ during fitting leads to substantially underestimated uncertainties in the recovered value of m
c
(by a factor of ∼4 for typical observations). If the IMFs of nearby dwarf galaxies are lognormal and do vary, observations must reach down to ∼m
c
in order to robustly detect these variations. The high-sensitivity, near-infrared capabilities of the James Webb Space Telescope and Wide-Field Infrared Survey Telescope have the potential to dramatically improve constraints on the low-mass IMF. We present an efficient observational strategy for using these facilities to measure the IMFs of Local Group dwarf galaxies.
Abstract
We present a suite of 15 cosmological zoom-in simulations of isolated dark matter haloes, all with masses of M
halo ≈ 1010 M⊙ at z = 0, in order to understand the relationship among halo ...assembly, galaxy formation and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (fire) project and are performed at extremely high resolution (m
baryon = 500 M⊙, m
dm = 2500 M⊙). The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way (M
*/M⊙ ≈ 105 − 107). This non-negligible spread in stellar mass at z = 0 in haloes within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity V
max, both of which are tightly linked to halo formation time. Much of this dependence of M
* on a second parameter (beyond M
halo) is a direct consequence of the M
halo ∼ 1010 M⊙ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming haloes remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Haloes that form galaxies with M
⋆ ≳ 2 × 106 M⊙ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius r
1/2. Lower-mass galaxies do not modify their host dark matter haloes at the mass scale studied here. This apparent stellar mass threshold of M
⋆ ≈ 2 × 106 − 2 × 10− 4 M
halo is broadly consistent with previous work and provides a testable prediction of fire feedback models in Λcold dark matter.