We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $\Delta \log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies.

We construct Tully-Fisher relations (TFRs), from large samples of galaxies with spatially-resolved H$\alpha$ emission maps from the K-band Multi-Object Spectrograph (KMOS) Redshift One Spectroscopic Survey (KROSS) at $z\approx1$. We compare these to data from the Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph (SAMI) Galaxy Survey at $z\approx0$. We stringently match the data quality of the latter to the former, and apply identical analysis methods and sub-sample selection criteria to both to conduct a direct comparison of the absolute $K$-band magnitude and stellar mass TFRs at $z\approx1$ and $z\approx0$. We find that matching the quality of the SAMI data to that of KROSS results in TFRs that differ significantly in slope, zero-point and (sometimes) scatter in comparison to the corresponding original SAMI relations. These differences are in every case as large or larger than the differences between the KROSS $z\approx1$ and matched SAMI $z\approx0$ relations. Accounting for these differences, we compare the TFRs at $z\approx1$ and $z\approx0$. For disk-like, star-forming galaxies we find no significant difference in the TFR zero-points between the two epochs. This suggests the growth of stellar mass and dark matter in these types of galaxies is intimately linked over this $\approx8$ Gyr period.

In order to determine the causes of kinematic asymmetry in the H$\alpha$ gas in the SAMI Galaxy Survey sample, we investigate the comparative influences of environment and intrinsic properties of galaxies on perturbation. We use spatially resolved H$\alpha$ velocity fields from the SAMI Galaxy Survey to quantify kinematic asymmetry ($\overline{v_{asym}}$) in nearby galaxies and environmental and stellar mass data from the GAMA survey. {We find that local environment, measured as distance to nearest neighbour, is inversely correlated with kinematic asymmetry for galaxies with $\mathrm{\log(M_*/M_\odot)}>10.0$, but there is no significant correlation for galaxies with $\mathrm{\log(M_*/M_\odot)}<10.0$. Moreover, low mass galaxies ($\mathrm{\log(M_*/M_\odot)}<9.0$) have greater kinematic asymmetry at all separations, suggesting a different physical source of asymmetry is important in low mass galaxies.} We propose that secular effects derived from gas fraction and gas mass may be the primary causes of asymmetry in low mass galaxies. High gas fraction is linked to high $\frac{\sigma_{m}}{V}$ (where $\sigma_m$ is H$\alpha$ velocity dispersion and $V$ the rotation velocity), which is strongly correlated with $\overline{v_{asym}}$, and galaxies with $\log(M_*/M_\odot)<9.0$ have offset $\overline{\frac{\sigma_{m}}{V}}$ from the rest of the sample. Further, asymmetry as a fraction of dispersion decreases for galaxies with $\log(M_*/M_\odot)<9.0$. Gas mass and asymmetry are also inversely correlated in our sample. We propose that low gas masses in dwarf galaxies may lead to asymmetric distribution of gas clouds, leading to increased relative turbulence.

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 investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use 'simulated slits' to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis.

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.

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K.

We present an analysis of the global stellar populations of galaxies in the SAMI Galaxy Survey. Our sample consists of 1319 galaxies spanning four orders of magnitude in stellar mass and includes all morphologies and environments. We derive luminosity-weighted, single stellar population equivalent stellar ages, metallicities and alpha enhancements from spectra integrated within one effective radius apertures. Variations in galaxy size explain the majority of the scatter in the age--mass and metallicity--mass relations. Stellar populations vary systematically in the plane of galaxy size and stellar mass, such that galaxies with high stellar surface mass density are older, more metal-rich and alpha-enhanced than less dense galaxies. Galaxies with high surface mass densities have a very narrow range of metallicities, however, at fixed mass, the spread in metallicity increases substantially with increasing galaxy size (decreasing density). We identify residual correlations with morphology and environment. At fixed mass and size, galaxies with late-type morphologies, small bulges and low Sersic n are younger than early-type, high n, high bulge-to-total galaxies. Age and metallicity both show small residual correlations with environment; at fixed mass and size, galaxies in denser environments or more massive halos are older and somewhat more metal rich than those in less dense environments. We connect these trends to evolutionary tracks within the size--mass plane.

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.

