We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying substructures within a galaxy, including embedded and counter-rotating discs, that are easily missed in unilateral studies of the photometry alone. Conversely, bars are rarely apparent in the stellar kinematics but are readily identified on images. Consequently, we distinguish the spin-ellipticity radial tracks of seven morpho-dynamical types: elliptical, lenticular, early spiral, late spiral, barred spiral, embedded disc, and 2-sigma galaxies. The importance of probing beyond the inner radii of galaxies is highlighted by the characteristics of galactic features in the spin-ellipticity radial tracks present at larger radii. The density of information presented through spin-ellipticity radial tracks emphasises a clear advantage to representing galaxies as a track, rather than a single point, in spin-ellipticity parameter space.

We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$\alpha$ gas velocity dispersion for our sample to be in the range $\bar{\sigma}_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $\Delta \bar{\sigma}_v / \bar{\sigma}_v = - 0.29 \pm 0.18$ for the H$\alpha$ flux-weighted mean velocity dispersion.

We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}{\delta}$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}{\delta}$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster.

The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm.

A primary goal of integral field spectroscopic (IFS) surveys is to provide a statistical census of galaxies classified by their internal kinematics. As a result, the observational spin parameter, $\lambda_R$, has become one of the most popular methods of quantifying the relative importance of velocity dispersion and rotation in supporting a galaxy's inner structure. The goal of this paper is to examine the relationship between the observationally deduced $\lambda_R$ and one of the most commonly used theoretical spin parameters in the literature, the Bullock et al. (2001) $\lambda'$. Using a set of $N$-body realisations of galaxies from which we construct mock IFS observations, we measure $\lambda_R$ as an observer would, incorporating the effects of beam smearing and seeing conditions. Assuming parameters typical of current IFS surveys, we confirm that there are strong positive correlations between $\lambda_R$ and measurement radius, and strong negative correlations between $\lambda_R$ and size of the PSF, for late-type galaxies; these biases can be reduced using a recently proposed empirical correction. Once observational biases are corrected for, we find that $\lambda_R$ provides a good approximation to $\sim \sqrt{3}/2 \; \lambda'(\rm R_{\rm eff})$, where $\lambda'$ is evaluated for the galactic stellar component within 1 R$_{\rm eff}$.

We present the second major release of data from the SAMI Galaxy Survey. Data Release Two includes data for 1559 galaxies, about 50% of the full survey. Galaxies included have a redshift range 0.004 < z < 0.113 and a large stellar mass range 7.5 < log (M_star/M_sun) < 11.6. The core data for each galaxy consist of two primary spectral cubes covering the blue and red optical wavelength ranges. For each primary cube we also provide three spatially binned spectral cubes and a set of standardised aperture spectra. For each core data product we provide a set of value-added data products. This includes all emission line value-added products from Data Release One, expanded to the larger sample. In addition we include stellar kinematic and stellar population value-added products derived from absorption line measurements. The data are provided online through Australian Astronomical Optics' Data Central. We illustrate the potential of this release by presenting the distribution of ~350,000 stellar velocity dispersion measurements from individual spaxels as a function of R/R_e, divided in four galaxy mass bins. In the highest stellar mass bin (log (M_star/M_sun)>11), the velocity dispersion strongly increases towards the centre, whereas below log (M_star/M_sun)<10 we find no evidence for a clear increase in the central velocity dispersion. This suggests a transition mass around log (M_star/M_sun) ~10 for galaxies with or without a dispersion-dominated bulge.

We introduce ProFound, a source finding and image analysis package. ProFound provides methods to detect sources in noisy images, generate segmentation maps identifying the pixels belonging to each source, and measure statistics like flux, size and ellipticity. These inputs are key requirements of ProFit, our recently released galaxy profiling package, where the design aim is that these two software packages will be used in unison to semi-automatically profile large samples of galaxies. The key novel feature introduced in ProFound is that all photometry is executed on dilated segmentation maps that fully contain the identifiable flux, rather than using more traditional circular or ellipse based photometry. Also, to be less sensitive to pathological segmentation issues, the de-blending is made across saddle points in flux. We apply ProFound in a number of simulated and real world cases, and demonstrate that it behaves reasonably given its stated design goals. In particular, it offers good initial parameter estimation for ProFit, and also segmentation maps that follow the sometimes complex geometry of resolved sources, whilst capturing nearly all of the flux. A number of bulge-disc decomposition projects are already making use of the ProFound and ProFit pipeline, and adoption is being encouraged by publicly releasing the software for the open source R data analysis platform under an LGPL-3 license on GitHub (github.com/asgr/ProFound).

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.

We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors.

