We analyze the kinematics as a function of stellar age for Andromeda (M31) mass analogs from the Illustris and Illustris TNG cosmological simulations. We divide the star particles into four age groups: < 1 Gyr, 1-5 Gyr, 5-10 Gyr, and >10 Gyr, and we compare the kinematics of these groups to that of the neutral gas cells. We calculate rotation curves for the stellar and gaseous components of each analog from 2 kpc to 20 kpc from the center of mass. We find that the lag, or asymmetric drift (AD), between the gas rotation curve and the stellar rotation curve on average increases with stellar age. This finding is consistent with observational measurements of AD in the disk of the Andromeda galaxy. When the M31 analogs are separated into groups based on merger history, we find that there is a difference in the AD of the analogs in the subgroup that corresponds to having had a 4:1 merger in a given time range compared to analogs that have not experienced a 4:1 merger in the same time frame. These differences are the most significant for recent 4:1 merger events (< 4 Gyr ago) and are therefore most visible in the 1-5 Gyr stellar age group where we find a difference in AD as large as 55 km/s. This subset of analogs that have had a 4:1 merger within the last 4 Gyr are also very consistent with AD measurements of stars in M31's disk, providing evidence that M31 may in fact have recently merged with a galaxy nearly 1/4 of its mass. Further work using high resolution zoom-in simulations is required to explore the contribution of internal heating to AD.

We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Str\"omgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

We present image-based evolution of galaxy mergers from the Illustris cosmological simulation at 12 time-steps over 0.5 < z < 5. To do so, we created approximately one million synthetic deep Hubble Space Telescope and James Webb Space Telescope images and measured common morphological indicators. Using the merger tree, we assess methods to observationally select mergers with stellar mass ratios as low as 10:1 completing within +/- 250 Myr of the mock observation. We confirm that common one- or two-dimensional statistics select mergers so defined with low purity and completeness, leading to high statistical errors. As an alternative, we train redshift-dependent random forests (RFs) based on 5-10 inputs. Cross-validation shows the RFs yield superior, yet still imperfect, measurements of the late-stage merger fraction, and they select more mergers in bulge-dominated galaxies. When applied to CANDELS morphology catalogs, the RFs estimate a merger rate increasing to at least z = 3, albeit two times higher than expected by theory. This suggests possible mismatches in the feedback-determined morphologies, but affirms the basic understanding of galaxy merger evolution. The RFs achieve completeness of roughly 70% at 0.5 < z < 3, and purity increasing from 10% at z = 0.5 to 60% at z = 3. At earlier times, the training sets are insufficient, motivating larger simulations and smaller time sampling. By blending large surveys and large simulations, such machine learning techniques offer a promising opportunity to teach us the strengths and weaknesses of inferences about galaxy evolution.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

We present late-time ($\sim$240-260 days after peak brightness) optical photometry and nebular (+236 and +264 days) spectroscopy of SN 2018oh, the brightest Type Ia supernova (SN Ia) observed by the Kepler telescope. The Kepler/K2 30-minute cadence observations started days before explosion and continued past peak brightness. For several days after explosion, SN 2018oh had blue "excess" flux in addition to a normal SN rise. The flux excess can be explained by the interaction between the SN and a Roche-lobe filling non-degenerate companion star. Such a scenario should also strip material from the companion star, that would emit once the SN ejecta become optically thin, imprinting relatively narrow emission features in its nebular spectrum. We search our nebular spectra for signs of this interaction, including close examination of wavelengths of hydrogen and helium transitions, finding no significant narrow emission. We place upper limits on the luminosity of these features of $2.6,\ 2.9\ \mathrm{and}\ 2.1\times10^{37}\ \mathrm{erg\ s^{-1}}$ for H$\alpha$, He I $\lambda$5875, and He I $\lambda$6678, respectively. Assuming a simple model for the amount of swept-up material, we estimate upper mass limits for hydrogen of $5.4\times10^{-4}\ \mathrm{M_{\odot}}$ and helium of $4.7\times10^{-4}\ \mathrm{M_{\odot}}$. Such stringent limits are unexpected for the companion-interaction scenario consistent with the early data. No known model can explain the excess flux, its blue color, and the lack of late-time narrow emission features.

We analyze the kinematics of Andromeda's disk as a function of stellar age by using photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and spectroscopy from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We use HI 21-cm and CO ($\rm J=1 \rightarrow 0$) data to examine the difference between the deprojected rotation velocity of the gas and that of the stars. We divide the stars into four stellar age bins, from shortest lived to longest lived: massive main sequence stars (0.03 Gyr), more luminous intermediate mass asymptotic giant branch (AGB) stars (0.4 Gyr), less luminous intermediate mass AGB stars (2 Gyr), and low mass red giant branch stars (4 Gyr). There is a clear correlation between the offset of the stellar and the gas rotation velocity, or the asymmetric drift: the longer lived populations lag farther behind the gas than short lived populations. We also examine possible causes of the substructure in the rotation curves and find that the most significant cause of scatter in the rotation curves comes from the tilted ring model being an imperfect way to account for the multiple warps in Andromeda's disk.