The fraction of galaxies supported by internal rotation compared to galaxies stabilized by internal pressure provides a strong constraint on galaxy formation models. In integral field spectroscopy surveys, this fraction is biased because survey instruments typically only trace the inner parts of the most massive galaxies. We present aperture corrections for the two most widely used stellar kinematic quantities $V/\sigma$ and $\lambda_{R}$. Our demonstration involves integral field data from the SAMI Galaxy Survey and the ATLAS$^{\rm{3D}}$ Survey. We find a tight relation for both $V/\sigma$ and $\lambda_{R}$ when measured in different apertures that can be used as a linear transformation as a function of radius, i.e., a first-order aperture correction. We find that $V/\sigma$ and $\lambda_{R}$ radial growth curves are well approximated by second order polynomials. By only fitting the inner profile (0.5$R_{\rm{e}}$), we successfully recover the profile out to one $R_{\rm{e}}$ if a constraint between the linear and quadratic parameter in the fit is applied. However, the aperture corrections for $V/\sigma$ and $\lambda_{R}$ derived by extrapolating the profiles perform as well as applying a first-order correction. With our aperture-corrected $\lambda_{R}$ measurements, we find that the fraction of slow rotating galaxies increases with stellar mass. For galaxies with $\log M_{*}/M_{\odot}>$ 11, the fraction of slow rotators is $35.9\pm4.3$ percent, but is underestimated if galaxies without coverage beyond one $R_{\rm{e}}$ are not included in the sample ($24.2\pm5.3$ percent). With measurements out to the largest aperture radius the slow rotator fraction is similar as compared to using aperture corrected values ($38.3\pm4.4$ percent). Thus, aperture effects can significantly bias stellar kinematic IFS studies, but this bias can now be removed with the method outlined here.

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 evidence for halo assembly bias as a function of geometric environment. By classifying GAMA galaxy groups as residing in voids, sheets, filaments or knots using a tidal tensor method, we find that low-mass haloes that reside in knots are older than haloes of the same mass that reside in voids. This result provides direct support to theories that link strong halo tidal interactions with halo assembly times. The trend with geometric environment is reversed at large halo mass, with haloes in knots being younger than haloes of the same mass in voids. We find a clear signal of halo downsizing - more massive haloes host galaxies that assembled their stars earlier. This overall trend holds independently of geometric environment. We support our analysis with an in-depth exploration of the L-Galaxies semi-analytic model, used here to correlate several galaxy properties with three different definitions of halo formation time. We find a complex relationship between halo formation time and galaxy properties, with significant scatter. We confirm that stellar mass to halo mass ratio, specific star-formation rate and mass-weighed age are reasonable proxies of halo formation time, especially at low halo masses. Instantaneous star-formation rate is a poor indicator at all halo masses. Using the same semi-analytic model, we create mock spectral observations using complex star-formation and chemical enrichment histories, that approximately mimic GAMA's typical signal-to-noise and wavelength range. We use these mocks to assert how well potential proxies of halo formation time may be recovered from GAMA-like spectroscopic data.

We describe the selection of galaxies targeted in eight low redshift clusters (APMCC0917, A168, A4038, EDCC442, A3880, A2399, A119 and A85; $0.029 < z < 0.058$) as part of the Sydney-AAO Multi-Object integral field Spectrograph Galaxy Survey (SAMI-GS). We have conducted a redshift survey of these clusters using the AAOmega multi-object spectrograph on the 3.9m Anglo-Australian Telescope. The redshift survey is used to determine cluster membership and to characterise the dynamical properties of the clusters. In combination with existing data, the survey resulted in 21,257 reliable redshift measurements and 2899 confirmed cluster member galaxies. Our redshift catalogue has a high spectroscopic completeness ($\sim 94\%$) for $r_{\rm petro} \leq 19.4$ and clustercentric distances $R< 2\rm{R}_{200}$. We use the confirmed cluster member positions and redshifts to determine cluster velocity dispersion, $\rm{R}_{200}$, virial and caustic masses, as well as cluster structure. The clusters have virial masses $14.25 \leq {\rm log }({\rm M}_{200}/\rm{M}_{\odot}) \leq 15.19$. The cluster sample exhibits a range of dynamical states, from relatively relaxed-appearing systems, to clusters with strong indications of merger-related substructure. Aperture- and PSF-matched photometry are derived from SDSS and VST/ATLAS imaging and used to estimate stellar masses. These estimates, in combination with the redshifts, are used to define the input target catalogue for the cluster portion of the SAMI-GS. The primary SAMI-GS cluster targets have $R< \rm{R}_{200}$, velocities $|v_{\rm pec}| < 3.5\sigma_{200}$ and stellar masses $9.5 \leq {\rm log(M}^*_{approx}/\rm{M}_{\odot}) \leq 12$. Finally, we give an update on the SAMI-GS progress for the cluster regions.

