Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Misaligned low mass galaxies typically have boosted BH luminosity, BH growth and have had more energy injected into the gas over the last 8 Gyr in comparison to aligned galaxies. These properties likely lead to outflows and gas loss, in agreement with active low mass galaxies in observations. We show that splitting on BH luminosity at $z=0$ produces statistically consistent distributions of kinematic misalignment at $z=0$, however, splitting on the maximum BH luminosity over the last 8 Gyrs does not. While instantaneous correlation at $z=0$ is difficult due to misalignment persisting on longer timescales, the relationship between BH activity and misalignment is clear. High mass quenched galaxies with misalignment typically have similar BH luminosities, show no overall gas loss, and have typically lower gas phase metallicity over the last 8 Gyrs in comparison to those aligned; suggesting external origin.

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_{\phi}\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.

The morphology of thin stellar streams can be used to test the nature of dark matter. It is therefore crucial to extend searches for globular cluster streams to other galaxies than the Milky Way. In this paper, we investigate the current and future prospects of detecting globular cluster streams in external galaxies in resolved stars (e.g. with WFIRST) and using integrated light (e.g. with HSC, LSST and Euclid). In particular, we inject mock-streams to data from the PAndAS M31 survey, and produce simulated M31 backgrounds mimicking what WFIRST will observe in M31. Additionally, we estimate the distance limit to which globular cluster streams will be observable. Our results demonstrate that for a 1 hour (1000 sec.) exposure, using conservative estimates, WFIRST should detect globular cluster streams in resolved stars in galaxies out to distances of ~3.5 Mpc (~2 Mpc). This volume contains 199 (122) galaxies of which >90% are dwarfs. With integrated light, thin streams can be resolved out to ~100 Mpc with HSC and LSST and to ~600 Mpc with WFIRST and Euclid. The low surface brightness of the streams (typically >30 mag/arcsec$^2$), however, will make them difficult to detect, unless the streams originate from very young clusters. We emphasize that if the external galaxies do not host spiral arms or galactic bars, gaps in their stellar streams provide an ideal test case for evidence of interactions with dark matter subhalos. Furthermore, obtaining a large samples of thin stellar streams can help constrain the orbital structure and hence the potentials of external halos.

In spatially resolved galaxy observations, star formation rate radial profiles are found to correlate with total specific star formation rates. A central depletion in star formation is thought to correlate with the globally depressed star formation rates of, for example, galaxies within the Green Valley. We present, for the first time, radial specific star formation rate profiles for a statistical sample of simulated galaxies from the Illustris and EAGLE large cosmological simulations. For galaxies on the star-forming sequence, simulated specific star formation rate profiles are in loose agreement with observations, although galaxies from the EAGLE simulation are too centrally peaked, and galaxies from Illustris have a steeper decline at large radii. However, both galaxy samples show centrally concentrated star formation for galaxies in the Green Valley at all galaxy stellar masses, indicating that quenching occurs from the outside-in, in strong conflict with observations of inside-out quenching. These results appear in spite of the different feedback models in these two simulations. We conclude that the distribution of star formation within galaxies is a strong additional constraint for star formation and feedback models in simulations, in particular related to the quenching of star formation.

Stars in galaxies form from the cold rotationally supported gaseous disks that settle at the center of dark matter halos. In the simplest models, such angular momentum is acquired early on at the time of collapse of the halo and preserved thereafter, implying a well-aligned spin for the stellar and gaseous component. Observations however have shown the presence of gaseous disks in counterrotation with the stars. We use the Illustris numerical simulations to study the origin of such counterrotation in low mass galaxies ($M_\star = 2 \times 10^9$ - $5 \times 10^{10}\; \rm M_\odot$), a sample where mergers have not played a significant role. Only ${\sim}1\%$ of our sample shows a counterrotating gaseous disk at $z=0$. These counterrotating disks arise in galaxies that have had a significant episode of gas removal followed by the acquisition of new gas with misaligned angular momentum. In our simulations, we identify two main channels responsible for the gas loss: a strong feedback burst and gas stripping during a fly-by passage through a more massive group environment. Once settled, counterrotation can be long-lived with several galaxies in our sample displaying misaligned components consistently for more than $2$ Gyr. As a result, no major correlation with the present day environment or structural properties might remain, except for a slight preference for early type morphologies and a lower than average gas content at a given stellar mass.

