The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

We explore the structure of the element abundance--age--orbit distribution of the stars in the Milky Way's low-$\alpha$ disk, by (re-)deriving precise [Fe/H], [X/Fe] and ages, along with orbits, for red clump stars from the APOGEE survey. There has been a long-standing theoretical expectation and observational evidence that metallicity ([Fe/H]) and age are informative about a star's orbit, e.g. about its angular momentum and the corresponding mean Galactocentric distance or its vertical motion. Indeed, our analysis of the APOGEE data confirms that [Fe/H] or age alone can predict the stars' orbits far less well than the combination of the two. Remarkably, we find and show explicitly, that for known [Fe/H] and age, the other abundances [X/Fe] of Galactic disk stars can be predicted well (on average to 0.02 dex) across a wide range of Galactocentric radii, and therefore provide little additional information, e.g. for predicting their orbit. While the age-abundance space for metal poor stars and potentially for stars near the Galactic center is rich or complex, for the bulk of the Galaxy's low-$\alpha$ disk it is simple: [Fe/H] and age contain most information, unless [X/Fe] can be measured to 0.02, or better. Consequently, we do not have the precision with current (and likely near-future) data to assign stars to their individual (coeval) birth clusters, from which the disk is presumably formed. We can, however, place strong constraints on future models of galactic evolution, chemical enrichment and mixing.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

The large-scale structure of the Universe should soon be measured at high redshift during the Epoch of Reionization (EoR) through line-intensity mapping. A number of ongoing and planned surveys are using the 21 cm line to trace neutral hydrogen fluctuations in the intergalactic medium (IGM) during the EoR. These may be fruitfully combined with separate efforts to measure large-scale emission fluctuations from galactic lines such as [CII], CO, H-$\alpha$, and Ly-$\alpha$ during the same epoch. The large scale power spectrum of each line encodes important information about reionization, with the 21 cm power spectrum providing a relatively direct tracer of the ionization history. Here we show that the large scale 21 cm power spectrum can be extracted using only cross-power spectra between the 21 cm fluctuations and each of two separate line-intensity mapping data cubes. This technique is more robust to residual foregrounds than the usual 21 cm auto-power spectrum measurements and so can help in verifying auto-spectrum detections. We characterize the accuracy of this method using numerical simulations and find that the large-scale 21 cm power spectrum can be inferred to an accuracy of within 5% for most of the EoR, reaching 0.6% accuracy on a scale of $k\sim0.1\,\text{Mpc}^{-1}$ at $\left< x_i \right> = 0.36$ ($z = 8.34$ in our model). An extension from two to $N$ additional lines would provide $N(N-1)/2$ cross-checks on the large-scale 21 cm power spectrum. This work strongly motivates redundant line-intensity mapping surveys probing the same cosmological volumes.

The orbital properties of stars in the disk are signatures of their formation, but they are also expected to change over time due to the dynamical evolution of the Galaxy. Stellar orbits can be quantified by three dynamical actions, J_r, L_z, and J_z, which provide measures of the orbital eccentricity, guiding radius, and non-planarity, respectively. Changes in these dynamical actions over time reflect the strength and efficiency of the evolutionary processes that drive stellar redistributions. We examine how dynamical actions of stars are correlated with their age using two samples of stars with well-determined ages: 78 solar twin stars (with ages to ~5%) and 4376 stars from the APOKASC2 sample (~20%). We compute actions using spectroscopic radial velocities from previous surveys and parallax and proper motion measurements from Gaia DR2. We find weak gradients in all actions with stellar age, of (7.51 +/- 0.52, -29.0 +/- 1.83, 1.54 +/- 0.18) kpc km/s/Gyr for J_r, L_z, and J_z, respectively. There is, however, significant scatter in the action-age relation. We caution that our results will be affected by the restricted spatial extent of our sample, particularly in the case of J_z. Nevertheless, these action-age gradients and their associated variances provide strong constraints on the efficiency of the mechanisms that drive the redistribution of stellar orbits over time and demonstrate that actions are informative as to stellar age. The shallow action-age gradients combined with the large dispersion in each action at a given age, however, renders the prospect of age inference from orbits of individual stars bleak. Using the precision measurements of [Fe/H] and [$\alpha$/Fe] for our stars we investigate the abundance-action relationship and find weak correlations. Similar to our stellar age results, dynamical actions afford little discriminating power between low- and high-$\alpha$ stars.

The amplitude of redshifted 21 cm fluctuations during the Epoch of Reionization (EoR) is expected to show a distinctive "rise and fall" behavior with decreasing redshift as reionization proceeds. On large scales (k <~ 0.1 Mpc^{-1}) this can mostly be characterized by evolution in the product of the mean 21 cm brightness temperature and a bias factor, b_21(z). This quantity evolves in a distinctive way that can help in determining the average ionization history of the intergalactic medium (IGM) from upcoming 21 cm fluctuation data sets. Here we consider extracting <T_21> b_21(z) using a combination of future redshifted 21 cm and [CII] line-intensity mapping data sets. Our method exploits the dependence of the 21 cm-[CII]-[CII] cross-bispectrum on the shape of triangle configurations in Fourier space. This allows one to determine <T_21> b_21(z) yet, importantly, is less sensitive to foreground contamination than the 21 cm auto-spectrum, and so can provide a valuable cross-check. We compare the results of simulated bispectra with second-order perturbation theory: on large scales the perturbative estimate of <T_21> b_21(z) matches the true value to within 10% for <x_i> <~ 0.8. We consider the 21 cm auto-bispectrum and show that this statistic may also be used to extract the 21 cm bias factor. Finally, we discuss the survey requirements for measuring the cross-bispectrum. Although we focus on the 21 cm-[CII]-[CII] bispectrum during reionization, our method may be of broader interest and can be applied to any two fields throughout cosmic history.