We investigate the properties of the stellar populations of model galaxies as a function of galaxy evolutionary history and angular momentum content. We use the new semi-analytic model presented in Tonini et al. (2016). This new model follows the angular momentum evolution of gas and stars, providing the base for a new star formation recipe, and treatment of the effects of mergers that depends on the central galaxy dynamical structure. We find that the new recipes have the effect of boosting the efficiency of the baryonic cycle in producing and recycling metals, as well as preventing minor mergers from diluting the metallicity of bulges and ellipticals. The model reproduces the stellar mass - stellar metallicity relation for galaxies above 1e10 solar masses, including Brightest Cluster Galaxies. Model disks, galaxies dominated by instability-driven components, and merger-driven objects each stem from different evolutionary channels. These model galaxies therefore occupy different loci in the galaxy mass-size relation, which we find to be in accord with the Atlas 3D classification of disk galaxies, fast rotators and slow rotators. We find that the stellar populations' properties depend on the galaxy evolutionary type, with more evolved stellar populations being part of systems that have lost or dissipated more angular momentum during their assembly history.

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.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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.

We present a theoretical model for the evolution of mass, angular momentum and size of galaxy disks and bulges, and we implement it into the semi-analytic galaxy formation code SAGE. The model follows both secular and violent evolutionary channels, including smooth accretion, disk instabilities, minor and major mergers. We find that the combination of our recipe with hierarchical clustering produces two distinct populations of bulges: merger-driven bulges, akin to classical bulges and ellipticals, and instability-driven bulges, akin to secular (or pseudo-)bulges. The model mostly reproduces the mass-size relation of gaseous and stellar disks, the evolution of the mass-size relation of ellipticals, the Faber-Jackson relation, and the magnitude-colour diagram of classical and secular bulges. The model predicts only a small overlap of merger-driven and instability-driven components in the same galaxy, and predicts different bulge types as a function of galaxy mass and disk fraction. Bulge type also affects the star formation rate and colour at a given luminosity. The model predicts a population of merger-driven red ellipticals that dominate both the low-mass and high-mass ends of the galaxy population, and span all dynamical ages; merger-driven bulges in disk galaxies are dynamically old and do not interfere with subsequent evolution of the star-forming component. Instability-driven bulges dominate the population at intermediate galaxy masses, especially thriving in massive disks. The model green valley is exclusively populated by instability-driven bulge hosts. Through the present implementation the mass accretion history is perceivable in the galaxy structure, morphology and colours.

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

Intensity mapping of the neutral hydrogen (HI) is a new observational tool that can be used to efficiently map the large-scale structure of the Universe over wide redshift ranges. The power spectrum of the intensity maps contains cosmological information on the matter distribution and probes galaxy evolution by tracing the HI content of galaxies at different redshifts and the scale-dependence of HI clustering. The cross-correlation of intensity maps with galaxy surveys is a robust measure of the power spectrum which diminishes systematics caused by instrumental effects and foreground removal. We examine the cross-correlation signature at redshift z=0.9 using a variant of the semi-analytical galaxy formation model SAGE (Croton et al. 2016) applied to the Millennium simulation in order to model the HI gas of galaxies as well as their optical magnitudes based on their star-formation history. We determine the clustering of the cross-correlation power for different types of galaxies determined by their colours, acting as a proxy for their star-formation activity. We find that the cross-correlation coefficient for red quiescent galaxies falls off more quickly on smaller scales k>0.2h/Mpc than for blue star-forming galaxies. Additionally, we create a mock catalogue of highly star-forming galaxies using a selection function to mimic the WiggleZ survey, and use this to predict existing and future cross-correlation measurements of the GBT and Parkes telescope. We find that the cross-power of highly star-forming galaxies shows a higher clustering on small scales than any other galaxy type and that this significantly alters the power spectrum shape on scales k>0.2h/Mpc. We show that the cross-correlation coefficient is not negligible when interpreting the cosmological cross-power spectrum. On the other hand, it contains information about the HI content of the optically selected galaxies.

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.