Using the stellar kinematic maps and ancillary imaging data from the Sydney AAO Multi Integral field (SAMI) Galaxy Survey, the intrinsic shape of kinematically-selected samples of galaxies is inferred. We implement an efficient and optimised algorithm to fit the intrinsic shape of galaxies using an established method to simultaneously invert the distributions of apparent ellipticities and kinematic misalignments. The algorithm output compares favourably with previous studies of the intrinsic shape of galaxies based on imaging alone and our re-analysis of the ATLAS3D data. Our results indicate that most galaxies are oblate axisymmetric. We show empirically that the intrinsic shape of galaxies varies as a function of their rotational support as measured by the "spin" parameter proxy Lambda_Re. In particular, low spin systems have a higher occurrence of triaxiality, while high spin systems are more intrinsically flattened and axisymmetric. The intrinsic shape of galaxies is linked to their formation and merger histories. Galaxies with high spin values have intrinsic shapes consistent with dissipational minor mergers, while the intrinsic shape of low-spin systems is consistent with dissipationless multi-merger assembly histories. This range in assembly histories inferred from intrinsic shapes is broadly consistent with expectations from cosmological simulations.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

We examine the kinematic morphology of early-type galaxies (ETGs) in eight galaxy clusters in the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The clusters cover a mass range of 14.2<log(M_200/M_odot)<15.2 and we measure spatially-resolved stellar kinematics for 315 member galaxies with stellar masses 10.0<log(M_*/M_odot)<11.7 within 1R_200 of the cluster centers. We calculate the spin parameter, lambda_R and use that to classify the kinematic morphology of the galaxies as fast or slow rotators. The total fraction of slow rotators in the early-type galaxy population, F_SR=0.14+/-0.02 and does not depend on host cluster mass. Across the eight clusters, the fraction of slow rotators increases with increasing local overdensity. We also find that the slow-rotator fraction increases at small clustercentric radii (R_cl<0.3R_200), and note that there is also an increase in slow-rotator fraction at R_cl~0.6R_200. The slow rotators at these larger radii reside in cluster substructure. We find the strongest increase in slow-rotator fraction occurs with increasing stellar mass. After accounting for the strong correlation with stellar mass, we find no significant relationship between spin parameter and local overdensity in the cluster environment. We conclude that the primary driver for the kinematic morphology--density relationship in galaxy clusters is the changing distribution of galaxy stellar mass with local environment. The presence of slow rotators in substructure suggests that the cluster kinematic morphology--density relationship is a result of mass segregation of slow-rotating galaxies forming in groups that later merge with clusters and sink to the cluster center via dynamical friction.

We present ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level C++ library (libprofit), accessible via a command-line interface and documented API, along with high-level R (ProFit) and Python (PyProFit) interfaces (available at github.com/ICRAR/ libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit respectively). R ProFit is also available pre-built from CRAN, however this version will be slightly behind the latest GitHub version. libprofit offers fast and accurate two- dimensional integration for a useful number of profiles, including Sersic, Core-Sersic, broken-exponential, Ferrer, Moffat, empirical King, point-source and sky, with a simple mechanism for adding new profiles. We show detailed comparisons between libprofit and GALFIT. libprofit is both faster and more accurate than GALFIT at integrating the ubiquitous Serrsic profile for the most common values of the Serrsic index n (0.5 < n < 8). The high-level fitting code ProFit is tested on a sample of galaxies with both SDSS and deeper KiDS imaging. We find good agreement in the fit parameters, with larger scatter in best-fit parameters from fitting images from different sources (SDSS vs KiDS) than from using different codes (ProFit vs GALFIT). A large suite of Monte Carlo-simulated images are used to assess prospects for automated bulge-disc decomposition with ProFit on SDSS, KiDS and future LSST imaging. We find that the biggest increases in fit quality come from moving from SDSS- to KiDS-quality data, with less significant gains moving from KiDS to LSST.

Recent cosmological hydrodynamical simulations suggest that integral field spectroscopy can connect the high-order stellar kinematic moments h3 (~skewness) and h4 (~kurtosis) in galaxies to their cosmological assembly history. Here, we assess these results by measuring the stellar kinematics on a sample of 315 galaxies, without a morphological selection, using 2D integral field data from the SAMI Galaxy Survey. A proxy for the spin parameter ($\lambda_{R_e}$) and ellipticity ($\epsilon_e$) are used to separate fast and slow rotators; there exists a good correspondence to regular and non-regular rotators, respectively, as also seen in earlier studies. We confirm that regular rotators show a strong h3 versus $V/\sigma$ anti-correlation, whereas quasi-regular and non-regular rotators show a more vertical relation in h3 and $V/\sigma$. Motivated by recent cosmological simulations, we develop an alternative approach to kinematically classify galaxies from their individual h3 versus $V/\sigma$ signatures. We identify five classes of high-order stellar kinematic signatures using Gaussian mixture models. Class 1 corresponds to slow rotators, whereas Classes 2-5 correspond to fast rotators. We find that galaxies with similar $\lambda_{R_e}-\epsilon_e$ values can show distinctly different h3-$V/\sigma$ signatures. Class 5 objects are previously unidentified fast rotators that show a weak h3 versus $V/\sigma$ anti-correlation. These objects are predicted to be disk-less galaxies formed by gas-poor mergers. From morphological examination, however, there is evidence for large stellar disks. Instead, Class 5 objects are more likely disturbed galaxies, have counter-rotating bulges, or bars in edge-on galaxies. Finally, we interpret the strong anti-correlation in h3 versus $V/\sigma$ as evidence for disks in most fast rotators, suggesting a dearth of gas-poor mergers among fast rotators.

