We present population-orbit superposition models for external galaxies based on Schwarzschild's orbit-superposition method, by tagging the orbits with age and metallicity. The models fit the density distributions, as well as kinematic, age and metallicity maps from Integral Field Unit (IFU) spectroscopy observations. We validate the method and demonstrate its power by applying it to mock data, similar to those obtained by the Multi-Unit Spectroscopic Explorer (MUSE) IFU on the Very Large Telescope (VLT). These mock data are created from Auriga galaxy simulations, viewed at three different inclination angles ($\vartheta=40^o, 60^o, 80^o$). Constrained by MUSE-like mock data, our model can recover the galaxy's stellar orbit distribution projected in orbital circularity $\lambda_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $\lambda_z$ and stellar age $t$. A physically motivated age-metallicity relation improves recovering the intrinsic stellar population distributions. We decompose galaxies into cold, warm and hot + counter-rotating components based on their orbit circularity distribution, and find that the surface density, mean velocity, velocity dispersion, age and metallicity maps of each component from our models well reproduce those from simulation, especially for projections close to edge-on. These galaxies exhibit strong global age vs. $\sigma_z$ relation, which is well recovered by our model. The method has the power to reveal the detailed build-up of stellar structures in galaxies, and offers a complement to local resolved, and high-redshift studies of galaxy evolution.

We present a combination of the Schwarzschild orbit-superposition dynamical modelling technique with the spatially-resolved mean stellar age and metallicity maps to uncover the formation history of galaxies. We apply this new approach to a remarkable 5-pointing mosaic of VLT/MUSE observations obtained by Gu\'erou et al. (2016) extending to a maximum galactocentric distance of $120''$ ($\sim 5.6\ {\rm kpc}$) along the major axis, corresponding to $\sim 2.5\ {\rm R}_e$. Our method first identifies `families' of orbits from the dynamical model that represent dynamically-distinct structures of the galaxy. Individual ages and metallicities of these components are then fit for using the stellar-population information. Our results highlight components of the galaxy that are distinct in the combined stellar dynamics/populations space, which implies distinct formation paths. We find evidence for a dynamically-cold, metal-rich disk, consistent with a gradual in-situ formation. This disk is embedded in a generally-old population of stars, with kinematics ranging from dispersion-dominated in the centre to an old, diffuse, metal-poor stellar halo at the extremities. We find also a direct correlation between the dominant dynamical support of these components, and their associated age, akin to the relation observed in the Milky Way. This approach not only provides a powerful model for inferring the formation history of external galaxies, but also paves the way to a complete population-dynamical model.

We apply the Jeans Anisotropic MGE (JAM) dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range $10^{10}<M_*/{\rm M}_{\odot}<10^{11.6}$ that cover a large radial range of $0.1-4.0$ effective radii. We combine SLUGGS and ATLAS$^{\rm 3D}$ datasets to model the total-mass profiles of a sample of 21 fast-rotator galaxies, utilising a hyperparameter method to combine the two independent datasets. The total-mass density profile slope values derived for these galaxies are consistent with those measured in the inner regions of galaxies by other studies. Furthermore, the total-mass density slopes ($\gamma_{\rm tot}$) appear to be universal over this broad stellar mass range, with an average value of $\gamma_{\rm tot}=-2.12\,\pm\,0.05$, i.e. slightly steeper than isothermal. We compare our results to model galaxies from the Magneticum and EAGLE cosmological hydrodynamic simulations, in order to probe the mechanisms that are responsible for varying total-mass density profile slopes. The simulated-galaxy slopes are shallower than the observed values by $\sim0.1-0.3$, indicating that the physical processes shaping the mass distributions of galaxies in cosmological simulations are still incomplete. For galaxies with $M_*>10^{10.7}{\rm M}_{\odot}$ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total mass density slopes for low-mass galaxies are less likely to be determined by merger activity.

We study trends in the slope of the total mass profiles and dark matter fractions within the central half-light radius of 258 early-type galaxies, using data from the volume-limited ATLAS$^{\mathrm{3D}}$ survey. We use three distinct sets of dynamical models, which vary in their assumptions and also allow for spatial variations in the stellar mass-to-light ratio, to test the robustness of our results. We confirm that the slopes of the total mass profiles are approximately isothermal, and investigate how the total-mass slope depends on various galactic properties. The most statistically-significant correlations we find are a function of either surface density, \(\Sigma_e\), or velocity dispersion, \(\sigma_e\). However there is evidence for a break in the latter relation, with a nearly universal logarithmic slope above \(\log_{10}[\sigma_e/(\si{km~s^{-1}})]\sim 2.1\) and a steeper trend below this value. For the 142 galaxies above that critical \(\sigma_e\) value, the total mass-density logarithmic slopes have a mean value \(\left\langle\gamma^\prime\right\rangle = -2.192 \pm 0.016\) (\(1\sigma\) error) with an observed rms scatter of only \(\sigma_{\gamma^\prime}=0.167 \pm 0.016\). Considering the observational errors, we estimate an intrinsic scatter of \(\sigma_{\gamma^\prime}^\mathrm{intr} \approx 0.15\). These values are broadly consistent with those found by strong lensing studies at similar radii and agree, within the tight errors, with values recently found at much larger radii via stellar dynamics or HI rotation curves (using significantly smaller samples than this work).

The Dark Energy Survey (DES) is currently undertaking an observational program imaging $1/4$ of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the DES will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts (GRBs) over five years. Once GRBs are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of GRB activity, collates useful information from archival DES data, and disseminates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that non-public DES data provide for relative photometry of the optical counterparts of GRBs, as well as for identifying key characteristics (e.g., photometric redshifts) of potential GRB host galaxies. We provide the functional details of the DESAlert software as it presently operates, as well as the data products that it produces, and we show sample results from the application of DESAlert to several previously-detected GRBs.

Integral field unit spectrographs allow the 2D exploration of the kinematics and stellar populations of galaxies, although they are generally restricted to small fields-of-view. Using the large field-of-view of the DEIMOS multislit spectrograph on Keck and our Stellar Kinematics using Multiple Slits (SKiMS) technique, we are able to extract sky-subtracted stellar light spectra to large galactocentric radii. Here we present a new DEIMOS mask design named SuperSKiMS that explores large spatial scales without sacrificing high spatial sampling. We simulate a set of observations with such a mask design on the nearby galaxy NGC 1023, showing that the kinematic and metallicity measurements can reach radii where the galaxy surface brightness is several orders of magnitude fainter than the sky. Such a technique is also able to reproduce the kinematic and metallicity 2D distributions obtained from literature integral field spectroscopy in the innermost galaxy regions. In particular, we use the simulated NGC 1023 kinematics to model its total mass distribution to large radii, obtaining comparable results with those from published integral field unit observation. Finally, from new spectra of NGC 1023 we obtain stellar 2D kinematics and metallicity distributions that show good agreement with integral field spectroscopy results in the overlapping regions. In particular, we do not find a significant offset between our SKiMS and the ATLAS$^\rm{3D}$ stellar velocity dispersion at the same spatial locations.