We present a methodology for the regularisation and combination of sparse sampled and irregularly gridded observations from fibre-optic multi-object integral-field spectroscopy. The approach minimises interpolation and retains image resolution on combining sub-pixel dithered data. We discuss the methodology in the context of the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey underway at the Anglo-Australian Telescope. The SAMI instrument uses 13 fibre bundles to perform high-multiplex integral-field spectroscopy across a one degree diameter field of view. The SAMI Galaxy Survey is targeting 3000 galaxies drawn from the full range of galaxy environments. We demonstrate the subcritical sampling of the seeing and incomplete fill factor for the integral-field bundles results in only a 10% degradation in the final image resolution recovered. We also implement a new methodology for tracking covariance between elements of the resulting datacubes which retains 90% of the covariance information while incurring only a modest increase in the survey data volume.

We take advantage of the first data from the Sydney-AAO Multi-object Integral field (SAMI) Galaxy Survey to investigate the relation between the kinematics of gas and stars, and stellar mass in a comprehensive sample of nearby galaxies. We find that all 235 objects in our sample, regardless of their morphology, lie on a tight relation linking stellar mass ($M_{*}$) to internal velocity quantified by the $S_{0.5}$ parameter, which combines the contribution of both dispersion ($\sigma$) and rotational velocity ($V_{rot}$) to the dynamical support of a galaxy ($S_{0.5}=\sqrt{0.5V_{rot}^{2}+\sigma^{2}}$). Our results are independent of the baryonic component from which $\sigma$ and $V_{rot}$ are estimated, as the $S_{0.5}$ of stars and gas agree remarkably well. This represents a significant improvement compared to the canonical $M_{*}$ vs. $V_{rot}$ and $M_{*}$ vs. $\sigma$ relations. Not only is no sample pruning necessary, but also stellar and gas kinematics can be used simultaneously, as the effect of asymmetric drift is taken into account once $V_{rot}$ and $\sigma$ are combined. Our findings illustrate how the combination of dispersion and rotational velocities for both gas and stars can provide us with a single dynamical scaling relation valid for galaxies of all morphologies across at least the stellar mass range 8.5$<log(M_{*}/M_{\odot})<$11. Such relation appears to be more general and at least as tight as any other dynamical scaling relation, representing a unique tool for investigating the link between galaxy kinematics and baryonic content, and a less biased comparison with theoretical models.

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

(Abridged) We investigate the morphology dependence of the Tully-Fisher relation, and the expansion of the relation into a three-dimensional manifold defined by luminosity, total circular velocity and a third dynamical parameter, to fully characterise spiral galaxies across all morphological types. We use a semi-analytic hierarchical model of galaxy evolution to build the theoretical TF relation. With this tool, we analyse a unique dataset of galaxies for which we cross-match luminosity with total circular velocity and central velocity dispersion. We provide a theoretical framework to calculate such measurable quantities from semi-analytic models. We establish the morphology dependence of the TF relation in both model and data. We analyse the dynamical properties of the model galaxies and determine that the parameter 'sigma/Vc', i.e. the ratio between random and total motions defined by velocity dispersion and circular velocity, accurately characterises the varying slope of the TF relation for different model galaxy types. We apply these dynamical cuts to the observed galaxies and find that such selection produces a differential slope of the TF relation. The TF slope in different ranges of 'sigma/Vc' is consistent with that for the photometric classification in Sa, Sb, Sc. We conclude that 'sigma/Vc' is a good parameter to classify galaxy type, and we argue that such classification based on dynamics more closely mirrors the physical properties of the observed galaxies, compared to visual classification. We also argue that dynamical classification is useful for samples where eye inspection is not reliable or impractical. We conclude that 'sigma/Vc' is a suitable parameter to characterise the hierarchical assembly history that determines the disk-to-bulge ratio, and to expand the TF relation into a three-dimensional manifold, defined by luminosity, circular velocity and 'sigma/Vc'.

(Abridged) In a recent work we explored the dependence of galaxy stellar population properties derived from broad-band spectral energy distribution fitting on the fitting parameters, e.g. SFHs, age grid, metallicity, IMF, dust reddening, reddening law, filter setup and wavelength coverage. In this paper we consider also redshift as a free parameter in the fit and study whether one can obtain reasonable estimates of photometric redshifts and stellar population properties at once. We use mock star-forming as well as passive galaxies placed at various redshifts (0.5 to 3) as test particles. Mock star-forming galaxies are extracted from a semi-analytical galaxy formation model. We show that for high-z star-forming galaxies photometric redshifts, stellar masses and reddening can be determined simultaneously when using a broad wavelength coverage and a wide template setup in the fit. Masses are similarly well recovered (median ~ 0.2 dex) as at fixed redshift. For old galaxies with little recent star formation masses are better recovered than in the fixed redshift case, such that the median recovered stellar mass improves by up to 0.3 dex whereas the uncertainty in the redshift accuracy increases by only ~ 0.05. However, a failure in redshift recovery also means a failure in mass recovery. As at fixed redshift mismatches in SFH and degeneracies between age, dust and now also redshift cause underestimated ages, overestimated reddening and underestimated masses. Stellar masses are best determined at low redshift without reddening in the fit (median underestimation ~ 0.1 dex for similarly well recovered redshifts). Not surprisingly, the recovery of properties is substantially better for passive galaxies. In all cases, the recovery of physical parameters is crucially dependent on the wavelength coverage adopted in the fitting. Scaling relations for the transformation of stellar masses are provided.