We report Giant Metrewave Radio Telescope (GMRT) , Very Large Telescope (VLT) and Spitzer Space Telescope observations of ESO 184$-$G82, the host galaxy of GRB 980425/SN 1998bw, that yield evidence of a companion dwarf galaxy at a projected distance of 13 kpc. The companion, hereafter GALJ193510-524947, is a gas-rich, star-forming galaxy with a star formation rate of $\rm0.004\,M_{\odot}\, yr^{-1}$, a gas mass of $10^{7.1\pm0.1} M_{\odot}$, and a stellar mass of $10^{7.0\pm0.3} M_{\odot}$. The interaction between ESO 184$-$G82 and GALJ193510-524947 is evident from the extended gaseous structure between the two galaxies in the GMRT HI 21 cm map. We find a ring of high column density HI gas, passing through the actively star forming regions of ESO 184$-$G82 and the GRB location. This ring lends support to the picture in which ESO 184$-$G82 is interacting with GALJ193510-524947. The massive stars in GALJ193510-524947 have similar ages to those in star-forming regions in ESO 184$-$G82, also suggesting that the interaction may have triggered star formation in both galaxies. The gas and star formation properties of ESO 184$-$G82 favour a head-on collision with GALJ193510-524947 rather than a classical tidal interaction. We perform state-of-the art simulations of dwarf--dwarf mergers and confirm that the observed properties of ESO 184$-$G82 can be reproduced by collision with a small companion galaxy. This is a very clear case of interaction in a gamma ray burst host galaxy, and of interaction-driven star formation giving rise to a gamma ray burst in a dense environment.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of star formation without apparent cause. We aim to characterize the signatures of minor mergers of dwarf galaxies with dark satellites to aid their observational identification. We explore and quantify a variety of structural, morphological, and kinematic indicators of merging dwarf galaxies and their remnants using a suite of hydrodynamical simulations. The most sensitive indicators of mergers with dark satellites are large asymmetries in the gaseous and stellar distributions, enhanced central surface brightness and starbursts, and velocity offsets and misalignments between the cold gas and stellar components. In general, merging systems span a wide range of values of the most commonly used indicators, while isolated objects tend to have more confined values. Interestingly, we find in our simulations that a significantly off-centered burst of star formation can pinpoint the location of the dark satellite. Observational systems with such characteristics are perhaps the most promising for unveiling the presence of the hitherto, missing satellites.

It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. In this paper we follow the evolution of gas-rich disky dwarf galaxies as they experience a minor merger with a dark satellite. We aim to characterize the effects of such an interaction on the dwarf's star formation, morphology, and kinematical properties. We performed a suite of carefully set-up hydrodynamical simulations of dwarf galaxies that include dark matter, gas, and stars merging with a satellite consisting solely of dark matter. For the host system we vary the gas fraction, disk size and thickness, halo mass, and concentration, while we explore different masses, concentrations, and orbits for the satellite. We find that the interactions cause strong starbursts of both short and long duration in the dwarfs. Their star formation rates increase by factors of a few to 10 or more. They are strongest for systems with extended gas disks and high gas fractions merging with a high-concentration satellite on a planar, radial orbit. In contrast to analogous simulations of Milky Way-mass galaxies, many of the systems experience strong morphological changes and become spheroidal even in the presence of significant amounts of gas. The simulated systems compare remarkably well with the observational properties of a large selection of irregular dwarf galaxies and blue compact dwarfs. This implies that mergers with dark satellites might well be happening but not be fully evident, and may thus play a role in the diversity of the dwarf galaxy population.

In the context of the current $\Lambda$CDM cosmological model small dark matter haloes are abundant and satellites of dwarf galaxies are expected to be predominantly dark. Since low mass galaxies have smaller baryon fractions interactions with these satellites may leave particularly dramatic imprints. We uncover the influence of the most massive of these dark satellites on disky dwarf galaxies and the possible dynamical and morphological transformations that result from these interactions. We use a suite of carefully set-up, controlled simulations of isolated dwarf galaxies. The primary dwarf galaxies have solely a stellar disk in the dark matter halo and the secundaries are completely devoid of baryons. We vary the disk mass, halo concentration, initial disk thickness and inclination of the satellite orbit. The disky dwarf galaxies are heated and disrupted due to the minor merger event, more extremely for higher satellite over disk mass ratios, and the morphology and kinematics are significantly altered. Moreover, for less concentrated haloes the minor merger can completely destroy the disk leaving a low-luminosity spheroidal-like galaxy instead. We conclude that dwarf galaxies are very much susceptible to being disturbed by dark galaxies and that even a minor merger event can significantly disrupt and alter the structure and kinematics of a dwarf galaxy. This process may be seen as a new channel for the formation of dwarf spheroidal galaxies.

One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.