We investigate the relationship between stellar and gas specific angular momentum $j$, stellar mass $M_{*}$ and optical morphology for a sample of 488 galaxies extracted from the SAMI Galaxy Survey. We find that $j$, measured within one effective radius, monotonically increases with $M_{*}$ and that, for $M_{*}>$10$^{9.5}$ M$_{\odot}$, the scatter in this relation strongly correlates with optical morphology (i.e., visual classification and S\'ersic index). These findings confirm that massive galaxies of all types lie on a plane relating mass, angular momentum and stellar light distribution, and suggest that the large-scale morphology of a galaxy is regulated by its mass and dynamical state. We show that the significant scatter in the $M_{*}-j$ relation is accounted for by the fact that, at fixed stellar mass, the contribution of ordered motions to the dynamical support of galaxies varies by at least a factor of three. Indeed, the stellar spin parameter (quantified via $\lambda_R$) correlates strongly with S\'ersic and concentration indices. This correlation is particularly strong once slow-rotators are removed from the sample, showing that late-type galaxies and early-type fast rotators form a continuous class of objects in terms of their kinematic properties.

Observations consistently show that elliptical galaxies follow a tight "fundamental plane" scaling relation between size, mean surface brightness and velocity dispersion, with the form $R \propto \sigma^{a}\mu^{b}$. This relation not only has very small (<0.05 dex) intrinsic scatter, but also has significantly different coefficients from the expected virial scaling (a "tilt"). We analyze hundreds of simulations of elliptical galaxies formed from mergers of spiral galaxies in groups to determine if the fundamental plane can emerge from multiple, mostly minor and hierarchical collisionless mergers. We find that these simulated ellipticals lie on a similar fundamental plane with $a \approx 1.7$ and $b \approx 0.3$. The scatter about this plane is not larger than observed, while the tilt is in the correct sense, although $a$ is larger than for typical observations. This supports the idea that collisionless mergers can contribute significantly to the tilt of the fundamental plane. The tilt is mainly driven by a mass-dependent dark matter fraction, such that more massive galaxies have larger dark matter fractions within $R_{e}$. We further discuss the origin of this mass-dependent dark matter fraction and its compatibility with strong lensing observations, as well as the links between the fundamental plane, dynamical masses and the virial theorem.

We investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent ``ring" in phase space. Using several zoom simulations of clusters we show that the coherent distribution of the poststarbursts can be reasonably well-reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if satellite quenching occurs on a rapid timescale (0.1 < tau_{Q} < 0.5 Gyr) after galaxies make their first passage of R ~ 0.5R_{200}, a process that takes a total time of ~ 1 Gyr after first infall. We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (tau_{Q} = 0.4^{+0.3}_{-0.4} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid, and occurs well within R_{200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ~ 1 make it difficult to determine if hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two.

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $\Lambda$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping.

We present collisionless simulations of dry mergers in groups of three to twenty-five galaxies to test the hypothesis that elliptical galaxies form at the centers of such groups. Mock observations of the central remnants confirm their similarity to ellipticals, despite having no dissipational component. We vary the profile of the original spiral's bulge and find that ellipticals formed from spirals with exponential bulges have too low Sersic indices. Mergers of spirals with de Vaucouleurs (classical) bulges produce remnants with larger Sersic indices correlated with luminosity, as with SDSS ellipticals. Exponential bulge mergers are better fits to faint ellipticals, whereas classical bulge mergers better match luminous ellipticals. Similarly, luminous ellipticals are better reproduced by remnants undergoing many (>5) mergers, and fainter ellipticals by those with fewer mergers. The remnants follow tight size-luminosity and velocity dispersion-luminosity (Faber-Jackson) relations (<0.12 dex scatter), demonstrating that stochastic merging can produce tight scaling relations if the merging galaxies also follow tight scaling relations. The slopes of the size-luminosity and Faber-Jackson relations are close to observations but slightly shallower in the former case. Both relations' intercepts are offset - remnants are too large but have too low dispersions at fixed luminosity. Some remnants show substantial (v/\sigma > 0.1) rotational support, although most are slow rotators and few are very fast rotators (v/\sigma > 0.5). These findings contrast with previous studies concluding that dissipation necessary to produce ellipticals from binary mergers of spirals. Multiple, mostly minor and dry mergers can produce bright ellipticals, whereas significant dissipation could be required to produce faint, rapidly-rotating ellipticals.

We present a series of hundreds of collisionless simulations of galaxy group mergers. These simulations are designed to test whether the properties of elliptical galaxies - including the key fundamental plane scaling relation, morphology and kinematics - can be simultaneously reproduced by dry multiple mergers in galaxy groups. Preliminary results indicate that galaxy group mergers can produce elliptical remnants lying on a tilted fundamental plane, even without a central dissipational component from a starburst. This suggests that multiple mergers in groups are an alternate avenue for the formation of elliptical galaxies which could well dominate for luminous ellipticals.