We address the fundamental question of matching the rest-frame K-band luminosity function (LF) of galaxies over the Hubble time using semi-analytic models, after modification of the stellar population modelling. We include the Maraston evolutionary synthesis models, that feature a higher contribution by the Thermally Pulsating - Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and the Menci model. We leave all other input physics and parameters unchanged. We find that the modification of the stellar population emission can solve the mismatch between models and the observed rest-frame K-band luminosity from the brightest galaxies derived from UKIDSS data at high redshift. For all explored semi-analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy was recently pointed out. The reason for the success is that at these cosmic epochs the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB phase which makes them redder without the need of changing their mass or age. At the same time, the known overestimation of the faint end is enhanced in the K-band when including the TP-AGB contribution. At lower redshifts (z<2) some of the explored models deviate from the data. This is due to too short merging timescales and inefficient 'radio-mode' AGN feedback. Our results show that a strong evolution in mass predicted by hierarchical models is compatible with no evolution on the bright-end of the K-band LF from z=3 to the local universe. This means that, at high redshifts and contrary to what is commonly accepted, K-band emission is not necessarily a good tracer of galaxy mass.

(Abridged) We build a theoretical model to study the origin of the globular cluster metallicity bimodality in the hierarchical galaxy assembly scenario, based on the observed galaxy mass-[O/H] relation and the galaxy stellar mass function up to z ~4, and on theoretical merger rates. We derive a new galaxy [Fe/H]-M(star) relation as a function of z, and by assuming that GCs share the metallicity of their parent galaxy when they form, we populate the merger tree with GCs. We perform a series of Monte-Carlo simulations of the galaxy assembly, and study the properties of the final GC population as a function of galaxy mass, assembly and star formation history, and under different assumptions for the evolution of the galaxy mass-[Fe/H] relation. The main results are: 1) The hierarchical clustering scenario naturally predicts a metallicity bimodality in the galaxy GC population: the metal-rich GCs are formed in the galaxy main progenitor around z~2, and the metal-poor GCs are accreted from satellites and formed at z~3-4. 2) The model reproduces the observed relations for the metallicity of the metal-rich and metal-poor GCs as a function of galaxy mass. The positions of the metal-poor and metal-rich peaks depend exclusively on the evolution of the galaxy mass-[Fe/H] relation and the [O/Fe], both of which can be constrained by this method. We find that the galaxy [O/Fe] evolves linearly with z from a value of ~0.5 at z~4 to a value of ~0.1 at z=0. 3) Given a galaxy mass, the relative strength of the metal-rich and metal-poor peaks depends exclusively on the galaxy assembly and star formation history: galaxies in denser environments and/or early types galaxies show a larger fraction of metal-poor GCs, while galaxies with a sparse merger history and/or late type galaxies are dominated by metal-rich GCs. 4) The GC metallicity bimodality disappears for galaxy masses below M(star)~1e9, and for z>2.

We investigate the evolution of Brightest Cluster Galaxies (BCGs) from redshift z~1.6 to z~0. We use the semi-analytic model of Croton et al. (2006) with a new spectro-photometric model based on the Maraston (2005) stellar populations and a new recipe for the dust extinction. We compare the model predictions of the K-band luminosity evolution and the J-K, V-I and I-K colour evolution with a series of datasets, including Collins et al. (Nature, 2009) who argued that semi-analytic models based on the Millennium simulation cannot reproduce the red colours and high luminosity of BCGs at z>1. We show instead that the model is well in range of the observed luminosity and correctly reproduces the colour evolution of BCGs in the whole redshift range up to z~1.6. We argue that the success of the semi-analytic model is in large part due to the implementation of a more sophisticated spectro-photometric model. An analysis of the model BCGs shows an increase in mass by a factor ~2 since z~1, and star formation activity down to low redshifts. While the consensus regarding BCGs is that they are passively evolving, we argue that this conclusion is affected by the degeneracy between star formation history and stellar population models used in SED-fitting, and by the inefficacy of toy-models of passive evolution to capture the complexity of real galaxies, expecially those with rich merger histories like BCGs. Following this argument, we also show that in the semi-analytic model the BCGs show a realistic mix of stellar populations, and that these stellar populations are mostly old. In addition, the age-redshift relation of the model BCGs follows that of the universe, meaning that given their merger history and star formation history, the ageing of BCGs is always dominated by the ageing of their stellar populations. In a LambdaCDM universe, we define such evolution as "passive in the hierarchical sense".

(Abridged) We explore the dependence of galaxy stellar population properties derived from broad-band SED-fitting - such as age, stellar mass, dust reddening, etc. - on a variety of parameters, such as SFHs, metallicity, IMF, dust reddening and reddening law, and wavelength coverage. Mock galaxies serve as test particles. We confirm our earlier results based on real z=2 galaxies, that usually adopted \tau-models lead to overestimate the SFR and to underestimate the stellar mass. Here, we show that - for star-forming galaxies - ages, masses and reddening, can be well determined simultaneously only when the correct SFH is identified. This is the case for inverted-\tau-models at high-z, for which we find that the mass recovery (at fixed IMF) is as good as ~0.04 dex. Since the right SFH is usually unknown we quantify offsets generated by adopting standard fitting setups. Stellar masses are generally underestimated resulting from underestimating ages. For fitting setups with a variety of SFHs the median mass recovery at z ~ 2-3 is as decent as ~0.1 dex, albeit with large scatter. The situation worsens towards lower redshifts because of the variety of possible SFHs and ages (~0.6 dex at z=0.5). A practical trick to improve upon this is to exclude reddening from the fitting to avoid unrealistically young and dusty solutions. Reddening and SFRs should then be determined by a separate fit. As expected, the recovery of properties is better for passive galaxies. For the two galaxy types the parameter recovery is optimal for a wavelength coverage from the rest-frame UV to the rest-frame near-IR. We quantify the effect of narrowing the wavelength coverage or adding/removing filters which can be useful for planning observational surveys. Finally, we provide scaling relations that allow the transformation of stellar masses obtained using different template fitting setups and stellar population models.

We use the semi-analytic model GalICS to predict the Tully-Fisher relation in the B, I and for the first time, in the K band, and its evolution with redshift, up to z~1. We refined the determination of the disk galaxies rotation velocity, with a dynamical recipe for the rotation curve, rather than a simple conversion from the total mass to maximum velocity. The new recipe takes into account the disk shape factor, and the angular momentum transfer occurring during secular evolution leading to the formation of bulges. This produces model rotation velocities that are lower by ~20-25% for the majority of the spirals. We implemented stellar population models with a complete treatment of the TP-AGB, which leads to a revision of the mass-to-light ratio in the near-IR. I/K band luminosities increase by ~0.3/0.5 mags at redshift z=0 and by ~0.5/1 mags at z=3. With these two new recipes in place, the comparison between the predicted Tully-Fisher relation with a series of datasets in the optical and near-IR, at redshifts between 0 and 1, is used as a diagnostics of the assembly and evolution of spiral galaxies in the model. At 0.4<z<1.2 the match between the new model and data is remarkably good, especially for later-type spirals (Sb/Sc). At z=0 the new model shows a net improvement in comparison with its original version of 2003, and in accord with recent observations in the K band, the model Tully-Fisher also shows a morphological differentiation. However, in all bands the z=0 model Tully-Fisher is too bright. We argue that this behaviour is caused by inadequate star formation histories in the model galaxies at low redshifts. The star-formation rate declines too slowly, due to continuous gas infall that is not efficiently suppressed. An analysis of the model disk scale lengths, at odds with observations, hints to some missing physics in the modeling of disk formation inside dark matter halos.

Fitting synthetic spectral energy distributions (SED) to the multi-band photometry of galaxies to derive their star formation rates (SFR), stellar masses, ages, etc. requires making a priori assumptions about their star formation histories (SFH). A widely adopted parameterization of the SFH, the so-called tau-models where SFR goes as e^{-t/tau) is shown to lead to unrealistically low ages when applied to star forming galaxies at z ~ 2, a problem shared by other SFHs when the age is left as a free parameter in the fitting. This happens because the SED of such galaxies, at all wavelengths, is dominated by their youngest stellar populations, which outshine the older ones. Thus, the SED of such galaxies conveys little information on the beginning of star formation. To cope with this problem, we explore a variety of SFHs, such as constant SFR and inverted-tau models - with SFR as e^{+t/tau) - along with various priors on age, including assuming that star formation started at high redshift in all the galaxies. We find that inverted-tau models with such latter assumption give SFRs and extinctions in excellent agreement with the values derived using only the UV part of the SED. These models are also shown to accurately recover the SFRs and masses of mock galaxies at z ~ 2 constructed from semi-analytic models. All other explored SFH templates do not fulfil these two test. In particular, direct-tau models with unconstrained age in the fitting procedure overstimate SFRs and underestimate stellar mass, and would exacerbate an apparent mismatch between the cosmic evolution of the volume densities of SFR and stellar mass. We conclude that for high-redshift star forming galaxies an exponentially increasing SFR with a high formation redshift is preferable to other forms of the SFH so far adopted in the literature.

In a recent paper we presented the first semi-analytic model of galaxy formation in which the Thermally-Pulsing Asymptotic Giant Branch phase of stellar evolution has been fully implemented. Here we address the comparison with observations, and show how the TP-AGB recipe affects the performance of the model in reproducing the colours and near-IR luminosities of high-redshift galaxies. We find that the semi-analytic model with the TP-AGB better matches the colour-magnitude and colour-colour relations at z ~ 2, both for nearly-passive and for star-forming galaxies. The model with TP-AGB produces star-forming galaxies with red V-K colours, thus revising the unique interpretation of high-redshift red objects as 'red & dead'. We also show that without the TP-AGB the semi-analytic model fails at reproducing the observed colours, a situation that cannot be corrected by dust reddening. We also explore the effect of nebular emission on the predicted colour-magnitude relation of star-forming galaxies, to conclude that it does not play a significant role in reddening their colours, at least in the range of star-formation rates covered by the model. Finally, the rest-frame K-band luminosity function at z ~ 2.5 is more luminous by almost 1 magnitude. This indicates that the AGN feedback recipe that is adopted to regulate the high-mass end of the luminosity function should be sophisticated to take the effect of the stellar populations into account at high redshifts.

The spectro-photometric properties of galaxies in galaxy formation models are obtained by combining the predicted history of star formation and mass accretion with the physics of stellar evolution through stellar population models. In the recent literature, significant differences have emerged regarding the implementation of the Thermally-Pulsing Asymptotic Giant Branch phase of stellar evolution. The emission in the TP-AGB phase dominates the bolometric and near-IR spectrum of intermediate-age (~1 Gyr) stellar populations, hence it is crucial for the correct modeling of the galaxy luminosities and colours. In this paper for the first time, we incorporate a full prescription of the TP-AGB phase in a semi-analytic model of galaxy formation. We find that the inclusion of the TP-AGB in the model spectra dramatically alters the predicted colour-magnitude relation and its evolution with redshift. When the TP-AGB phase is active, the rest-frame V-K galaxy colours are redder by almost 2 magnitudes in the redshift range z~2-3 and by 1 magnitude at z~1. Very red colours are produced in disk galaxies, so that the V-K colour distributions of disk and spheroids are virtually undistinguishable at low redshifts. We also find that the galaxy K-band emission is more than 1 magnitude higher in the range z~1-3. This may alleviate the difficulties met by the hierarchical clustering scenario in predicting the red galaxy population at high redshifts. The comparison between simulations and observations have to be revisited in the light of our results.

The structure and dynamics of dark matter halos, as predicted by the hierarchical clustering scenario, are at odds with the properties inferred from the observations at galactic scales. My Thesis addresses this problem by taking an evolutionary approach. I analysed in detail the many and different observational evidences of a discrepancy the predicted halo equilibrium state and the one inferred from the measurable properties of disk galaxies, as well as of the scaling relations existing between the angular momentum, geometry and mass distribution of the luminous and dark components, and realized that they all seem to point towards the same conclusion: the baryons hosted inside the halo, by collapsing and assembling to form the galaxy, perturb the halo equilibrium structure and made it evolve into new configurations. From the theoretical point of view, the behaviour of dark matter halos as collisionless systems of particles makes their equilibrium structure and mass distribution extremely sensitive to perturbations of their inner dynamics. The galaxy formation occurring inside the halos is a tremendous event, and the dynamical coupling between the baryons and the dark matter during the protogalaxy collapse represents a perturbation of the halo dynamical structure large enough to trigger a halo evolution, according to the relative mass and angular momentum of the two components. My conclusion is that the structure and dynamics of dark matter halos, as well as the origin of the connection between the halo and galaxy properties, are to be understood in in terms of a joint evolution of the baryonic and dark components, originating at the epoch of the collapse and formation of the galaxy.

In the current LambdaCDM cosmological scenario, N-body simulations provide us with a Universal mass profile, and consequently a Universal equilibrium circular velocity of the virialized objects, as galaxies. In this paper we obtain, by combining kinematical data of their inner regions with global observational properties, the Universal Rotation Curve (URC) of disk galaxies and the corresponding mass distribution out to their virial radius. This curve extends the results of Paper I, concerning the inner luminous regions of Sb-Im spirals, out to the edge of the galaxy halos.

Aims. We present an analysis of 37 high-quality extended rotation curves that highlights the existence of a new discrepancy (or a new aspect of an old discrepancy) between the density profiles predicted by the Lambda Cold Dark Matter (LambdaCDM) theory and the actual distribution of dark matter in galaxies. Methods. We compare the predicted face-value density vs. enclosed mass relationship, at large distances, to the observational data at the last measured radii of the rotation curves and in two whole rotation curves of high quality. A further analysis is performed by studying a relation, inbuilt in LambdaCDM, that links at radius R, the enclosed halo mass M(R) and its density rho(R) in a way that is independent of the mass of the virialised object. Results. We find that the predicted density vs. enclosed mass relationship has a systematic offset with respect to the observational data. In test case extended rotation curves, at their last measured point, the predicted NFW densities are up to a factor 3 lower than those derived from the kinematics. Moreover, the abovementioned relation, inbuilt in LambdaCDM, does not hold for the objects of our sample. Such a new outer discrepancy is different and maybe complementary with respect to the core/cusp issue, for which the NFW densities turn out to be higher than those observed and it seems to imply a global mass rearrangement of a pristine NFW-LambdaCDM halo.

We propose that angular momentum transfer from the baryons to the Dark Matter (DM) during the early stages of galaxy formation can flatten the halo inner density profile and modify the halo dynamics. We compute the phase-space distribution function of DM halos, that corresponds to the density and anisotropy profiles obtained from N-body simulations in the concordance cosmology. We then describe an injection of angular momentum into the halo by modifying the distribution function, and show that the system evolves into a new equilibrium configuration; the latter features a constant central density and a tangentially-dominated anisotropy profile in the inner regions, while the structure is nearly unchanged beyond 10% of the virial radius. Then we propose a toy model to account for such a halo evolution, based on the angular momentum exchange due to dynamical friction; at the epoch of galaxy formation this is efficiently exerted by the DM onto the gas clouds spiralling down the potential well. The comparison between the angular momentum profile gained by the halo through dynamical friction and that provided by the perturbed distribution function reveals a surprising similarity, hinting at the reliability of the process.

We compute the angular momentum, the spin parameter and the related distribution function for Dark Matter halos hosting a spiral galaxy. We base on scaling laws, inferred from observations, that link the properties of the galaxy to those of the host halo; we further assume that the Dark Matter has the same total specific angular momentum of the baryons. Our main results are: (i) we find that the gas component of the disk significantly contributes to the total angular momentum of the system; (ii) by adopting for the Dark Matter the observationally supported Burkert profile, we compute the total angular momentum of the disk and its correlation with the rotation velocity; (iii) we find that the distribution function of the spin parameter $\lambda$ peaks at a value of about 0.03, consistent with a no-major-merger scenario for the late evolution of spiral galaxies.

We construct a generic mass model of a disk galaxy by combining the usual two mass components, an exponential thin stellar disk and a dark halo for which we alternatively assume one of the two most common mass distributions. We construct its generic rotation curve, with two free parameters linked to its shape and one to its amplitude, namely the halo characteristic scale lenght, the fractional amount of Dark Matter and the total mass at one disk scale lenght. We find that there is no need of a fine-tuning in the above parameters to produce smooth, featureless rotation curves like the ones actually observed. The rotation curves obtained from different mass distributions are significanlty affected by the variation of the parameters, and show substantial difference in their slope; however, they do not feature irregularities marking the transition from a disk-dominated to a halo-dominated regime. On the other hand we show, by direct analysis of a large number of plausible and different mass models, that an observed rotation curve fulfilling well defined quality requirements can be uniquely and properly decomposed in its dark and luminous components; through a suitable mass modelling together with a maximum-likelihood fitting method, in fact, the disentangling of the mass components is accomplished with a very high resolution.