We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), are possible infant GC candidates. Objects belonging to the second group have a high stellar fraction ($f_{\rm star} \geq 0.95$) and show a potential resemblance to infant ultra-faint dwarf galaxies. The high baryon fraction objects identified in this study are characterised by a stellar content similar to the one observed in present-day GCs, but they still contain a high gas fraction ($f_{\rm gas} \sim 0.95$) and a relatively low amount of dark matter. They are compact, dense systems. Their sizes are consistent with recent estimates based on the first observations of possible proto-GCs at high redshifts. These types of infant GC candidates appear to be more massive and more abundant in massive host galaxies, indicating that the assembly of galaxies via mergers may play an important role in building several GC-host scaling relations. Specifically, we express the relation between the mass of the most massive infant GC and its host stellar mass as $\log(M_{\rm cl}) = (0.31\pm0.15)\log(M_{\rm *,gal} + (4.17\pm1.06)$. We also report a new relation between the most massive infant GC and the parent specific star formation rate of the form $\log(M_{\rm cl}) = (0.85\pm0.30)\log(sSFR) + \alpha$ that describes the data at both low and high redshift. Finally, we assess the present-day GC mass (GC number) -- halo mass relation offers a satisfactory description of the behaviour of our infant GC candidates at high redshift, suggesting that such a relation may be set at formation.

Earliest quasars at the cosmic dawn are powered by mass accretion onto supermassive black holes of a billion solar masses. Massive black hole seeds forming through the direct collapse mechanism are considered the most promising candidates but how do they grow and co-evolve with their host galaxies at early cosmic times remains unknown. We here present results from a cosmological radiation hydrodynamical simulation including self-consistent modeling of both Pop III and Pop II star formation, their radiative and supernova feedback in the host galaxy along with X-ray feedback from an accreting massive black hole (MBH) of $\rm 10^5 ~M_{\odot}$ in a halo of $\rm 2 \times 10^9~M_{\odot}$ from $z=26$ down to $z=16$. Our results show that energy deposition from X-rays in the proximity of MBH suppresses Pop III star formation for about 12 Myr while at the same time these X-rays catalyze $\rm H_2$ formation which leads to the formation of a Pop III star cluster of 500 $\rm M_{\odot}$ in the close vicinity of the MBH. We find that mode of star formation for Pop III is episodic and bursty due to the clumpy accretion while for Pop II it is continuous. The stellar mass of the host galaxy at $z \sim 16$ is $\rm 2 \times 10^7~M_{\odot}$ with a star formation rate (SFR) of $\rm \sim 0.1-1~M_{\odot}/yr$. In total, the MBH accretes $\rm 1.5 \times 10^6~M_{\odot}$ during 120 Myr with the mean accretion rate of $\rm \sim 0.01~M_{\odot}/yr$ corresponding to an average Eddington fraction of 50\%.

Supermassive primordial stars forming during catastrophic baryon collapse in atomically-cooling halos at $z \sim$ 15 - 20 may be the origin of the first quasars in the universe. However, no simulation to date has followed the evolution of these halos at resolutions that are high enough or for times that are long enough to determine if collapse actually produces SMSs. Here we report new cosmological simulations of baryon collapse in atomically-cooled halos for times that are long enough for SMSs to form and die as direct-collapse black holes (DCBHs). We find that the high infall rates required to build up such stars do persist until the end of their lives and could fuel the rapid growth of their BHs thereafter. Our simulations also demonstrate that binary and even small multiples of SMSs can form in low-spin and high-spin halos, respectively. This discovery raises the exciting prospect of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the {\em James Webb Space Telescope} and ground-based telescopes in the coming decade.

We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(\epsilon\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions.

Recent observations have successfully detected [O III] $88.3\,{\rm \mu m}$ and [C II] $157.6\,{\rm \mu m}$ lines from galaxies in the early Universe with the Atacama Large Millimeter Array (ALMA). Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present relations between the metal line emission and galaxy evolution at $z=6-15$. We find that galaxies during their starburst phases have high [O III] luminosity of $\sim 10^{42}~\rm erg~s^{-1}$. Once supernova feedback quenches star formation, [O III] luminosities rapidly decrease and continue to be zero for $\sim 100\,{\rm Myr}$. The slope of the relation between $\log{(\rm SFR/M_{\odot}~yr^{-1})}$ and $\log{(L_{\rm [O_{III}]}/{\rm L_{\odot}})}$ at $z=6-9$ is 1.03, and 1.43 for $\log{(L_{\rm [C_{II}]}/{\rm L_{\odot}})}$. As gas metallicity increases from sub-solar to solar metallicity by metal enrichment from star formation and feedback, the line luminosity ratio $L_{\rm [O_{III}]} / L_{\rm [C_{II}]}$ decreases from $\sim 10$ to $\sim 1$ because the O/C abundance ratio decreases due to carbon-rich winds from AGB stars. Therefore, we suggest that the combination of [O III] and [C II] lines is a good probe to investigate the relative distribution of ionized and neutral gas in high-$z$ galaxies. In addition, we show that deep [C II] observations with a sensitivity of $\sim 10^{-2}~{\rm mJy~arcsec^{-2}}$ can probe the extended neutral gas disks of high-$z$ galaxies.

The first galaxies forming a few hundred million years after the Big Bang are the key drivers of cosmic evolution and ideal laboratories to study theories of galaxy formation. We here study the role of UV radiation in suppressing star formation in primordial galaxies by destroying molecular hydrogen, the main coolant in primordial gas and provide estimates of cold dense gas at the onset of star formation. To accomplish this goal, we perform three dimensional cosmological simulations of minihalos in different environments forming at $ z\sim 25$ by varying strength of background UV flux below the Lyman limit between 0.01-1000 in units of $\rm J_{21}=10^{-21}~erg/cm^2/s/Hz/sr$. Particularly, we include photo-detachment of $\rm H^-$, the self-shielding of $\rm H_2$ which both were neglected in previous studies and use updated reaction rates. Our results show that depending on the background level $\rm H_2$ formation is suppressed, delaying gravitational collapse until halos reach the atomic cooling limit. We find that the formation of cold dense molecular gas and subsequently star formation gets delayed by 100 to 230 Myr depending on the level of the background radiation and the growth history of the dark matter halos. The fraction of dense self-shielded gas is a strong function of the background flux and exponentially declines with the strength of incident UV flux above $\rm J_{21} \geq 1$. We find that taking into account $\rm H_2$ self-shielding is crucial for accurately estimating the amount of cold dense gas available for star formation.

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

Super-massive black holes (SMBHs) observed to have masses of $M_\bullet \sim 10^9 \, \mathrm{M_\odot}$ at $z\gtrsim6$, $<1$ Gyr after the Big Bang, are thought to have been seeded by massive black holes which formed before growing concurrently with the formation of their host galaxies. We model analytically the idealised growth of seed black holes, fed through gas inflow from growing proto-galaxy discs. The inflow depends on the disc gravitational stability and thus varies with black hole and disc mass. We find that for a typical host halo, the efficiency of angular momentum transport, as parametrised by the disc viscosity, is the limiting factor in determining the inflow rate and the black hole accretion rate. For our fiducial case we find an upper black hole mass estimate of $M_\bullet \sim 1.8 \times 10^7 \, \mathrm{M_{\odot}}$ at $z=6$. Only in the extreme case of $\sim 10^{16}$ M$_{\odot}$ haloes at $z=6$ produces SMBH masses of $\sim 10^9$ M$_{\odot}$. However, the number density of such haloes is many orders of magnitude below the estimated 1 Gpc$^{-3}$ of SMBHs at $z=6$, indicating that viscosity driven accretion is too inefficient to feed the growth of seeds into $M_\bullet \sim 10^9 \, \mathrm{M_\odot}$ SMBHs by $z \sim 6$. We demonstrate that major mergers are capable of resolving the apparent discrepancy in black hole mass at $z=6$, with some dependence on the exact choice of orbital parameters of the merger.

Massive $10^{4-5}\rm\ M_{\odot}$ black hole seeds resulting from the \textit{direct} collapse of pristine gas require a metal-free atomic cooling halo with extremely low H$_2$ fraction, allowing the gas to cool isothermally in the presence of atomic hydrogen. In order to achieve this chemo-thermodynamical state, the gas needs to be irradiated by both: Lyman-Werner (LW) photons in the energy range $11.2-13.6$ eV capable of photodissociating H$_2$, and $0.76$ eV photons capable of photodetaching H$^-$. Employing cosmological simulations capable of creating the first galaxies in high resolution, we explore if there exists a subset that favours DCBH formation in their vicinity. We find a fundamental relation between the maximum distance at which a galaxy can cause DCBH formation and its star formation rate (SFR), which automatically folds in the chemo-thermodynamical effects of both H$_2$ photo-dissociation and H$^-$ photo-detachment. This is in contrast to the scatter in the LW flux parameter seen at the maximum distance. It shows up to a 3 order of magnitude scatter, which can be interpreted as a scatter in `J$_{crit}$'. Thus, computing the rates and/or the LW flux from a galaxy is no longer necessary to identify neighbouring sites of DCBH formation, as our relation allows one to distinguish regions where DCBH formation could be triggered in the vicinity of a galaxy of a given SFR.

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

The dark ages of the Universe end with the formation of the first generation of stars residing in primeval galaxies. These objects were the first to produce ultraviolet ionizing photons in a period when the cosmic gas changed from a neutral state to an ionized one, known as Epoch of Reionization (EoR). A pivotal aspect to comprehend the EoR is to probe the intertwined relationship between the fraction of ionizing photons capable to escape dark haloes, also known as the escape fraction ($f_{esc}$), and the physical properties of the galaxy. This work develops a sound statistical model suitable to account for such non-linear relationships and the non-Gaussian nature of $f_{esc}$. This model simultaneously estimates the probability that a given primordial galaxy starts the ionizing photon production and estimates the mean level of the $f_{esc}$ once it is triggered. The model was employed in the First Billion Years simulation suite, from which we show that the baryonic fraction and the rate of ionizing photons appear to have a larger impact on $f_{esc}$ than previously thought. A naive univariate analysis of the same problem would suggest smaller effects for these properties and a much larger impact for the specific star formation rate, which is lessened after accounting for other galaxy properties and non-linearities in the statistical model.

Super-massive black holes (SMBHs) with $M_{\bullet} \sim 10^9 $ M$_{\odot}$ at $z>6$ likely originate from massive seed black holes (BHs). We investigate the consequences of seeding SMBHs with direct collapse BHs (DCBHs) ($M_{\bullet}=10^{4-6}\, \mathrm{M}_\odot$) on proto-galactic disc growth. We show that even in the absence of direct feedback effects, the growth of seed BHs reduces the development of gravitational instabilities in host galaxy discs, suppressing star formation and confining stars to a narrow ring in the disc and leading to galaxies at $z \sim 6$ which lie above the local BH-stellar mass relation. The relative magnitude of cosmic and BH accretion rates governs the evolution of the BH-stellar mass relation. For typical DCBH formation epochs, $z_{\rm{i}} \sim 10$, we find star formation is inhibited in haloes growing at the average rate predicted by $\Lambda$CDM which host BHs capable of reaching $M_{\bullet}\sim 10^9 \, \mathrm{M}_{\odot}$ by $z\gtrsim6$. Slower growing BHs cause a delay in the onset of star formation; a $M_{\bullet} \sim 10^6 $ M$_{\odot}$ seed growing at $0.25$ times the Eddington limit will delay star formation by $\sim100$ Myr. This delay is reduced by a factor of $\sim10$ if the halo growth rate is increased by $\sim 0.6\, \sigma$. Our results suggest that SMBHs seeded by DCBHs and their host galaxies form in separate progenitor haloes. In the absence of subsequent mergers, higher than average cosmic accretion or earlier seed formation ($z_{\rm i} \sim 20$) are required to place the evolving BH on the local BH-stellar mass relation by $z=6$.

Recent observations have successfully detected UV-bright and infrared-bright galaxies in the epoch of reionization. However, the origin of their radiative properties has not been understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present predictions of multi-wavelength radiative properties of the first galaxies at $z\sim 6-15$. Using zoom-in initial conditions, we investigate three massive galaxies and their satellites in different environment and halo masses: $M_{\rm h} = 2.4\times 10^{10}{\rm M_{\odot}}$ (Halo-10), $1.6\times 10^{11}{\rm M_{\odot}}$ (Halo-11) and $0.7\times 10^{12}{\rm M_{\odot}}$ (Halo-12) at $z=6$. We find that most of gas and dust are ejected from star-forming regions by supernova feedback, which allows UV photons to escape. We show that the peak of the spectral energy distribution (SED) rapidly changes between UV and infrared wavelengths on a time-scale of $\sim$\,100\,Myrs due to intermittent star formation and feedback, and the escape fraction of UV photons fluctuates in the range of $0.2-0.8$ at $z<10$ with a time-averaged value of 0.3. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. We predict the detectability of high-$z$ galaxies with the Atacama Large Millimeter Array (ALMA). For a sensitivity limit of $0.1\,{\rm mJy}$ at $850\,{\rm \mu m}$, the detection probability of galaxies in halos $M_{\rm h} \gtrsim 10^{11}\,{\rm M_{\odot}}$ at $z\lesssim 7$ exceeds fifty per cent. We argue that supernova feedback can produce the observed diversity of SEDs for high-$z$ galaxies.

Massive Ultracompact Galaxies (MUGs) are common at z=2-3, but very rare in the nearby Universe. Simulations predict that the few surviving MUGs should reside in galaxy clusters, whose large relative velocities prevent them from merging, thus maintaining their original properties (namely stellar populations, masses, sizes and dynamical state). We take advantage of the high-completeness, large-area spectroscopic GAMA survey, complementing it with deeper imaging from the KiDS and VIKING surveys. We find a set of 22 bona-fide MUGs, defined as having high stellar mass (>8x10^10 M_Sun) and compact size (R_e<2 Kpc) at 0.02 < z < 0.3. An additional set of 7 lower-mass objects (6x10^10 < M_star/M_Sun < 8x10^10) are also potential candidates according to typical mass uncertainties. The comoving number density of MUGs at low redshift (z < 0.3) is constrained at $(1.0\pm 0.4)x 10^-6 Mpc^-3, consistent with galaxy evolution models. However, we find a mixed distribution of old and young galaxies, with a quarter of the sample representing (old) relics. MUGs have a predominantly early/swollen disk morphology (Sersic index 1<n<2.5) with high stellar surface densities (<Sigma_e> ~ 10^10 M_Sun Kpc^-2). Interestingly, a large fraction feature close companions -- at least in projection -- suggesting that many (but not all) reside in the central regions of groups. Halo masses show these galaxies inhabit average-mass groups. As MUGs are found to be almost equally distributed among environments of different masses, their relative fraction is higher in more massive overdensities, matching the expectations that some of these galaxies fell in these regions at early times. However, there must be another channel leading some of these galaxies to an abnormally low merger history because our sample shows a number of objects that do not inhabit particularly dense environments. (abridged)

We present mass estimates of supermassive black holes in six nearby fast rotating early-type galaxies (NGC4339, NGC4434, NGC4474, NGC4551, NGC4578 and NGC4762) with effective stellar velocity dispersion around 100 km/s. We use near-infrared laser-guide adaptive optics observations with the GEMINI/NIFS to derive stellar kinematics in the galactic nuclei, and SAURON observations from the ATLAS3D Survey for large-scale kinematics. We build axisymmetric Jeans anisotropic models and axisymmetric Schwarzschild dynamical models. Both modelling approaches recover consistent orbital anisotropies and black hole masses within 1-2sigma confidence level, except for one galaxy for which the difference is just above the 3sigma level. Two black holes (NGC4339 and NGC4434) are amongst the largest outliers from the current black hole mass - velocity dispersion relation, with masses of $(4.3^{+4.8}_{-2.3})\times10^7$ and $(7.0^{+2.0}_{-2.8})\times10^7$ M$_\odot$, respectively ($3\sigma$ confidence level). The black holes in NGC4578 and NGC4762 lie on the scaling relation with masses of $(1.9^{+0.6}_{-1.4})\times10^7$ and $(2.3^{+0.9}_{-0.6})\times10^7$ M$_\odot$, respectively (3sigma confidence level). For two galaxies (NGC4474 and NGC4551) we are able to place upper limits on their black holes masses ($<7\times10^6$ and $<5\times10^6$ M$_\odot$, respectively, $3\sigma$ confidence level). The kinematics for these galaxies clearly indicate central velocity dispersion drops within a radius of 35 pc and 80 pc, respectively. These drops cannot be associated with cold stellar structures and our data do not have the resolution to exclude black holes with masses an order of magnitude smaller than the predictions. Parametrizing the orbital distribution in spherical coordinates, the vicinity of the black holes is characterized by isotropic or mildly tangential anisotropy.

We present the results of a new study of the relationship between infrared excess (IRX), UV spectral slope (beta) and stellar mass at redshifts 2<z<3, based on a deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field (HUDF). Excluding the most heavily-obscured sources, we use a stacking analysis to show that z~2.5 star-forming galaxies in the mass range 9.25 <= log(M/Msun) <= 10.75 are fully consistent with the IRX-beta relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass complete, sample of 2 <= z <= 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between beta and stellar mass, making it possible to predict UV attenuation (A_1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z~2.5 star-forming galaxies follow A_1600-mass and IRX-mass relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX-beta relation may have been biased toward low values of IRX at red values of beta, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than beta. Although the situation at lower stellar masses remains uncertain, we conclude that for 2<z<3 star-forming galaxies with log(M/Msun) >= 9.75, both the IRX-beta and IRX-mass relations are well described by a Calzetti-like attenuation law.

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.

The diffuse soft X-ray emissivity from galactic winds is computed during the Epoch of Reionization (EoR). We consider two analytic models, a pressure-driven wind and a superbubble model, and a 3D cosmological simulation including gas dynamics from the First Billion Years (FiBY) project. The analytic models are normalized to match the diffuse X-ray emissivity of star-forming galaxies in the nearby Universe. The cosmological simulation uses physically motivated star formation and wind prescriptions, and includes radiative transfer corrections. The models and the simulation all are found to produce sufficient heating of the Intergalactic Medium to be detectable by current and planned radio facilities through 21 cm measurements during the EoR. While the analytic models predict a 21 cm emission signal relative to the Cosmic Microwave Background sets in by $z_{\rm trans} \simeq 8 - 10$, the predicted signal in the FiBY simulation remains in absorption until reionization completes. The 21 cm absorption differential brightness temperature reaches a minimum of $\Delta T \simeq -130$ to $-200$ mK, depending on model. Allowing for additional heat from high mass X-ray binaries pushes the transition to emission to $z_{\rm trans} \simeq 10 - 12$, with shallower absorption signatures having a minimum of $\Delta T \simeq -110$ to $-140$ mK. The 21 cm signal may be a means of distinguishing between the wind models, with the superbubble model favouring earlier reheating. While an early transition to emission may indicate X-ray binaries dominate the reheating, a transition to emission as early as $z_{\rm trans} > 12$ would suggest the presence of additional heat sources.

We present the results of a study investigating the dust attenuation law at $z\simeq 5$, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N=498 galaxies drawn from the First Billion Years (FiBY) simulation project. The simulated galaxies at $z\simeq 5$, which have M$_{1500} \leq -18.0$ and $7.5 \leq \rm{log(M/M}_{\odot}\rm{)} \leq 10.2$, display a mass-dependent $\alpha$-enhancement, with a median value of $[\alpha/\rm{Fe}]_{z=5}~\simeq~4~\times~[\alpha/\rm{Fe}]_{Z_{\odot}}$. The median Fe/H ratio of the simulated galaxies is $0.14\pm0.05$ which, even including the effects of nebular continuum, produces steep intrinsic UV continuum slopes; $\langle \beta_{i} \rangle = -2.4 \pm 0.05$. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form $A(\lambda) \propto \lambda^{n}$, we explore the parameter values which best reproduce the observed $z=5$ luminosity function (LF) and colour-magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass, and the dust attenuation slope ($n$) is within the range $-0.7 \leq n \leq-0.3$, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model (SPS) assumptions. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV ($n=-0.55$), and other similarly `grey' star-forming galaxy attenuation curves recently derived at $z\simeq2$. In contrast, attenuation curves as steep as the Small Magellanic Cloud (SMC) extinction curve ($n=-1.24$) are formally ruled out. Finally, we show that our models are consistent with recent 1.3mm ALMA observations of the Hubble Ultra Deep Field (HUDF), and predict the form of the $z\simeq5$ IRX$-\beta$ relation.

In this study we focus on the $z\sim6.6$ Lyman-$\alpha$ CR7 consisting of clump A that is host to a potential direct collapse black hole (DCBH), and two metal enriched star forming clumps B and C. In contrast to claims that signatures of metals rule out the existence of DCBHs, we show that metal pollution of A from star forming clumps clumps B and C is inevitable, and that A can form a DCBH well before its metallicity exceeds the critical threshold of $10^{-5}-10^{-6}\ \rm Z_{\odot}$. Assuming metal mixing happens instantaneously, we derive the metallicity of A based on the star formation history of B and C. We find that treating a final accreting black hole of $10^6-10^7\ \rm M_{\odot}$ in A for nebular emission already pushes its $H_{160}$ - [3.6] and [3.6]-[4.5] colours into the 3$\sigma$ limit of observations. Hence, we show that the presence of metals in DCBH hosts is inevitable, and that it is the coevolution of the LW radiation field and metals originating from neighbouring galaxies that governs DCBH formation in a neighbouring {initially} pristine atomic cooling haloes.

The critical Lyman--Werner flux required for direct collapse blackholes (DCBH) formation, or J$_{crit}$, depends on the shape of the irradiating spectral energy distribution (SED). The SEDs employed thus far have been representative of {{realistic}} single stellar populations. We study the effect of binary stellar populations on the formation of DCBH, as a result of their contribution to the Lyman--Werner radiation field. Although binary populations with ages $>$ 10 Myr yield a larger LW photon output, we find that the corresponding values of J$_{crit}$ can be up to 100 times higher than single stellar populations. We attribute this to the shape of the binary SEDs as they produce a sub--critical rate of H$^-$ photodetaching 0.76 eV photons as compared to single stellar populations, reaffirming the role that H$^-$ plays in DCBH formation. This further corroborates the idea that DCBH formation is better understood in terms of a critical region in the H$_2$--H$^-$ photo--destruction rate parameter space, rather than a single value of LW flux.

We investigate the environment in which direct-collapse black holes may form by analysing a cosmological, hydrodynamical simulation that is part of the First Billion Years project. This simulation includes the most relevant physical processes leading to direct collapse of haloes, most importantly, molecular hydrogen depletion by dissociation of $H_2$ and $H^-$ from the evolving Lyman-Werner radiation field. We selected a sample of pristine atomic cooling haloes that have never formed stars in their past, have not been polluted with heavy elements and are cooling predominantly via atomic hydrogen lines. Amongst them we identified six haloes that could potentially harbour massive seed black holes formed via direct collapse (with masses in the range of $10^{4-6} M_{sun}$). These potential hosts of direct-collapse black holes form as satellites and are found within 15 physical kpc of proto-galaxies, with stellar masses in the range $10^{5-7} M_{sun}$ and maximal star formation rates of 0.1 Msun/yr over the past 5 Myr, and are exposed to the highest flux of Lyman-Werner radiation emitted from the neighbouring galaxies. It is the proximity to these proto-galaxies that differentiates these haloes from rest of the sample.

We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$\alpha$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We here derive upper limits on the mass and luminosity of Population III (POPIII) dominated proto-galaxies based on the collapse of primordial gas under the effect of angular momentum loss via Ly$\alpha$ radiation drag and the gas accretion onto a galactic centre. Our model predicts that POPIII-dominated galaxies at z ~ 7 are hosted by haloes with $M_{\rm halo} \sim 1.5 \times 10^{8} - 1.1 \times 10^{9} \rm ~M_{\odot}$, that they have Ly$\alpha$ luminosities of $L_{\rm Ly\alpha} \sim 3.0 \times 10^{42} - 2.1 \times 10^{43}$ erg/s, stellar mass of $M_{\rm star} \sim 0.8 \times 10^{5} - 2.5 \times 10^{6} \rm ~M_{\odot}$, and outflowing gas with velocities $V_{\rm out} \sim 40$ km/s due to Ly$\alpha$ radiation pressure. We show that the POPIII galaxy candidate CR7 violates the derived limits on stellar mass and Ly$\alpha$ luminosity and thus is unlikely to be a POPIII galaxy. POPIII-dominated galaxies at z ~ 7 have HeII line emission that is ~1- 3 orders of magnitude lower then that of Ly$\alpha$, they have high Ly$\alpha$ equivalent width of > 300 $\AA$ and should be found close to bright star forming galaxies. The HeII 1640 $\AA$ line is in comfortable reach of next generation telescopes, like the JWST or TMT.

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $\Lambda$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of assembling supermassive black holes by $z>6$ to explain the bright quasars observed at this epoch. The formation of a DCBH seed with $\rm M_{BH}\sim10^{4-5}\ \rm M_{\odot}$ requires a pristine atomic-cooling halo to be illuminated by an external radiation field that is sufficiently strong to entirely suppress H$_{2}$ cooling in the halo. Many previous studies have attempted to constrain the critical specific intensity that is likely required to suppress H$_{2}$ cooling, denoted as $J_{\rm crit}$. However, these studies have typically assumed that the incident external radiation field can be modeled with a black-body spectrum. Under this assumption, it is possible to derive a {unique} value for $J_{\rm crit}$ that depends only on the temperature of the black-body. In this study we consider a more realistic spectral energy distribution (SED) for the external source of radiation that depends entirely on its star formation history and age. The rate of destruction of the species responsible for suppressing molecular hydrogen cooling depends on the detailed shape of the SED. Therefore the value of $J_{\rm crit}$ is tied to the shape of the incident SED of the neighbouring galaxy. We fit a parametric form to the rates of destruction of H$_2$ and H$^-$ that permit direct collapse. Owing to this, we find that $J_{\rm crit}$ is not a fixed threshold but can lie anywhere in the range $J_{\rm crit} \sim 0.5$--$10^{3}$, depending on the details of the source stellar population.

We present the results of a high-resolution, 5 GHz, Karl G. Jansky Very Large Array study of the nuclear radio emission in a representative subset of the Atlas3D survey of early-type galaxies (ETGs). We find that 51 +/- 4% of the ETGs in our sample contain nuclear radio emission with luminosities as low as 10^18 W/Hz. Most of the nuclear radio sources have compact (< 25-110 pc) morphologies, although < 10% display multi-component core+jet or extended jet/lobe structures. Based on the radio continuum properties, as well as optical emission line diagnostics and the nuclear X-ray properties, we conclude that the majority of the central 5 GHz sources detected in the Atlas3D galaxies are associated with the presence of an active galactic nucleus (AGN). However, even at sub-arcsecond spatial resolution, the nuclear radio emission in some cases appears to arise from low-level nuclear star formation rather than an AGN, particularly when molecular gas and a young central stellar population is present. This is in contrast to popular assumptions in the literature that the presence of a compact, unresolved, nuclear radio continuum source universally signifies the presence of an AGN. Additionally, we examine the relationships between the 5 GHz luminosity and various galaxy properties including the molecular gas mass and - for the first time - the global kinematic state. We discuss implications for the growth, triggering, and fueling of radio AGNs, as well as AGN-driven feedback in the continued evolution of nearby ETGs.

We present calculations on the formation of massive black holes with 10^5 Msun at z > 6 that can be the seeds of supermassive black holes at z > 6. Under the assumption of compact star cluster formation in merging galaxies, star clusters in haloes of 10^8 ~ 10^9 Msun can undergo rapid core-collapse leading to the formation of very massive stars (VMSs) with ~1000 Msun which directly collapse into black holes with similar masses. Star clusters in halos of > 10^9 Msun experience type-II supernovae before the formation of VMSs due to long core-collapse time scales. We also model the subsequent growth of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small radii is the main driver for accretion of stars onto black holes. As a result, more than ninety percent of stars in the initial cluster are swallowed by the central black holes before z=6. Using dark matter merger trees we derive black hole mass functions at z=6-20. The mass function ranges from 10^3 to 10^5 Msun at z <~ 15. Major merging of galaxies of >~ 4*10^8 Msun at z ~ 20 successfully leads to the formation of >~ 10^5 Msun BHs by z >~ 10 which can be the potential seeds of supermassive black holes seen today.

Approximately 20% of early-type galaxies host small nuclear stellar discs that are tens to a few hundred parsecs in size. Such discs are expected to be easily disrupted during major galactic encounters, hence their age serve to constrain their assembly history. We use VIMOS integral-field spectroscopic observations for the intermediate-mass E0 galaxy NGC4458 and age-date its nuclear disc via high-resolution fitting of various model spectra. We find that the nuclear disc is at least 6 Gyr old. A clue to gain narrow limits to the stellar age is our knowledge of the nuclear disc contribution to the central surface brightness. The presence of an old nuclear disk, or the absence of disruptive encounters since z~0.6, for a small galaxy such as NGC4458 which belongs to the Virgo cluster, may be consistent with a hierarchical picture for galaxy formation where the smallest galaxies assembles earlier and the crowded galactic environments reduce the incidence of galaxy mergers. On the other hand, NGC4458 displays little or no bulk rotation except for a central kpc-scale kinematically-decoupled core. Slow rotation and decoupled core are usually explained in terms of mergers. The presence and age of the nuclear disc constraint these mergers to have happened at high redshift.

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.

We present the stellar population content of early-type galaxies from the Atlas3D survey. Using spectra integrated within apertures covering up to one effective radius, we apply two methods: one based on measuring line-strength indices and applying single stellar population (SSP) models to derive SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and one based on spectral fitting to derive non-parametric star-formation histories, mass-weighted average values of age, metallicity, and half-mass formation timescales. Using homogeneously derived effective radii and dynamically-determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_JAM, Sigma_e, R_maj), showing that at fixed mass, compact early-type galaxies are on average older, more metal-rich, and more alpha-enhanced than their larger counterparts. From non-parametric star-formation histories, we find that the duration of star formation is systematically more extended in lower mass objects. Assuming that our sample represents most of the stellar content of today's local Universe, approximately 50% of all stars formed within the first 2 Gyr following the big bang. Most of these stars reside today in the most massive galaxies (>10^10.5 M_sun), which themselves formed 90% of their stars by z~2. The lower-mass objects, in contrast, have formed barely half their stars in this time interval. Stellar population properties are independent of environment over two orders of magnitude in local density, varying only with galaxy mass. In the highest-density regions of our volume (dominated by the Virgo cluster), galaxies are older, alpha-enhanced and have shorter star-formation histories with respect to lower density regions.

Using varying models for the star formation rate (SFR) of Population (Pop) III and II stars at z>6 we derive the expected redshift history of the global 21 cm signal from the inter-galactic medium (IGM). To recover the observed Thomson scattering optical depth of the cosmic microwave background (CMB) requires SFRs at the level of ~ 10^-3 Msun/yr/Mpc^3 at z ~ 15 from Pop III stars, or ~ 10^-1Msun/yr/Mpc^3 at z ~ 7 from Pop II stars. In the case the SFR is dominated by Pop III stars, the IGM quickly heats above the CMB at z > 12 due to heating from supernovae. In addition, Lyman-alpha photons from haloes hosting Pop III stars couple the spin temperature to that of the gas, resulting in a deep absorption signal. If the SFR is dominated by Pop II stars, the IGM slowly heats and exceeds the CMB temperature at z ~ 10. However, the larger and varying fraction of Pop III stars is able to break this degeneracy. We find that the impact of the initial mass function (IMF) of Pop III stars on the 21 cm signal results in an earlier change to a positive signal if the IMF slope is ~ -1.2. Measuring the 21 cm signal at z > 10 with next generation radio telescopes such as the Square Kilometre Array will be able to investigate the contribution from Pop III and Pop II stars to the global star formation rate.

Long gamma-ray burst's (LGRB's) association to the death of massive stars suggest they could be used to probe the cosmic star formation history (CSFH) with high accuracy, due to their high luminosities. We utilise cosmological simulations from the First Billion Years project to investigate the biases between the CSFH and the LGRB rate at z>5, assuming various different models and constraints on the progenitors of LGRBs. We populate LGRBs using a selection based on environmental properties and demonstrate that the LGRB rate should trace the CSFH to high redshifts. The measured LGRB rate suggests that LGRBs have opening angles of theta_jet=0.1 deg, although the degeneracy with the progenitor model cannot rule out an underlying bias. We demonstrate that proxies that relate the LGRB rate with global LGRB host properties do not reflect the underlying LGRB environment, and are in fact a result of the host galaxy's spatial properties, such that LGRBs can exist in galaxies of solar metallicity. However, we find a class of host galaxies that have low stellar mass and are metal-rich, and that their metallicity dispersions would not allow low-metallicity environments. Detection of hosts with this set of properties would directly reflect the progenitor's environment. We predict that 10% of LGRBs per year are associated with this set of galaxies that would have forbidden line emission that could be detected by instruments on the James Webb Space Telescope. Such a discovery would place strong constraints on the collapsar model and suggest other avenues to be investigated.

Galactic archeology based on star counts is instrumental to reconstruct the past mass assembly of Local Group galaxies. The development of new observing techniques and data-reduction, coupled with the use of sensitive large field of view cameras, now allows us to pursue this technique in more distant galaxies exploiting their diffuse low surface brightness (LSB) light. As part of the Atlas3D project, we have obtained with the MegaCam camera at the Canada-France Hawaii Telescope extremely deep, multi--band, images of nearby early-type galaxies. We present here a catalog of 92 galaxies from the Atlas3D sample, that are located in low to medium density environments. The observing strategy and data reduction pipeline, that achieve a gain of several magnitudes in the limiting surface brightness with respect to classical imaging surveys, are presented. The size and depth of the survey is compared to other recent deep imaging projects. The paper highlights the capability of LSB--optimized surveys at detecting new prominent structures that change the apparent morphology of galaxies. The intrinsic limitations of deep imaging observations are also discussed, among those, the contamination of the stellar halos of galaxies by extended ghost reflections, and the cirrus emission from Galactic dust. The detection and systematic census of fine structures that trace the present and past mass assembly of ETGs is one of the prime goals of the project. We provide specific examples of each type of observed structures -- tidal tails, stellar streams and shells --, and explain how they were identified and classified. We give an overview of the initial results. The detailed statistical analysis will be presented in future papers.

We report on empirical trends between the dynamically determined stellar initial mass function (IMF) and stellar population properties for a complete, volume-limited sample of 260 early-type galaxies from the Atlas3D project. We study trends between our dynamically-derived IMF normalisation and absorption line strengths, and interpret these via single stellar population- (SSP-) equivalent ages, abundance ratios (measured as [alpha/Fe]), and total metallicity, [Z/H]. We find that old and alpha-enhanced galaxies tend to have on average heavier (Salpeter-like) mass normalisation of the IMF, but stellar population does not appear to be a good predictor of the IMF, with a large range of normalisation at a given population parameter. As a result, we find weak IMF-[alpha/Fe] and IMF-age correlations, and no significant IMF-[Z/H] correlation. The observed trends appear significantly weaker than those reported in studies that measure the IMF normalisation via low-mass star demographics inferred through stellar spectral analysis.

We use the Atlas3D sample to perform a study of the intrinsic shapes of early-type galaxies, taking advantage of the available combined photometric and kinematic data. Based on our ellipticity measurements from the Sloan Digital Sky Survey Data Release 7, and additional imaging from the Isaac Newton Telescope, we first invert the shape distribution of fast and slow rotators under the assumption of axisymmetry. The so-obtained intrinsic shape distribution for the fast rotators can be described with a Gaussian with a mean flattening of q=0.25 and standard deviation sigma_q = 0.14, and an additional tail towards rounder shapes. The slow rotators are much rounder, and are well described with a Gaussian with mean q = 0.63 and sigma_q =0.09. We then checked that our results were consistent when applying a different and independent method to obtain intrinsic shape distributions, by fitting the observed ellipticity distributions directly using Gaussian parametrisations for the intrinsic axis ratios. Although both fast and slow rotators are identified as early-type galaxies in morphological studies, and in many previous shape studies are therefore grouped together, their shape distributions are significantly different, hinting at different formation scenarios. The intrinsic shape distribution of the fast rotators shows similarities with the spiral galaxy population. Including the observed kinematic misalignment in our intrinsic shape study shows that the fast rotators are predominantly axisymmetric, with only very little room for triaxiality. For the slow rotators though there are very strong indications that they are (mildly) triaxial.

We study the effect of baryons on the inner dark matter profile of the first galaxies using the First Billion Years simulation between z=16-6 before secular evolution sets in. Using a large statistical sample from two simulations of the same volume and cosmological initial conditions, one with and one without baryons, we are able to directly compare haloes with their baryon-free counterparts, allowing a detailed study of the modifications to the dark matter density profile due to the presence of baryons during the first billion years of galaxy formation. For each of the ~ 5000 haloes in our sample we quantify the impact of the baryons using eta, defined as the ratio of dark matter mass enclosed in 100 pc in the baryonic run to its counterpart without baryons. During this epoch of rapid growth of galaxies, we find that many haloes of these first galaxies show an enhancement of dark matter in the halo centre compared to the baryon-free simulation, while many others show a deficit. We find that the mean value of eta is close to unity, but there is a large dispersion, with a standard deviation of 0.677. The enhancement is cyclical in time and tracks the star formation cycle of the galaxy; as gas falls to the centre and forms stars, the dark matter moves in as well. Supernova feedback then removes the gas, and the dark matter again responds to the changing potential. We study three physical models relating the motion of baryons to that of the dark matter: adiabatic contraction, dynamical friction, and rapid outflows. Abridged, see text for full abstract

Understanding the processes that can destroy H$_2$ and H$^-$ species is quintessential in governing the formation of the first stars, black holes and galaxies. In this study we compute the reaction rate coefficients for H$_2$ photo--dissociation by Lyman--Werner photons ($11.2 - 13.6$ eV), and H$^-$ photo--detachment by 0.76 eV photons emanating from self-consistent stellar populations that we model using publicly available stellar synthesis codes. So far studies that include chemical networks for the formation of molecular hydrogen take these processes into account by assuming that the source spectra can be approximated by a power-law dependency or a black-body spectrum at 10$^4$ or $10^5$ K. We show that using spectra generated from realistic stellar population models can alter the reaction rates for photo-dissociation, $\rm k_{\rm{di}}$, and photo-detachment, $\rm k_{\rm{de}}$, significantly. In particular, $\rm k_{\rm{de}}$ can be up to $\sim 2-4$ orders of magnitude lower in the case of realistic stellar spectra suggesting that previous calculations have over-estimated the impact that radiation has on lowering H$_2$ abundances. In contrast to burst modes of star formation, we find that models with continuous star formation predict increasing $\rm k_{\rm{de}}$ and $\rm k_{\rm{di}}$, which makes it necessary to include the star formation history of sources to derive self-consistent reaction rates, and that it is not enough to just calculate J$_{21}$ for the background. For models with constant star formation rate the change in shape of the spectral energy distribution leads to a non-negligible late-time contribution to $\rm k_{\rm{de}}$ and $\rm k_{\rm{di}}$, and we present self-consistently derived cosmological reaction rates based on star formation rates consistent with observations of the high redshift Universe.

We constrain the assembly history of galaxies in the projected central 0.5 Mpc of the Coma cluster by performing structural decomposition on 69 massive (M_star >= 10^9 M_sun) galaxies using high-resolution F814W images from the HST Treasury Survey of Coma. Each galaxy is modeled with up to three Sersic components having a free Sersic index n. After excluding the two cDs in the projected central 0.5 Mpc of Coma, 57% of the galactic stellar mass in the projected central 0.5 Mpc of Coma resides in classical bulges/ellipticals while 43% resides in cold disk-dominated structures. Most of the stellar mass in Coma may have been assembled through major (and possibly minor) mergers. Hubble types are assigned based on the decompositions, and we find a strong morphology-density relation; the ratio of (E+S0):spirals is (91.0%):9.0%. In agreement with earlier work, the size of outer disks in Coma S0s/spirals is smaller compared with lower-density environments captured with SDSS (Data Release 2). Among similar-mass clusters from a hierarchical semi-analytic model, no single cluster can simultaneously match all the global properties of the Coma cluster. The model strongly overpredicts the mass of cold gas and underpredicts the mean fraction of stellar mass locked in hot components over a wide range of galaxy masses. We suggest that these disagreements with the model result from missing cluster physics (e.g., ram-pressure stripping), and certain bulge assembly modes (e.g., mergers of clumps). Overall, our study of Coma underscores that galaxy evolution is not solely a function of stellar mass, but also of environment.

We present measurements of the star formation rate (SFR) in the early-type galaxies (ETGs) of the ATLAS3D sample, based on Wide-field Infrared Survey Explorer (WISE) 22um and Galaxy Evolution Explorer (GALEX) far-ultraviolet emission. We combine these with gas masses estimated from 12CO and HI data in order to investigate the star formation efficiency (SFE) in a larger sample of ETGs than previously available. We first recalibrate (based on WISE data) the relation between old stellar populations (traced at Ks-band) and 22um luminosity, allowing us to remove the contribution of 22um emission from circumstellar dust. We then go on to investigate the position of ETGs on the Kennicutt-Schmidt (KS) relation. Molecular gas-rich ETGs have comparable star formation surface densities to normal spiral galaxy centres, but they lie systematically offset from the KS relation, having lower star formation efficiencies by a factor of ~2.5 (in agreement with other authors). This effect is driven by galaxies where a substantial fraction of the molecular material is in the rising part of the rotation curve, and shear is high. We show here for the first time that although the number of stars formed per unit gas mass per unit time is lower in ETGs, it seems that the amount of stars formed per free-fall time is approximately constant. The scatter around this dynamical relation still correlates with galaxy properties such as the shape of the potential in the inner regions. This leads us to suggest that dynamical properties (such as shear or the global stability of the gas) may be important second parameters that regulate star formation and cause much of the scatter around star-formation relations.

We present results on the radiation drag exerted by an isotropic and homogenous background of Lya photons on neutral gas clouds orbiting within HII regions around Population III stars of different masses. The Doppler shift causes a frequency difference between photons moving in the direction of the cloud and opposite to it resulting in a net momentum loss of the cloud in the direction of motion. We find that half of the angular momentum of gas with v_theta <~ 20 km/s near (r <~ 3 kpc) a Population III star of 120 Msun at z=20 is lost within ~ 10^6 yr. The radiation drag is a strong function of cloud velocity that peaks at v ~ 20 km/s reflecting the frequency dependence of the photon cross section. Clouds moving with velocities larger than ~ 100 km/s loose their angular momentum on time scales of ~ 10^8 yr. At lower redshifts radiation drag becomes inefficient as the Lya photon density in HII regions decreases by a factor (1+z)^3 and angular momentum is lost on time scales > 10^8 yr even for low velocity clouds. Our results suggest that a sweet spot exists for the loss of angular momentum by radiation drag for gas clouds at z > 10 and with v ~ 20 km/s. Comparison to dynamical friction forces acting on typical gas clouds suggest that radiation drag is the dominant effect impacting the orbit. We propose that this effect can suppress the formation of extended gas discs in the first galaxies and help gas accretion near galactic centres and central black holes.

One quarter of all nearby early-type galaxies (ETGs) outside Virgo host a disc/ring of HI with size from a few to tens of kpc and mass up to ~1e+9 solar masses. Here we investigate whether this HI is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of HI masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much HI and have a similar rate of HI discs/rings as FRs. Accretion of HI is therefore not always linked to the growth of an inner stellar disc. The weak relation between HI and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the HI is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that LCDM hydrodynamical simulations have difficulties reproducing the HI properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating HI. This establishes the HI properties of ETGs as a novel challenge to simulations.

We present a study of the cold gas contents of the Atlas3D early-type galaxies, in the context of their optical colours, near-UV colours, and H\beta\ absorption line strengths. Early-type (elliptical and lenticular) galaxies are not as gas-poor as previously thought, and at least 40% of local early-type galaxies are now known to contain molecular and/or atomic gas. This cold gas offers the opportunity to study recent galaxy evolution through the processes of cold gas acquisition, consumption (star formation), and removal. Molecular and atomic gas detection rates range from 10% to 34% in red sequence early-type galaxies, depending on how the red sequence is defined, and from 50% to 70% in blue early-type galaxies. Notably, massive red sequence early-type galaxies (stellar masses > 5e10 Msun, derived from dynamical models) are found to have HI masses up to M(HI)/Mstar ~ 0.06 and H_2 masses up to M(H$_2$)/Mstar ~ 0.01. Some 20% of all massive early-type galaxies may have retained atomic and/or molecular gas through their transition to the red sequence. However, kinematic and metallicity signatures of external gas accretion (either from satellite galaxies or the intergalactic medium) are also common, particularly at stellar masses <= 5e10 Msun, where such signatures are found in ~ 50% of H$_2$-rich early-type galaxies. Our data are thus consistent with a scenario in which fast rotator early-type galaxies are quenched former spiral galaxies which have undergone some bulge growth processes, and in addition, some of them also experience cold gas accretion which can initiate a period of modest star formation activity. We discuss implications for the interpretation of colour-magnitude diagrams.

We present new Spectrographic Areal Unit for Research on Optical Nebulae (SAURON) integral-field spectroscopy and Swift Ultraviolet Optical Telescope (UVOT) observations of molecular outflow host galaxy NGC 1266 that indicate NGC 1266 has experienced a rapid cessation of star formation. Both the SAURON maps of stellar population age and the Swift UVOT observations demonstrate the presence of young ($< 1$ Gyr) stellar populations within the central 1 kpc, while existing Combined Array for Research in Millimeter-wave Astronomy (CARMA) CO(1--0) maps indicate that the sites of current star formation are constrained to the inner few hundred parsecs of the galaxy only. The optical spectrum of NGC 1266 from Moustakas & Kennicutt (2006) reveal a characteristic post-starburst (K+A) stellar population and Davis et al. (2012) confirm that ionized gas emission in the system originate from a shock. Galaxies with K+A spectra and shock-like ionized gas line ratios may comprise an important, overlooked segment of the post-starburst population, containing exactly those objects in which the AGN is actively expelling the star-forming material. While AGN activity is not the likely driver of the post-starburst event that occurred 500 Myr ago, the faint spiral structure seen in the Hubble Space Telescope (HST) Wide-field Camera 3 (WFC3) Y-, J- and H-band imaging seems to point to the possibility of gravitational torques being the culprit. If the molecular gas were driven into the center at the same time as the larger scale galaxy disk underwent quenching, the AGN might be able to sustain the presence of molecular gas for $\gtrsim 1$ Gyr by cyclically injecting turbulent kinetic energy into the dense molecular gas via a radio jet, inhibiting star formation.

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies and their satellites. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments $h_3$ and $h_4$ are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the $\lambda_{\mathrm{R}}$-parameter. The velocity, velocity dispersion, $h_3$, and $h_4$ fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS$^{\rm{3D}}$ survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of $h_3$ and $v/\sigma$. An additional formation path for fast rotators includes gas poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated $h_3$ and $v/\sigma$. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since $z \sim 2$, rotate slower and have older stellar populations. (shortened)

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

[Abridged and Edited] We investigate nuclear light profiles in 135 ATLAS3D galaxies for which the Hubble Space Telescope (HST) imaging is available and compare them to the large scale kinematics obtained with the SAURON integral-field spectrograph. Specific angular momentum, lambda_R, correlates with the shape of nuclear light profiles, where cores are typically found in slow rotators and core-less galaxies are fast rotators. Cores are found only in massive galaxies and only in systems with the stellar mass M>8x10^10 Msun. Based on our sample, we, however, see no evidence for a bimodal distribution of nuclear slopes. The best predictor for finding a core is based on the stellar velocity dispersion within an effective radius, sigma_e, and specific angular momentum, where cores are found for lambda_R<0.25 and sigma_e>160 km/s. We estimate that only about 10% of nearby early-type galaxies contain cores. Furthermore, we show that there is a genuine population of fast rotators with cores. We also show that core fast rotators are morphologically, kinematically and dynamically different from core slow rotators. The cores of fast rotators could harbour black holes of similar masses to those in core slow rotators, but typically more massive than those found in core-less fast rotators. Core-less galaxies, and especially core-less fast rotators, are under-luminous in the diffuse X-ray emission, but the presence of a core does not imply high X-ray luminosities. We postulate a possible population of core-less galaxies among slow rotators, which can not be explained as face-on discs, but comprise a genuine sub-population of slow rotators. These galaxies are typically less massive and flatter than core slow rotators, and show evidence for dynamical cold structures and exponential photometric components. We discuss possible processes for formation of cores and their subsequent preservation.

We predict the existence and observational signatures of a new class of objects that assembled early, during the first billion years of cosmic time: Obese Black-hole Galaxies (OBGs). OBGs are objects in which the mass of the central black hole initially exceeds that of the stellar component of the host galaxy, and the luminosity from black-hole accretion dominates the starlight. From a cosmological simulation, we demonstrate that there are sites where star formation is initially inhibited and direct-collapse black holes (DCBHs) form due to the photo-dissociating effect of Lyman-Werner radiation on molecular hydrogen. We show that the formation of OBGs is inevitable, because the probability of finding the required extra-galactic environment and the right physical conditions in a halo conducive to DCBH formation is quite high in the early universe. We estimate an OBG number density of 0.009/Mpc^3 at z~8 and 0.03/Mpc^3 at z~6. Extrapolating from our simulation volume, we infer that the most luminous quasars detected at z~6 likely transited through an earlier OBG phase. We find that these primordial galaxies start off with an over-massive BH and acquire their stellar component from subsequent merging as well as in-situ star formation. In doing so, they inevitably go through an OBG phase dominated by the accretion luminosity at the Eddington rate or below, released from the growing BH. The OBG phase is characterised by an ultra-violet (UV) spectrum with slope of beta ~ -2.3 and the absence of a Balmer Break. OBGs should also be spatially unresolved, and are expected to be brighter than the majority of known high-redshift galaxies. OBGs could potentially be revealed via HST follow-up imaging of samples of brighter Lyman-break galaxies provided by wide-area ground-based surveys such as UltraVISTA, and should be easily uncovered and studied with instruments aboard JWST...(abridged)

We study the (M,sigma) and (M,Re) projections of the thin Mass Plane (MP) (M,sigma,Re) which describes the distribution of the galaxy population. The distribution of galaxy properties on the MP is characterized by: (i) a zone of exclusion described by two power-laws joined by a break at M 3e10 Msun. This results in a break in the mean M-sigma relation with Msigma^2.3 and M sigma^4.7 at small and large sigma respectively; (ii) a mass M 2e11 Msun which separates a population dominated by flat fast rotator with disks and spiral galaxies at lower masses, from one dominated by quite round slow rotators at larger masses; (iii) below that mass the distribution of ETGs properties tends to be constant along lines of constant sigma; (iv) it forms a parallel sequence with the distribution of spiral galaxies; (v) at even lower masses, the distribution of fast rotator ETGs and late spirals naturally extends to that of dwarf ETGs (Sph) and dwarf irregulars (Im) respectively. We show via dynamical models that sigma traces the bulge fraction, which drives the observed trends in M/L, Hbeta, colour, IMF and molecular gas fraction. We interpret this as due to a combination of two main effects: (i) an increase of the bulge fraction which increases sigma and greatly enhances the likelihood for a galaxy to have its star formation quenched, and (ii) dry merging, increasing galaxy mass along lines of nearly constant sigma, while leaving the population unchanged. [Abriged]

We study the volume-limited and nearly mass selected (stellar mass M>6*10^9 Msun) Atlas3D sample of 260 early-type galaxies. We construct detailed axisymmetric dynamical models (JAM), which allow for orbital anisotropy, include a dark matter halo, and reproduce in detail both the galaxy images and the high-quality integral-field stellar kinematics. We derive accurate total M/L and dark matter fractions f_DM, within a sphere of radius r=Re. We also measure the stellar M/L and derive a median dark matter fraction f_DM=13%. We find that the thin two-dimensional subset spanned by galaxies in the (M_JAM,sigma_e,R_e) coordinates system, which we call the Mass Plane (MP) has an observed rms scatter of 19% and an intrinsic one of 11%. The MP satisfies the scalar virial relation M_JAM sigma_e^2 R_e within our tight errors. However, the details of how both Re and sigma_e are determined are critical in defining the precise deviation from the virial exponents. We revisit the (M/L)-sigma_e relation, which describes most of the deviations between the MP and the FP. The best-fitting relation is (M/L) sigma_e^0.72 (r-band). It provides an upper limit to any systematic increase of the IMF mass normalization with sigma_e. We study of the link between sigma_e and the galaxies circular velocity V_circ within 1Re (where stars dominate) and find the relation max(V_circ)~1.76*sigma_e, which has an observed scatter of 7%. The accurate parameters described in this paper are used in the companion Paper XX of this series to explore the variation of global galaxy properties, including the IMF, on the projections of the MP. [Abridged]

For early-type galaxies, the ability to sustain a corona of hot, X-ray emitting gas could have played a key role in quenching their star-formation history. Yet, it is still unclear what drives the precise amount of hot gas around these galaxies. By combining photometric and spectroscopic measurements for the early-type galaxies observed during the Atlas3D integral-field survey with measurements of their X-ray luminosity based on X-ray data of both low and high spatial resolution we conclude that the hot-gas content of early-type galaxies can depend on their dynamical structure. Specifically, whereas slow rotators generally have X-ray halos with luminosity L_X,gas and temperature T values that are in line with what is expected if the hot-gas emission is sustained by the thermalisaton of the kinetic energy carried by the stellar-mass loss material, fast rotators tend to display L_X,gas values that fall consistently below the prediction of this model, with similar T values that do not scale with the stellar kinetic energy as observed in the case of slow rotators. Considering that fast rotators are likely to be intrinsically flatter than slow rotators, and that the few L_X,gas-deficient slow rotators also happen to be relatively flat, the observed L_X,gas deficiency in these objects would support the hypothesis whereby flatter galaxies have a harder time in retaining their hot gas. We discuss the implications that a different hot-gas content could have on the fate of both acquired and internally-produced gaseous material, considering in particular how the L_X,gas deficiency of fast rotators would make them more capable to recycle the stellar-mass loss material into new stars than slow rotators. This is consistent with the finding that molecular gas and young stars are detected only in fast rotators in the Atlas3D sample, and that fast rotators tend to dustier than slow rotators. [Abridged]

[Abridged] We present a detailed study of the physical properties of the molecular gas in a sample of 18 molecular gas-rich early-type galaxies (ETGs) from the ATLAS$ 3D sample. Our goal is to better understand the star formation processes occurring in those galaxies, starting here with the dense star-forming gas. We use existing integrated $^{12}$CO(1-0, 2-1), $^{13}$CO(1-0, 2-1), HCN(1-0) and HCO$^{+}$(1-0) observations and present new $^{12}$CO(3-2) single-dish data. From these, we derive for the first time the average kinetic temperature, H$_{2}$ volume density and column density of the emitting gas, this using a non-LTE theoretical model. Since the CO lines trace different physical conditions than of those the HCN and HCO$^{+}$ lines, the two sets of lines are treated separately. We also compare for the first time the predicted CO spectral line energy distributions (SLEDs) and gas properties of our molecular gas-rich ETGs with those of a sample of nearby well-studied disc galaxies. The gas excitation conditions in 13 of our 18 ETGs appear analogous to those in the centre of the Milky Way. Such results have never been obtained before for ETGs and open a new window to explore further star-formation processes in the Universe. The conclusions drawn should nevertheless be considered carefully, as they are based on a limited number of observations and on a simple model. In the near future, with higher CO transition observations, it should be possible to better identify the various gas components present in ETGs, as well as more precisely determine their associated physical conditions. To achieve these goals, we show here from our theoretical study, that mid-J CO lines (such as the $^{12}$CO(6-5) line) are particularly useful.

We study the global efficiency of star formation in high resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars 2 to 5 times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (above 10^4 cm-3) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and HI surface densities and their star formation rates as determined by their non-stellar 8um emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC524, a lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole Red Sequence, but can contribute to the reddening of some galaxies.

We explore the connection between the local escape velocity, V_esc, and the stellar population properties in the ATLAS3D survey, a complete, volume-limited sample of nearby early-type galaxies. We make use of ugriz photometry to construct Multi-Gaussian Expansion models of the surface brightnesses of our galaxies. We are able to fit the full range of surface brightness profiles found in our sample, and in addition we reproduce the results of state-of-the-art photometry in the literature with residuals of 0.04 mags. We utilise these photometric models and SAURON integral-field spectroscopy, combined with Jeans dynamical modelling, to determine the local V_esc derived from the surface brightness. We find that the local V_esc is tightly correlated with the Mgb and Fe5015 linestrengths and optical colours, and anti-correlated with the Hbeta linestrength. In the case of the Mgb and Colour - V_esc relations we find that the relation within individual galaxies follows the global relation between different galaxies. We intentionally ignored any uncertain contribution due to dark matter since we are seeking an empirical description of stellar population gradients in early-type galaxies that is ideal for quantitative comparison with model predictions. We identify a population of galaxies that occur only at low V_esc that exhibit negative gradients in the Mgb - and Colour -V_esc relations. These galaxies typically have young central stellar populations and contain significant amounts of molecular gas and dust. Combining these results with N-body simulations of binary mergers we use the Mgb-V_esc relation to constrain the possible number of dry mergers experienced by the local early-type galaxy population - a typical massive ETG can have experienced only ~1.5 major mergers before becoming a significant outlier in the Mgb-V_esc relation. [Abridged]

We detail the rich molecular story of NGC 1266, its serendipitous discovery within the ATLAS3D survey (Cappellari et al. 2011) and how it plays host to an AGN-driven molecular outflow, potentially quenching all of its star formation (SF) within the next 100 Myr. While major mergers appear to play a role in instigating outflows in other systems, deep imaging of NGC 1266 as well as stellar kinematic observations from SAURON, have failed to provide evidence that NGC 1266 has recently been involved in a major interaction. The molecular gas and the instantaneous SF tracers indicate that the current sites of star formation are located in a hypercompact disk within 200 pc of the nucleus (Fig. 1; SF rate ~ 2 Msuns/yr). On the other hand, tracers of recent star formation, such as the H{\beta} absorption map from SAURON and stellar population analysis show that the young stars are distributed throughout a larger area of the galaxy than current star formation. As the AGN at the center of NGC 1266 continues to drive cold gas out of the galaxy, we expect star formation rates to decline as the star formation is ultimately quenched. Thus, NGC 1266 is in the midst of a key portion of its evolution and continued studies of this unique galaxy may help improve our understanding of how galaxies transition from the blue to the red sequence (Alatalo et al. 2011).

NGC 1266 is a nearby field galaxy observed as part of the ATLAS3D survey (Cappellari et al. 2011). NGC 1266 has been shown to host a compact (< 200 pc) molecular disk and a mass-loaded molecular outflow driven by the AGN (Alatalo et al. 2011). Very Long Basline Array (VLBA) observations at 1.65 GHz revealed a compact (diameter < 1.2 pc), high bright- ness temperature continuum source most consistent with a low-level AGN origin. The VLBA continuum source is positioned at the center of the molecular disk and may be responsible for the expulsion of molecular gas in NGC 1266. Thus, the candidate AGN-driven molecular outflow in NGC 1266 supports the picture in which AGNs do play a significant role in the quenching of star formation and ultimately the evolution of the red sequence of galaxies.

[Abridged] We analyse the morphological structures in galaxies of the ATLAS3D sample by fitting a single Sersic profile and decomposing all non-barred objects (180 of 260 objects) in two components parameterised by an exponential and a general Sersic function. The aim of this analysis is to look for signatures of discs in light distributions of nearby early-type galaxies and compare them to kinematic properties. Using Sersic index from single component fits for a distinction between slow and fast rotators, or even late- and early-type galaxies, is not recommended. Assuming that objects with n>3 are slow rotators (or ellipticals), there is only a 22 per cent probability to correctly classify objects as slow rotators (or 37 per cent of previously classified as ellipticals). We show that exponential sub-components, as well as light profiles fitted with only a single component of a low Sersic index, can be linked with the kinematic evidence for discs in early-type galaxies. The median disk-to-total light ratio for fast and slow rotators is 0.41 and 0.0, respectively. Similarly, the median Sersic indices of the bulge (general Sersic component) are 1.7 and 4.8 for fast and slow rotators, respectively. Overall, discs or disc-like structures, are present in 83 per cent of early-type galaxies which do not have bars, and they show a full range of disk-to-total light ratios. Discs in early-type galaxies contribute with about 40 per cent to the total mass of the analysed (non-barred) objects. The decomposition into discs and bulges can be used as a rough approximation for the separation between fast and slow rotators, but it is not a substitute, as there is only a 59 per cent probability to correctly recognise slow rotators. Kinematics (i.e. projected angular momentum) remains the best approach to mitigate the influence of the inclination effects.

The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.

We use interferometric CO observations to compare the extent, surface brightness profiles and kinematics of the molecular gas in CO-rich Atlas3D early-type galaxies (ETGs) and spiral galaxies. We find that the molecular gas extent is smaller in absolute terms in ETGs than in late-type galaxies, but that the size distributions are similar once scaled by the galaxies optical/stellar characteristic scale-lengths. Virgo cluster ETGs have less extended molecular gas reservoirs than field counterparts. Approximately half of ETGs have molecular gas surface brightness profiles that follow the stellar light profile. These systems often have relaxed gas out to large radii, suggesting they are unlikely to have had recent merger/accretion events. A third of the sample galaxies show molecular gas surface brightness profiles that fall off slower than the light, and sometimes show a truncation. We suggest that ram pressure stripping and/or the presence of hot gas has compressed/truncated the gas in these systems. The remaining galaxies have rings, or composite profiles, that we argue can be caused by the effects of bars. We investigated the kinematics of the molecular gas using position-velocity diagrams, and compared the observed kinematics with dynamical model predictions, and the observed stellar and ionised gas velocities. We confirm that the molecular gas reaches beyond the turnover of the circular velocity curve in ~70% of our CO-rich Atlas3D ETGs. In general we find that in most galaxies the molecular gas is relaxed and dynamically cold. Molecular gas is a better direct tracer of the circular velocity than the ionised gas, justifying its use as a kinematic tracer for Tully-Fisher and similar analyses. (abridged)

We present the Combined Array for Research in Millimeter Astronomy (CARMA) ATLAS3D molecular gas imaging survey, a systematic study of the distribution and kinematics of molecular gas in CO-rich early-type galaxies. Our full sample of 40 galaxies (30 newly mapped and 10 taken from the literature) is complete to a 12CO(1-0) integrated flux of 18.5 Jy km/s, and it represents the largest, best-studied sample of its type to date. A comparison of the CO distribution of each galaxy to the g-r color image (representing dust) shows that the molecular gas and dust distributions are in good agreement and trace the same underlying interstellar medium. The galaxies exhibit a variety of CO morphologies, including discs (50%), rings (15%), bars+rings (10%), spiral arms (5%), and mildly (12.5%) and strongly (7.5%) disrupted morphologies. There appear to be weak trends between galaxy mass and CO morphology, whereby the most massive galaxies in the sample tend to have molecular gas in a disc morphology. We derive a lower limit to the total accreted molecular gas mass across the sample of 2.48x10^10 Msuns, or approximately 8.3x10^8 Msuns per minor merger within the sample, consistent with minor merger stellar mass ratios.

We report the discovery of a giant HI tail in the intra-group medium of HCG 44 as part of the Atlas3D survey. The tail is ~300 kpc long in projection and contains ~5x10^8 M_sun of HI. We detect no diffuse stellar light at the location of the tail down to ~28.5 mag/arcsec^2 in g band. We speculate that the tail might have formed as gas was stripped from the outer regions of NGC 3187 (a member of HCG 44) by the group tidal field. In this case, a simple model indicates that about 1/3 of the galaxy's HI was stripped during a time interval of <1 Gyr. Alternatively, the tail may be the remnant of an interaction between HCG 44 and NGC 3162, a spiral galaxy now ~650 kpc away from the group. Regardless of the precise formation mechanism, the detected HI tail shows for the first time direct evidence of gas stripping in HCG 44. It also highlights that deep HI observations over a large field are needed to gather a complete census of this kind of events in the local Universe.

The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue that these stars, {as well as many more abundant bound OB stars in the innermost few parsecs,} were generated by the interactions of an AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium.

We use an empirical approach to model the stellar mass of galaxies according to their host dark-matter haloes and subhaloes ('HASH'), where each galaxy resides in a subhalo taken from a large N-body cosmological simulation. This approach allows us to study the mass relation between subhaloes and galaxies (MR) using: the weak lensing signal (WL), the two point auto-correlation function (CF), and the stellar mass function of galaxies (SMF). Our method is based on modeling the lensing signal directly from the cosmological simulation, and should thus be more accurate than other methods based on analytic halo models. We find that the WL does not provide a strong constraint on the MR, allowing for more than a factor of 10 freedom in the subhalo mass of central galaxies, for a given stellar mass. The freedom for satellite galaxies is much larger, providing a very poor constraint on the fraction of satellite galaxies (0.05-0.8). These results are not significantly modified when using both the SMF and WL as constraints. We show that for the most massive galaxies, observational constraints based on the CF with 0.1 dex errors, are equivalent to 0.05 dex error in the WL. For intermediate and low mass galaxies the WL and CF constrain the MR in a different way. Although the WL is currently not adding much information at these masses, it has the potential of being important using future, more accurate measurements. The models found here do not match simultaneously the observed CF and WL signals, and show a limited ability to match the WL & SMF. We suspect that this is partially due to the cosmological model assumed here and we therefore adopt a mock WL signal through most of this work. In comparison to previous models in the literature, the method presented here is probably more general, as it reveals a larger range of solutions for a given set of observational constraints [abridged].

We study for the first time the environment of massive black hole (BH) seeds (~10^4-5 Msun) formed via the direct collapse of pristine gas clouds in massive haloes (>10^7 Msun) at z>6. Our model is based on the evolution of dark matter haloes within a cosmological N-body simulation, combined with prescriptions for the formation of BH along with both Pop III and Pop II stars. We calculate the spatially-varying intensity of Lyman Werner (LW) radiation from stars and identify the massive pristine haloes in which it is high enough to shut down molecular hydrogen cooling. In contrast to previous BH seeding models with a spatially constant LW background, we find that the intensity of LW radiation due to local sources, J_local, can be up to 10^6 times the spatially averaged background in the simulated volume and exceeds the critical value, J_crit, for the complete suppression of molecular cooling, in some cases by 4 orders of magnitude. Even after accounting for possible metal pollution in a halo from previous episodes of star formation, we find a steady rise in the formation rate of direct collapse (DC) BHs with decreasing redshift from 10^{-3}/Mpc^3/z at z=12 to 10^{-2}/Mpc^3/z at z=6. The onset of Pop II star formation at z~16 simultaneously marks the onset of the epoch of DCBH formation, as the increased level of LW radiation from Pop II stars is able to elevate the local levels of the LW intensity to J_local > J_crit while Pop III stars fail to do so at any time. The number density of DCBHs is sensitive to the number of LW photons and can vary by an order of magnitude at z=6 after accounting for reionisation feedback. Haloes hosting DCBHs are more clustered than similar massive counterparts that do not host DCBHs, especially at redshifts z>10. We also show that planned surveys with JWST should be able to detect the supermassive stellar precursors of DCBHs.

We use the SAURON and GMOS integral field spectrographs to observe the active galactic nucleus (AGN) powered outflow in NGC 1266. This unusual galaxy is relatively nearby (D=30 Mpc), allowing us to investigate the process of AGN feedback in action. We present maps of the kinematics and line strengths of the ionised gas emission lines Halpha, Hbeta, [OIII], [OI], [NII] and [SII], and report on the detection of Sodium D absorption. We use these tracers to explore the structure of the source, derive the ionised and atomic gas kinematics and investigate the gas excitation and physical conditions. NGC 1266 contains two ionised gas components along most lines of sight, tracing the ongoing outflow and a component closer to the galaxy systemic, the origin of which is unclear. This gas appears to be disturbed by a nascent AGN jet. We confirm that the outflow in NGC 1266 is truly multiphase, containing radio plasma, atomic, molecular and ionised gas and X-ray emitting plasma. The outflow has velocities up to \pm900 km/s away from the systemic velocity, and is very likely to be removing significant amounts of cold gas from the galaxy. The LINER-like line-emission in NGC 1266 is extended, and likely arises from fast shocks caused by the interaction of the radio jet with the ISM. These shocks have velocities of up to 800 km/s, which match well with the observed velocity of the outflow. Sodium D equivalent width profiles are used to set constraints on the size and orientation of the outflow. The ionised gas morphology correlates with the nascent radio jets observed in 1.4 GHz and 5 GHz continuum emission, supporting the suggestion that an AGN jet is providing the energy required to drive the outflow.

With the first metal enrichment by Population (Pop) III supernovae (SNe), the formation of the first metal-enriched, Pop II stars becomes possible. In turn, Pop III star formation and early metal enrichment are slowed by the high energy radiation emitted by Pop II stars. Thus, through the SNe and radiation they produce, Populations II and III coevolve in the early Universe, one regulated by the other. We present large (4 Mpc)^3, high resolution cosmological simulations in which we self-consistently model early metal enrichment and the stellar radiation responsible for the destruction of the coolants (H2 and HD) required for Pop III star formation. We find that the molecule-dissociating stellar radiation produced both locally and over cosmological distances reduces the Pop III star formation rate at z > 10 by up to an order of magnitude compared to the case in which this radiation is not included. However, we find that the effect of LW feedback is to enhance the amount of Pop II star formation. We attribute this to the reduced rate at which gas is blown out of dark matter haloes by SNe in the simulation with LW feedback, which results in larger reservoirs for metal-enriched star formation. Even accounting for metal enrichment, molecule-dissociating radiation and the strong suppression of low-mass galaxy formation due to reionization at z < 10, we find that Pop III stars are still formed at a rate of ~ 10^-5 M_sun yr^-1 Mpc^-3 down to z ~ 6. This suggests that the majority of primordial pair-instability SNe that may be uncovered in future surveys will be found at z < 10. We also find that the molecule-dissociating radiation emitted from Pop II stars may destroy H2 molecules at a high enough rate to suppress gas cooling and allow for the formation of supermassive primordial stars which collapse to form ~ 100,000 solar mass black holes.

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts.

(Abridged) We present here a detailed analysis of the star formation history (SFH) of FW4871, a massive galaxy at z=1.893+-0.002. We compare rest-frame optical and NUV slitless grism spectra from the Hubble Space Telescope with a large set of composite stellar populations to constrain the underlying star formation history. Even though the morphology features prominent tidal tails, indicative of a recent merger, there is no sign of on-going star formation within an aperture encircling one effective radius, which corresponds to a physical extent of 2.6 kpc. A model assuming truncation of an otherwise constant SFH gives a formation epoch zF~10, with a truncation after 2.7 Gyr, giving a mass-weighted age of 1.5 Gyr and a stellar mass of 0.8-3E11Msun, implying star formation rates of 30-110 Msun/yr. A more complex model including a recent burst of star formation places the age of the youngest component at 145 Myr, with a mass contribution lower than 20%, and a maximum amount of dust reddening of E(B-V)<0.4 mag (95% confidence levels). This low level of dust reddening is consistent with the low emission observed at 24 micron, corresponding to rest-frame 8 micron, where PAH emission should contribute significantly if a strong formation episode were present. The colour profile of FW4871 does not suggest a significant radial trend in the properties of the stellar populations out to 3Re. We suggest that the recent merger that formed FW4871 is responsible for the quenching of its star formation.

Feedback from active galactic nuclei (AGN) has become a major component in simulations of galaxy evolution, in particular for massive galaxies. AGN jets have been shown to provide a large amount of energy and are capable of quenching cooling flows. Their impact on the host galaxy, however, is still not understood. Subgrid models of AGN activity in a galaxy evolution context so far have been mostly focused on the quenching of star formation. To shed more light on the actual physics of the "radio mode" part of AGN activity, we have performed simulations of the interaction of a powerful AGN jet with the massive gaseous disc (10^11 solar masses) of a high-redshift galaxy. We spatially resolve both the jet and the clumpy, multi-phase interstellar medium (ISM) and include an explicit star formation model in the simulation. Following the system over more than 10^7 years, we find that the jet activity excavates the central region, but overall causes a significant change to the shape of the density probability distribution function and hence the star formation rate due to the formation of a blast wave with strong compression and cooling in the ISM. This results in a ring- or disc-shaped population of young stars. At later times, the increase in star formation rate also occurs in the disc regions further out since the jet cocoon pressurizes the ISM. The total mass of the additionally formed stars may be up to 10^10 solar masses for one duty cycle. We discuss the details of this jet-induced star formation (positive feedback) and its potential consequences for galaxy evolution and observable signatures.

We investigate the accuracy in the recovery of the stellar dynamics of barred galaxies when using axisymmetric dynamical models. We do this by trying to recover the mass-to-light ratio (M/L) and the anisotropy of realistic galaxy simulations using the Jeans Anisotropic Multi-Gaussian Expansion (JAM) method. However, given that the biases we find are mostly due to an application of an axisymmetric modeling algorithm to a non-axisymmetric system and in particular to inaccuracies in the de-projected mass model, our results are relevant for general axisymmetric modelling methods. We run N-body collisionless simulations to build a library with various luminosity distribution, constructed to mimic real individual galaxies, with realistic anisotropy. The final result of our evolved library of simulations contains both barred and unbarred galaxies. The JAM method assumes an axisymmetric mass distribution, and we adopt a spatially constant M/L and anisotropy beta_z=1-sigma_z^2/sigma_R^2 distributions. The models are fitted to two-dimensional maps of the second velocity moments V_rms=sqrt(V^2+sigma^2) of the simulations for various viewing angles (position angle of the bar and inclination of the galaxy). We find that the inclination is generally well recovered by the JAM models, for both barred and unbarred simulations. For unbarred simulations the M/L is also accurately recovered, with negligible median bias and with a maximum one of just Delta(M/L)<1.5% when the galaxy is not too close to face on. At very low inclinations (i<30 deg) the M/L can be significantly overestimated (9% in our tests, but errors can be larger for very face-on views). For barred simulations the M/L is on average (when PA=45 deg) essentially unbiased, but we measure an over/under estimation of up to Delta(M/L)=15% in our tests. The sign of the M/L bias depends on the position angle of the bar as expected. [Abridged]

Much of our knowledge of galaxies comes from analysing the radiation emitted by their stars. It depends on the stellar initial mass function (IMF) describing the distribution of stellar masses when the population formed. Consequently knowledge of the IMF is critical to virtually every aspect of galaxy evolution. More than half a century after the first IMF determination, no consensus has emerged on whether it is universal in different galaxies. Previous studies indicated that the IMF and the dark matter fraction in galaxy centres cannot be both universal, but they could not break the degeneracy between the two effects. Only recently indications were found that massive elliptical galaxies may not have the same IMF as our Milky Way. Here we report unambiguous evidence for a strong systematic variation of the IMF in early-type galaxies as a function of their stellar mass-to-light ratio, producing differences up to a factor of three in mass. This was inferred from detailed dynamical models of the two-dimensional stellar kinematics for the large Atlas3D representative sample of nearby early-type galaxies spanning two orders of magnitude in stellar mass. Our finding indicates that the IMF depends intimately on a galaxy's formation history.

To answer questions on the start and duration of the epoch of reionisation, periods of galaxy mergers and properties of other cosmological encounters, the cosmic star formation history (CSFH), is of fundamental importance. Using the association of long gamma-ray bursts (LGRBs) with the death of massive stars and their ultra-luminous nature, the CSFH can be probed to higher redshifts than current conventional methods. Unfortunately, no consensus has been reached on the manner in which the LGRB rate (LGRBR) traces the CSFH, leaving many of the questions mentioned mostly unexplored by this method. Observations by the GRB NIR detector (GROND) over the past 4 years have, for the first time, acquired highly complete LGRB samples. Driven by these completeness levels and new evidence of LGRBs also occurring in more massive and metal rich galaxies than previously thought, the possible biases of the LGRBR-CSFH connection are investigated over a large range of galaxy properties. The CSFH is modelled using empirical fits to the galaxy mass function and galaxy star formation rates. Biasing the CSFH by metallicity cuts, mass range boundaries, and other unknown redshift dependencies, a LGRBR is generated and compared to the highly complete GROND sample. It is found that there is no strong preference for a metallicity cut or fixed galaxy mass boundaries and that there are no unknown redshift effects, in contrast to previous work which suggest values of Z/Z_sun~0.1-0.3. From the best-fit models, we predict that ~1.2% of the LGRB burst sample exists above z=6. The linear relationship between the LGRBR and the CSFH suggested by our results implies that redshift biases present in previous LGRB samples significantly affect the inferred dependencies of LGRBs on their host galaxy properties. Such biases can lead to, e.g., an interpretation of metallicity limitations and evolving LGRB luminosity functions.

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmological simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feedback given by the HYD, as functions of the halo mass and redshift. Surprisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with deviations at the level of 0.1 dex, on an object-by-object basis. This is true for all redshifts, and for the mass of stars and gas components, although the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same level of accuracy is obtained even in case the SAM uses only one phase of gas within each galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides sufficient information for the SAM to reproduce the population of galaxies within the entire cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the evolution within merger-trees agrees on average. b) For a given galaxy, efficiencies fluctuate around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galaxies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We compare the efficiencies found here to standard SAM recipes and find that they often deviate significantly. For example, here the HYD shows smooth accretion that is less effective for low mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high redshift; and star formation changes only mildly with cosmic time. The method developed here can be applied in general to any HYD, and can thus serve as a common language for both HYDs and SAM (Abridged).

We study the formation of early-type galaxies (ETGs) through mergers with a sample of 70 high-resolution numerical simulations of binary mergers of disc galaxies. These simulations encompass various mass ratios, initial conditions and orbital parameters. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion: they preserve the structure of the input fast rotating spiral progenitors. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Slow Rotators hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in binary mergers, but it also requires a retrograde spin for the progenitor galaxies with respect to the orbital angular momentum. The importance of the initial spin of the progenitors is also investigated in the library of galaxy mergers of the GalMer project.

Surveying eighteen 12CO-bright galaxies from the ATLAS3D early-type galaxy sample with the Institut de Radio Astronomie Millim\'etrique (IRAM) 30m telescope, we detect 13CO(1-0) and 13CO(2-1) in all eighteen galaxies, HCN(1-0) in 12/18 and HCO+(1-0) in 10/18. We find that the line ratios 12CO(1-0)/13CO(1-0) and 12CO(1-0)/HCN(1-0) are clearly correlated with several galaxy properties: total stellar mass, luminosity-weighted mean stellar age, molecular to atomic gas ratio, dust temperature and dust morphology. We suggest that these correlations are primarily governed by the optical depth in the 12CO lines; interacting, accreting and/or starbursting early-type galaxies have more optically thin molecular gas while those with settled dust and gas discs host optically thick molecular gas. The ranges of the integrated line intensity ratios generally overlap with those of spirals, although we note some outliers in the 12CO(1- 0)/13CO(1-0), 12CO(2-1)/13CO(2-1) and HCN/HCO+(1-0) ratios. In particular, three galaxies are found to have very low 12CO(1-0)/13CO(1-0) and 12CO(2-1)/13CO(2-1) ratios. Such low ratios may signal particularly stable molecular gas which creates stars less efficiently than 'normal' (i.e. below Schmidt-Kennicutt prediction), consistent with the low dust temperatures seen in these galaxies.

We study via numerical N-body/SPH chemistry simulations the effects of primordial non-Gaussianities on the formation of the first stars and galaxies, and investigate the impact of supernova feedback in cosmologies with different fnl. Density distributions are biased to higher values, so star formation and the consequent feedback processes take place earlier in high-fnl models and later in low-fnl ones. Mechanical feedback is responsible for shocking and evacuating the gas from star forming sites earlier in the highly non-Gaussian cases, because of the larger bias at high densities. Chemical feedback translates into high-redshift metal filling factors that are larger by some orders of magnitude for larger fnl, but that converge within one Gyr, for both population III and population II-I stellar regimes. The efficient enrichment process, though, leads to metallicities > 0.01 Zsun by redshift ~9, almost independently from fnl. The impact of non-Gaussianities on the formation of dark-matter haloes at high redshift is directly reflected in the properties of the gas in these haloes, as models with larger fnl show more concentrated gas profiles at early times. Non-Gaussian signatures in the gas behaviour are lost after the first feedback takes place and introduces a significant degree of turbulence and chaotic motions.

We present the Atlas3D HI survey of 166 nearby early-type galaxies (ETGs) down to M(HI)~10^7 M_sun. We detect HI in ~40% of all ETGs outside the Virgo cluster and in ~10% of all ETGs inside it. This demonstrates that it is common for non-cluster ETGs to host HI. The HI morphology varies from regular discs/rings (the majority of the detections) to unsettled gas distributions. The former are either small discs (M(HI)<10^8 M_sun) confined within the stellar body and sharing the same kinematics of the stars, or large discs/rings (M(HI) up to 5x10^9 M_sun) extending to tens of kpc from the host galaxy and frequently kinematically decoupled from the stars. Neutral hydrogen provides material for star formation in ETGs. Galaxies with central HI exhibit signatures of star formation in ~70% of the cases, ~5 times more frequently than galaxies without central HI. The central ISM is dominated by molecular gas. In ETGs with a small gas disc the conversion of HI into H_2 is as efficient as in spirals. The ETG HI mass function has M*~2x10^9 M_sun and slope=-0.7. ETGs host much less HI than spirals as a family. However, a significant fraction of them is as HI-rich as spirals. The main difference between ETGs and spirals is that the former lack the high-column-density HI typical of the bright stellar disc of the latter. We find an envelope of decreasing M(HI) with increasing environment density. The gas-richest ETGs live in the poorest environments (where star-formation is more common), galaxies in the centre of Virgo have the lowest HI content, and the cluster outskirts are a transition region. We find an HI morphology-density relation. At low environment density HI is mostly distributed on large discs/rings. More disturbed HI morphologies dominate environment densities typical of rich groups, confirming the importance of processes occurring on a galaxy-group scale for the evolution of ETGs.

While stars are widely discussed as the source of the high energy photons which reionized the universe, an additional source of ionizing photons that must also contribute to reionization in this scenario is the supernovae (SNe) which mark the end of the life of massive stars. Here we estimate the relative contributions of SNe and stars to reionization. While the rate at which ionizing photons are produced in SN shocks is well below that at which they are produced by stars, the harder spectra of radiation emitted from SNe leads to an enhanced escape fraction of SN generated photons relative to that of stellar photons. In particular, along a given line of sight out of a galaxy, we find that for neutral hydrogen column densities N_H >~10^18 cm^-2 the contribution to reionization from SNe is greater than that from stars. Drawing on the results of simulations presented in the literature, we find that the overall (line of sight-averaged) SNe shock-generated ionizing photon escape fraction is larger than the stellar photon escape fraction by a factor of \simeq 4 to \simeq 7, depending on the metallicity of the stellar population. Overall, our results suggest that the effect of SNe is an enhancement of up to ~ 10 percent in the fraction of hydrogen reionized by stellar sources. We briefly discuss the implications of our results for the population of galaxies responsible for reionization.

We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical $\Lambda$CDM scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via $\lambda_R$. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the \atl survey, assuming that fast rotating ETGs have at least 10% of their total stellar mass in a disc component. In agreement with \atl observations we find that slow rotators are predominantly galaxies with $ M_* > 10^{10.5}$ M$_{\odot}$ contributing $\sim 20%$ to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between $\sim 50% -90%$ of their stellar mass from satellites and their most massive progenitors have on average up to 3 major mergers during their evolution. Fast rotators in contrast, accrete less than 50% and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stops, galaxies grow via infall from satellites. Frequent minor mergers thereby, destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering $\lambda_R$ into the slow rotator regime. Abridged...

We make use of interferometric CO and HI observations, and optical integral-field spectroscopy to probe the origin of the molecular and ionised interstellar medium (ISM) in local early-type galaxies (ETGs). We find that 36\pm5% of our sample of fast rotating ETGs have their ionised gas kinematically misaligned with respect to the stars, setting a strong lower limit on the importance of externally acquired gas (e.g. from mergers and cold accretion). Slow rotators have a flat distribution of misalignments, indicating that the dominant source of gas is external. The molecular, ionised and atomic gas in all the detected galaxies are always kinematically aligned, even when they are misaligned from the stars, suggesting that all these three phases of the ISM share a common origin. In addition, we find that the origin of the cold and warm gas in fast-rotating ETGs is strongly affected by environment, despite the molecular gas detection rate and mass fractions being fairly independent of group/cluster membership. Galaxies in dense groups and the Virgo cluster nearly always have their molecular gas kinematically aligned with the stellar kinematics, consistent with a purely internal origin. In the field, however, kinematic misalignments between the stellar and gaseous components indicate that >46% of local fast-rotating ETGs have their gas supplied from external sources. We discuss several scenarios which could explain the environmental dichotomy, but find it difficult to simultaneously explain the kinematic misalignment difference and the constant detection rate. Furthermore, our results suggest that galaxy mass may be an important independent factor associated with the origin of the gas, with the most massive fast-rotating galaxies in our sample (M_K<-24 mag; stellar mass of >8x10^10 Msun) always having kinematically aligned gas. (abridged)

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

The mass assembly of galaxies leaves imprints in their outskirts, such as shells and tidal tails. The frequency and properties of such fine structures depend on the main acting mechanisms - secular evolution, minor or major mergers - and on the age of the last substantial accretion event. We use this to constrain the mass assembly history of two apparently relaxed nearby Early-Type Galaxies (ETGs) selected from the Atlas-3D sample, NGC 680 and NGC 5557. Our ultra deep optical images obtained with MegaCam on the Canada-France-Hawaii Telescope reach 29 mag/arcsec^2 in the g-band. They reveal very low-surface brightness (LSB) filamentary structures around these ellipticals. Among them, a gigantic 160 kpc long tail East of NGC 5557 hosts gas-rich star-forming objects. NGC 680 exhibits two major diffuse plumes apparently connected to extended HI tails, as well as a series of arcs and shells. Comparing the outer stellar and gaseous morphology of the two ellipticals with that predicted from models of colliding galaxies, we argue that the LSB features are tidal debris, the star-forming objects near NGC 5557, long lived Tidal Dwarf Galaxies and that each of the two ETGs was assembled during a relatively recent, major wet merger, which likely occurred at a redshift below z = 0.5. The inner kinematics of NGC 680 is typical for fast rotators which make the bulk of nearby ETGs in the Atlas-3D sample. NGC 5557 belongs to the poorly populated class of massive, round, slow rotators that are predicted by semi-analytic models and cosmological simulations to be the end-product of a complex mass accretion history, involving ancient major mergers and more recent minor mergers. Our observations suggest that under specific circumstances a single binary merger may dominate the formation history of such objects and thus that at least some massive ETGs may form at relatively low redshift (abridged).

We study the formation of early-type galaxies through mergers with a sample of 70 high-resolution (softening length < 60 pc and 12*10^6 particles) numerical simulations of binary mergers of disc galaxies and 16 simulations of ETG remergers. These simulations, designed to accompany observations and models conducted within the Atlas3D project, encompass various mass ratios (from 1:1 to 6:1), initial conditions and orbital parameters. The progenitor disc galaxies are spiral-like with bulge to disc ratios typical of Sb and Sc galaxies. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion (based on the lambda_R parameter): they preserve the structure of the input fast rotating spiral progenitors. They have intrinsic ellipticities larger than 0.5, cover intrinsic lambda_R values between 0.2 and 0.6, within the range of observed Fast Rotators. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Fast Rotators produced in major mergers have intermediate flattening, with ellipticities between 0.4 and 0.6. Most Slow Rotators formed in these binary disc mergers hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. Besides a handful of specific observed systems -- the counter-rotating discs (2-sigma galaxies) -- these therefore cannot reproduce the observed population of Slow Rotators in the nearby Universe. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in these binary mergers, but it also requires a retrograde spin for the earlier-type (Sb) progenitor galaxy with respect to the orbital angular momentum. (Abridged)

We look at the morphology of fast and slow rotator early-type galaxies. Edge-on fast rotators are lenticular galaxies. They appear like spiral galaxies with the gas and dust removed, and in some cases are flat ellipticals with disky isophotes. Fast rotators are often barred and span the same full range of bulge fractions as spiral galaxies. The slow rotators are rounder and are generally consistent with being genuine, namely spheroidal-like, elliptical galaxies. We propose a revision to the tuning-fork diagram by Hubble as it gives a misleading description of ETGs. We study for the first time the kinematic morphology-density T-Sigma relation using fast and slow rotators to replace lenticulars and ellipticals. We find that our relation is cleaner than using classic morphology. Slow rotators are nearly absent at the lowest density environments [f(SR)<2%] and generally constitute a small fraction [f(SR)~4%] of the total galaxy population in the environments explored by our survey, with the exception of the densest core of the Virgo cluster [f(SR)~20%]. We find a clean log-linear relation between the fraction f(Sp) of spiral galaxies and the local galaxy surface density. The existence of a smooth kinematic T-Sigma relation in the field excludes processes related to the cluster environment as main contributors to the apparent conversion of spirals into fast-rotators in low-density environments. It shows that the segregation is driven by local effects at the small-group scale. Only at the largest densities in the Virgo core does the f(Sp) relation break down and steepens sharply, while the fraction of slow-rotators starts to significantly increase. This suggests that a different mechanism is at work there. (Abridged)

We report the discovery of a powerful molecular wind from the nucleus of the non-interacting nearby S0 field galaxy NGC 1266. The single-dish CO profile exhibits emission to +/- 400 km/s and requires a nested Gaussian fit to be properly described. Interferometric observations reveal a massive, centrally- concentrated molecular component with a mass of 1.1x10^9 Msuns and a molecular outflow with a molecular mass of 2.4x10^7 Msuns. The molecular gas close to the systemic velocity consists of a rotating, compact nucleus with a mass of about 4.1x108 Msuns within a radius of approximately 60 pc. This compact molecular nucleus has a surface density of \approx 2.7 \times 10^4 Msuns/pc^2, more than two orders of magnitude larger than that of giant molecular clouds in the disk of the Milky Way, and it appears to sit on the Kennicutt-Schmidt relation despite its extreme kinematics and energetic activity. We interpret this nucleus as a disk that confines the outflowing wind. A mass outflow rate of 13 Msuns/yr leads to a depletion timescale of about 85 Myr. The star formation in NGC 1266 is insufficient to drive the outflow, and thus it is likely driven by the active galactic nucleus (AGN). The concentration of the majority of the molecular gas in the central 100 pc requires an extraordinary loss of angular momentum, but no obvious companion or interacting galaxy is present to enable the transfer. NGC 1266 is the first known outflowing molecular system that does not show any evidence of a recent interaction.

[Abridged] We use the ATLAS3D sample of 260 early-type galaxies (ETGs) to study the apparent kinematic misalignment angle, Psi, defined as the angle between the photometric and kinematic major axis. We find that 71% of nearby ETGs are strictly aligned systems (Psi > 5 deg), an additional 14% have 5 < Psi < 10 deg and 90% of galaxies have Psi < 15 deg. Taking into account measurement uncertainties, 90% of galaxies can be considered aligned to better than 5 deg, suggesting that only a small fraction of early-type galaxies (~10%) are not consistent with axisymmetry within the projected half-light radius. We identify morphological features such as bars and rings (30%), dust structures (16%), blue nuclear colours (6%) and evidence of interactions (8%). We use kinemetry to analyse the mean velocity maps and separate galaxies in two broad types of regular and non-regular rotators. We find 82% of regular rotators and 17% non-regular rotators, with 2 galaxies that we were not able to classify due to data quality. The non-regular rotators are typically found in dense regions and are massive. The majority of galaxies do not have any specific kinemetric features, but we highlight here the frequency of the kinematically distinct cores (7%) and the aligned double peaks in the velocity dispersion maps (4%). Most of the galaxies that are misaligned have complex kinematics and are non-regular rotators. While the trends are weak, there is a tendency that large values of Psi are found in galaxies at intermediate environmental densities and among the most massive galaxies in the sample. We suggest that the most common formation mechanism for ETGs preserves the axisymmetry of the disk progenitors. Less commonly, the formation process results in a triaxial galaxy with much lower net angular momentum.

We have carried out a survey for 12CO J=1-0 and J=2-1 emission in the 260 early-type galaxies of the volume-limited Atlas3D sample, with the goal of connecting their star formation and assembly histories to their cold gas content. This is the largest volume-limited CO survey of its kind and is the first to include many Virgo Cluster members. Sample members are dynamically hot galaxies with a median stellar mass 3\times 10^{10} Msun; they are selected by morphology rather than colour, and the bulk of them lie on the red sequence. The overall CO detection rate is 56/259 = 0.22 \error 0.03, with no dependence on K luminosity and only a modest dependence on dynamical mass. There are a dozen CO detections among the Virgo Cluster members; statistical analysis of their H_2 mass distributions and their dynamical status within the cluster shows that the cluster's influence on their molecular masses is subtle at best, even though (unlike spirals) they seem to be virialized within the cluster. We suggest that the cluster members have retained their molecular gas through several Gyr residences in the cluster. There are also a few extremely CO-rich early-type galaxies with H_2 masses >= 10^9 Msun, and these are in low density environments. We do find a significant trend between molecular content and the stellar specific angular momentum. The galaxies of low angular momentum also have low CO detection rates, suggesting that their formation processes were more effective at destroying molecular gas or preventing its re-accretion. We speculate on the implications of these data for the formation of various sub-classes of early-type galaxies.

We demonstrate here using both single-dish and interferometric observations that CO molecules are an excellent kinematic tracer, even in high-mass galaxies, allowing us to investigate for the first time the CO Tully-Fisher relation of early-type galaxies. We compare the Tully-Fisher relations produced using both single-dish and interferometric data and various inclination estimation methods, and evaluate the use of the velocity profile shape as a criterion for selecting galaxies in which the molecular gas extends beyond the peak of the rotation curve. We show that the gradient and zero-point of the best-fit relations are robust, independent of the velocity measure and inclination used, and agree with those of relations derived using stellar kinematics. We also show that the early-type CO Tully-Fisher (CO-TF) relation is offset from the CO-TF of spirals by 0.98 \pm 0.22 magnitude at Ks-band, in line with other results. The intrinsic scatter of the relation is found to be \approx 0.4 magnitudes, similar to the level found in the spiral galaxy population. Next generation facilities such as the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/Sub-millimeter Array (ALMA) should allow this technique to be used in higher-redshift systems, providing a simple new tool to trace the mass-to-light ratio evolution of the most massive galaxies over cosmic time.

We provide a census of the apparent stellar angular momentum within 1 Re of a volume-limited sample of 260 early-type galaxies (ETGs) in the nearby Universe, using integral-field spectroscopy obtained in the course of the ATLAS3D project. We exploit the LambdaR parameter to characterise the existence of two families of ETGs: Slow Rotators which exhibit complex stellar velocity fields and often include stellar kinematically Distinct Cores (KDCs), and Fast Rotators which have regular velocity fields. Our complete sample of 260 ETGs leads to a new criterion to disentangle Fast and Slow Rotators which now includes a dependency on the apparent ellipticity (Epsilon). It separates the two classes significantly better than the previous prescription, and than a criterion based on V/Sigma: Slow Rotators and Fast Rotators have LambdaR lower and larger than kFSxSQRT(Epsilon), respectively, where kFS=0.31 for measurements made within 1 Re. We show that the vast majority of early-type galaxies are Fast Rotators: these have regular stellar rotation, with aligned photometric and kinematic axes (Paper II, Krajnovic et al. 2011}, include discs and often bars and represent 86% (224/260) of all early-type galaxies in the volume-limited ATLAS3D sample. Fast Rotators span the full range of apparent ellipticities from 0 to 0.85, and we suggest that they cover intrinsic ellipticities from about 0.35 to 0.85, the most flattened having morphologies consistent with spiral galaxies. Only a small fraction of ETGs are Slow Rotators representing 14% (36/260) of the ATLAS3D sample of ETGs. Of all Slow Rotators, 11% (4/36) exhibit two counter-rotating stellar disc-like components and are rather low mass objects (Mdyn<10^10.5 M_Sun). All other Slow Rotators (32/36) appear relatively round on the sky (Epsilon_e<0.4), tend to be massive (Mdyn>10^10.5 M_Sun), and often (17/32) exhibit KDCs.

The molecular gas content of local early-type galaxies is constrained and discussed in relation to their evolution. First, as part of the Atlas3D survey, we present the first complete, large (260 objects), volume-limited single-dish survey of CO in normal local early-type galaxies. We find a surprisingly high detection rate of 22%, independent of luminosity and at best weakly dependent on environment. Second, the extent of the molecular gas is constrained with CO synthesis imaging, and a variety of morphologies is revealed. The kinematics of the molecular gas and stars are often misaligned, implying an external gas origin in over a third of the systems, although this behaviour is drastically diffferent between field and cluster environments. Third, many objects appear to be in the process of forming regular kpc-size decoupled disks, and a star formation sequence can be sketched by piecing together multi-wavelength information on the molecular gas, current star formation, and young stars. Last, early-type galaxies do not seem to systematically obey all our usual prejudices regarding star formation, following the standard Schmidt-Kennicutt law but not the far infrared-radio correlation. This may suggest a greater diversity in star formation processes than observed in disk galaxies. Using multiple molecular tracers, we are thus starting to probe the physical conditions of the cold gas in early-types.

We study cosmological simulations of early structure formation, including non-equilibrium molecular chemistry, metal pollution from stellar evolution, transition from population III (popIII) to population II (popII) star formation, regulated by a given critical metallicity, and feedback effects. We investigate the properties of early metal spreading from the different stellar populations and its interplay with primordial molecular gas. We find that, independently of the details about popIII modeling, after the onset of star formation, regions enriched below the critical level are mostly found in isolated environments, while popII star formation regions are much more clumped. Typical star forming haloes show average SN driven outflow rates of up to 10^{-4} Msun/yr in enriched gas, initially leaving the original star formation regions almost devoid of metals. The polluted material, which is gravitationally incorporated in over-dense environments on timescales of 10^7 yr, is mostly coming from external, nearby star forming sites ("gravitational enrichment"). In parallel, the pristine-gas inflow rates are between 10^{-3} - 10^{-1} Msun/yr. However, thermal feedback from SN generates turbulence and destroys molecules within the pristine gas, and only the polluted material, incorporated via gravitational enrichment, can continue to cool by atomic metal fine-structure transitions on time scales short enough to end the initial popIII regime within less than 10^8 yr.

The mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry / wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra-deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the Atlas-3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec-2 in the g-band. Initial results for two galaxies are presented here.

We present results on the heating of the inter-cluster medium (ICM) by gravitational potential energy from in-falling satellites. We calculate the available excess energy of baryons once they are stripped from their satellite and added to the ICM of the hosting environment. this excess energy is a strong function of environment and we find that it can exceed the contribution from AGNs or supernovae (SN) by up to two orders of magnitude in the densest environments/haloes. Cooling by radiative losses is in general fully compensated by gravitational heating in massive groups and clusters with hot gas temperature > 1 keV. The reason for the strong environment dependence is the continued infall of substructure onto dense environments during their formation in contrast to field-like environments. We show that gravitational heating is able to reduce the number of too luminous galaxies in models and to produce model luminosity functions in agreement with observations.

Current semi-analytic models (SAMs) of galaxy formation over-predict the fraction of passive small late-type satellite galaxies in dense environments by a factor of two to three. We hypothesize that this is due to inaccurate prescriptions on cold gas evolution. In the hope of solving this problem we apply detailed prescriptions on the evolution of diffuse hot gases in satellites and on stellar mass loss, both of which are critical to model cold gas evolution. We replace the conventional shock-heating motivated instant stripping with a realistic gradual prescription based on ram pressure and tidal stripping. We also carefully consider stellar mass loss in our model. When both mechanisms are included, the fraction of passive late types matches the data much more closely. The satellite over-quenching problem is still present in small galaxies in massive haloes, however. In terms of the detectable residual star formation rates, gradual diffuse gas stripping appears to be much more important than stellar mass loss in our model. The implications of these results and other possibilities, such as redshift-dependent merging geometry and tidal disruption, are also discussed.

Motivated by recent theoretical work suggesting that a substantial fraction of Population (Pop) III stars may have had masses low enough for them to survive to the present day, we consider the role that the accretion of metal-enriched gas may have had in altering their surface composition, thereby disguising them as Pop II stars. We demonstrate that if weak, Solar-like winds are launched from low-mass Pop III stars formed in the progenitors of the dark matter halo of the Galaxy, then such stars are likely to avoid significant enrichment via accretion of material from the interstellar medium. We find that at early times accretion is easily prevented if the stars are ejected from the central regions of the haloes in which they form, either by dynamical interactions with more massive Pop III stars, or by violent relaxation during halo mergers. While accretion may still take place during passage through sufficiently dense molecular clouds at later times, we find that the probability of such a passage is generally low (< 0.1), assuming that stars have velocities of order the maximum circular velocity of their host haloes and accounting for the orbital decay of merging haloes. In turn, due to the higher gas density required for accretion onto stars with higher velocities, we find an even lower probability of accretion (~ 0.01) for the subset of Pop III stars formed at z > 10, which are more quickly incorporated into massive haloes than stars formed at lower redshift. While there is no a priori reason to assume that low-mass Pop III stars do not have Solar-like winds, without them surface enrichment via accretion is likely to be inevitable. We briefly discuss the implications that our results hold for stellar archaeology.

The Atlas3D project is a multi-wavelength survey combined with a theoretical modeling effort. The observations span from the radio to the millimeter and optical, and provide multi-colour imaging, two-dimensional kinematics of the atomic (HI), molecular (CO) and ionized gas (Hbeta, [OIII] and [NI]), together with the kinematics and population of the stars (Hbeta, Fe5015 and Mgb), for a carefully selected, volume-limited (1.16*10^5 Mpc^3) sample of 260 early-type (elliptical E and lenticular S0) galaxies (ETGs). The models include semi-analytic, N-body binary mergers and cosmological simulations of galaxy formation. Here we present the science goals for the project and introduce the galaxy sample and the selection criteria. The sample consists of nearby (D<42 Mpc) morphologically-selected ETGs extracted from a parent sample of 871 galaxies (8% E, 22% S0 and 70% spirals) brighter than M_K<-21.5 mag (stellar mass M_Star>6*10^9 M_Sun). We analyze possible selection biases and we conclude that the parent sample is essentially complete and statistically representative of the nearby galaxy population. We present the size-luminosity relation for the spirals and ETGs and show that the ETGs in the Atlas3D sample define a tight red sequence in a colour-magnitude diagram, with few objects in the transition from the blue cloud. We describe the strategy of the SAURON integral-field observations and the extraction of the stellar kinematics with the pPXF method. We give an overview of the characteristics of the other main datasets already available for our sample and of the ongoing modelling projects.

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

We study the specific star formation rate (SSFR) and its evolution at $z\gtsim 4$, in models of galaxy formation, where the star formation is driven by cold accretion flows. We show that constant star formation and feedback efficiencies cannot reproduce the observed trend of SSFR with stellar mass and its observed lack of evolution at $z>4$. Model galaxies with $\log(M_*) \ltsim 9.5$ M$_{\odot}$ show systematically lower specific star formation rates by orders of magnitudes, while massive galaxies with M$_* \gtsim 5 \times 10^{10}$ M$_{\odot}$ have up to an order of magnitude larger SSFRs, compared to recent observations by Stark et al.. To recover these observations we apply an empirical star formation efficiency in galaxies that scales with the host halo velocity dispersion as $\propto 1/\sigma^3$ during galaxy mergers. We find that this modification needs to be of stochastic nature to reproduce the observations, i.e. only applied during mergers and not during accretion driven star formation phases. Our choice of star formation efficiency during mergers allows us to capture both, the boost in star formation at low masses and the quenching at high masses, and at the same time produce a constant SSFR-stellar mass relation at $z\gtsim 4$ under the assumption that most of the observed galaxies are in a merger triggered star formation phase. Our results suggest that observed high-z low mass galaxies with high SSFRs are likely to be frequently interacting systems, which experienced bursts in their star formation rate and efficiency (mode 1), in contrast to low redshift $z \ltsim 3$ galaxies which are cold accretion-regulated star forming systems with lower star formation efficiencies (mode 2).

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

Observational and theoretical studies of extragalactic radio sources have suggested that an inhomogeneous environment may be responsible for observed arm length asymmetries of jets and the properties of extended emission line regions in high redshift radio galaxies. We perform 3D hydrodynamic simulations of the interaction of a powerful extragalactic bipolar jet with a disc-shaped clumpy interstellar medium of log-normal density distribution and analyze the asymmetry. Furthermore, we compute the relation between jet asymmetry and the ISM properties by means of Monte Carlo simulations based on a 1D propagation model for the jet through the dense medium. We find that the properties of the ISM can be related to a probability distribution of jet arm length asymmetries: Disc density and height are found to have the largest effect on the asymmetry for realistic parameter ranges, while the Fourier energy spectrum of the ISM and turbulent Mach number only have a smaller effect. The hydrodynamic simulations show that asymmetries generally may be even larger than expected from the 1D model due to the complex interaction of the jet and its bow shock with gaseous clumps, which goes much beyond simple energy disposal. From our results, observed asymmetries of medium-sized local radio galaxies may be explained by gas masses of 10^9 to 10^10 solar masses in massive elliptical galaxies. Furthermore, the simulations provide a theoretical basis for the observed correlation that emission line nebulae are generally found to be brighter on the side of the shorter lobe in high redshift radio galaxies (McCarthy et al. 1991). This interaction of jets with the cold gas phase suggests that star formation in evolving high redshift galaxies may be affected considerably by jet activity.

One possible scenario for the formation of massive black holes (BHs) in the early Universe is from the direct collapse of primordial gas in atomic-cooling dark matter haloes in which the gas is unable to cool efficiently via molecular transitions. We study the formation of such BHs, as well as the accretion of gas onto these objects and the high energy radiation emitted in the accretion process, by carrying out cosmological radiation hydrodynamics simulations. In the absence of radiative feedback, we find an upper limit to the accretion rate onto the central object which forms from the initial collapse of hot (~ 10^4 K) gas of the order of 0.1 MSun per year. This is high enough for the formation of a supermassive star, the immediate precursor of a BH, with a mass of the order of 10^5 MSun. Assuming that a fraction of this mass goes into a BH, we track the subsequent accretion of gas onto the BH self-consistently with the high energy radiation emitted from the accretion disk. Using a ray-tracing algorithm to follow the propagation of ionizing radiation, we model in detail the evolution of the photoionized region which forms around the accreting BH. We find that BHs with masses of the order of 10^4 MSun initially accrete at close to the Eddington limit, but that the accretion rate drops to of order 1 percent of the Eddington limit after ~ 10^6 yr, due to the expansion of the gas near the BH in response to strong photoheating and radiation pressure. One signature of the accretion of gas onto BHs formed by direct collapse, as opposed to massive Pop III star formation, is an extremely high ratio of the luminosity emitted in He II 1640 to that emitted in H_alpha (or Ly_alpha); this could be detected by the James Webb Space Telescope. Finally, we briefly discuss implications for the coevolution of BHs and their host galaxies.

The star formation history of nearby early-type galaxies is investigated via numerical modelling. Idealized hydrodynamical N-body simulations with a star formation prescription are used to study the minor merger process between a giant galaxy (host) and a less massive spiral galaxy (satellite) with reasonable assumptions for the ages and metallicities of the merger progenitors. We find that the evolution of the star formation rate is extended over several dynamical times and shows peaks which correspond to pericentre passages of the satellite. The newly formed stars are mainly located in the central part of the satellite remnant while the older stars of the initial disk are deposited at larger radii in shell-like structures. After the final plunge of the satellite, star formation in the central part of the remnant can continue for several Gyrs depending on the star formation efficiency. Although the mass fraction in new stars is small, we find that the half-mass radius differs from the half-light radius in the V and H bands. Moreover synthetic 2D images in J, H, NUV, Hb and V bands, using the characteristic filters of the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST), reveal that residual star formation induced by gas-rich minor mergers can be clearly observed during and after the final plunge, especially in the NUV band, for interacting systems at (z<0.023) over moderate numbers of orbits (~2 orbits correspond to typical exposure times of ~3600 sec). This suggests that WFC3 has the potential to resolve these substructures, characterize plausible past merger episodes, and give clues to the formation of early-type galaxies.

Nuclear Stellar Disks (NSDs), of a few tens to hundreds of parsec across, are a common and yet poorly studied feature of early-type galaxies. Still, such small disks represent a powerful tool to constrain the assembling history of galaxies, since they can be used to trace to the epoch when galaxies experienced their last major merger event. By studying the fraction and stellar age of NSDs it is thus possible to test the predictions for the assembly history of early-type galaxies according the current hierarchical paradigm for galaxy formation. In this paper we have produced the most comprehensive census of NSDs in nearby early-type galaxies by searching for such disks in objects within 100 Mpc and by using archival images from the Hubble Space Telescope. We found that NSDs are present in approximately 20% of early-type galaxies, and that the fraction of galaxies with NSDs does not depend on their Hubble type nor on their galactic environment, whereas the incidence of NSDs appears to decline in the most massive systems. Furthermore, we have separated the light contribution of twelve such disks from that of their surrounding stellar bulge in order to extract their physical properties. This doubles the number of decomposed NSDs and although the derived values for their central surface brightness and scale-length are consistent with previous studies they also give a hint of possible different characteristics due to different formation scenario between nuclear disks and other kinds of large galactic disks.

We present results on the evolution of the intrinsic scatter of black hole masses considering different implementations of a model in which black holes only grow via mergers. We demonstrate how merger driven growth affects the correlations between black hole mass and host bulge mass. The simple case of an initially log-normal distributed scatter in black hole and bulge masses combined with random merging within the galaxy population results in a decreasing scatter with merging generation/number as predicted by the Central-limit theorem. In general we find that the decrease in scatter {\sigma} is well approximated by {\sigma}merg(m) = {\sigma}ini \times (m + 1)^(-a/2) with a = 0.42 for a range of mean number of mergers m < 50. For a large mean number of mergers (m > 100) we find a convergence to a = 0.61. This is valid for a wide range of different initial distributions, refill-scenarios or merger mass-ratios. Growth scenarios based on halo merger trees of a (100 Mpc)^3 dark matter LambdaCDM-simulation show a similar behaviour with a scatter decrease of a = 0.30 with typical number of mergers m < 50 consistent with random merging (best matching model: a = 0.34). Assuming a present day scatter of 0.3 dex in black hole mass and a mean number of mergers not exceeding m = 50 our results imply a scatter of 0.6 dex at z = 3 and thus a possible scenario in which overmassive (and undermassive) black holes at high redshift are a consequence of a larger intrinsic scatter in black hole mass. A simple toy model connecting the growth of black holes to the growth of LambdaCDM dark matter halos via mergers, neglecting any contribution from accretion, yields a consistent M\cdot -MBulge relation at z = 0 - if we assume the correct initial relation.

We study resolution effects in numerical simulations of gas-rich and gas-poor major mergers, and show that the formation of slowly-rotating elliptical galaxies often requires a resolution that is beyond the present-day standards to be properly modelled. Our sample of equal-mass merger models encompasses various masses and spatial resolutions, ranging from about 200pc and 10^5 particles per component, typical of some recently published major merger simulations, to up to 32pc and 10^3 M_sun in simulations using 2.4 x 10^7 collisionless particles and 1.2 x 10^7 gas particles, among the highest resolutions reached so far for gas-rich major merger of massive disc galaxies. We find that the formation of fast-rotating early-type galaxies, that are flattened by a significant residual rotation, is overall correctly reproduced at all such resolutions. However, the formation of slow-rotating early-type galaxies, which have a low residual angular momentum and are supported mostly by anisotropic velocity dispersions, is strongly resolution-dependent. The evacuation of angular momentum from the main stellar body is largely missed at standard resolution, and systems that should be slow rotators are then found to be fast rotators. The effect is most important for gas-rich mergers, but is also witnessed in mergers with an absent or modest gas component. The effect is robust with respect to our initial conditions and interaction orbits, and originates in the physical treatment of the relaxation process during the coalescence of the galaxies. Our findings show that a high-enough resolution is required to accurately model the global properties of merger remnants and the evolution of their angular momentum. The role of gas-rich mergers of spiral galaxies in the formation of slow-rotating ellipticals may therefore have been underestimated.

We present results from the first cosmological simulations which study the onset of primordial, metal-free (population III), cosmic star formation and the transition to the present-day, metal-rich star formation (population II-I), including molecular (H$_2$, HD, etc.) evolution, tracing the injection of metals by supernov{\ae} into the surrounding intergalactic medium and following the change in the initial stellar mass function (IMF) according to the metallicity of the corresponding stellar population. Our investigation addresses the role of a wide variety of parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and the associated transition in the star formation mode. All simulations present common trends. Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, inducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a result of the rapid pollution within high-density regions due to the first SN/pair-instability SN, local metallicity is quickly boosted above the critical metallicity for the transition. The population III regime lasts for a very short period during the first stages of star formation ($\sim 10^7\,\rm yr$), and its average contribution to the total star formation rate density drops rapidly below $\sim 10^{-3}-10^{-2}$.

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

We present a preliminary analysis of the HI properties of early-type galaxies in the ATLAS3D sample. Using WSRT data for ~100 galaxies outside the Virgo cluster and data from the Alfalfa project for galaxies inside Virgo, we discuss the dependence of HI properties on environment. We detect HI in about half of the galaxies outside Virgo. For these systems, the HI morphology and kinematics change as a function of environment, going from regular, rotating systems around isolated galaxies to progressively more disturbed structures for galaxies with neighbours or in groups. In denser environment, inside Virgo, nearly none of the galaxies contains HI. We discuss future work in this field which will be enabled by next-generation, pre-SKA radio instruments. We present a simulated Apertif HI observation of an ATLAS3D early-type galaxy, showing how its appearance and detection level vary as a function of redshift.

We study the origin of the scaling relations of early-type galaxies (ETGs) by constructing detailed models of the stellar dynamics for the K-band selected, volume-limited ATLAS3D sample of 263 nearby ETGs, spanning a large range of masses and stellar velocity dispersions (60 < sigma < 350 km/s).

We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3e7 galaxies, obtained by post-processing the virtual galaxy catalog produced by (De Lucia et al. 2007) on the Millennium Simulation of cosmic structure (Springel et al. 2005). Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local HI- and H2-mass functions, the CO-luminosity function, the cold gas mass--diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local Universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift.

We present results of a semi-analytic model (SAM) that includes cold accretion and a porosity-based prescription for star formation. We can recover the puzzling observational results of low $V/\sigma$ seen in various massive disk or disk-like galaxies, if we allow 18 % of the accretion energy from cold flows to drive turbulence in gaseous disks at $z=2$. The increase of gas mass through cold flows is by itself not sufficient to increase the star formation rate sufficiently to recover the number density of $\dot{M}_*>120$ M$_{\odot}$ yr$^{-1}$ galaxies in our model. In addition, it is necessary to increase the star formation efficiency. This can be achieved naturally in the porosity model, where star formation efficiency scales $\propto \sigma$, which scales as cloud velocity dispersion. As cold accretion is the main driver for gas velocity dispersion in our model, star formation efficiency parallels cold accretion rates, and allows fast conversion into stars. At $z\sim 2$, we find a space density $~10^{-4}$ Mpc$^{-3}$ in star-forming galaxies with $\dot{M}_*>120$ M$_{\odot}$ yr$^{-1}$, in better agreement than earlier estimates form SAMs. However, the fundamental relation between $\dot{M}_*$ and $M_*$ is still offset from the observed relation, indicating the need for possibly more efficient star formation at high z perhaps associated with a role for AGN triggering.

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of $ \sim 6 \times 10^{-5}$ Mpc$^{-3}$ Gyr$^{-1}$ at $z \leq 1$ and a steep decline at larger z. Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early-type and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of $M_{*,c} \sim 6.3 \times 10^{10}$ M$_{\odot}$ by mergers only. As a consequence of this transition, we find that around $M_{*,c}$, the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggest that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z. We however find that only $10% -20%$ of galaxies more massive than $M_{*,c}$ experience dry major mergers within their last Gyr at any given redshift $z \le 1$.

We use the bulge Sersic index n and bulge-to-total ratio (B/T) to explore the fundamental question of how bulges form. We perform 2D bulge-disk-bar decomposition on H-band images of 143 bright, high stellar mass (>1.0e10 solar masses) low-to-moderately inclined (i<70 degrees) spirals. Our results are: (1) Our H-band bar fraction (~58%) is consistent with that from ellipse fits. (2) 70% of the stellar mass is in disks, 10% in bars, and 20% in bulges. (3) A large fraction (~69%) of bright spirals have B/T<0.2, and ~76% have low n<2 bulges. These bulges exist in barred and unbarred galaxies across a wide range of Hubble types. (4) About 65% (68%) of bright spirals with n<2 (B/T<0.2) bulges host bars, suggesting a possible link between bars and bulges. (5) We compare the results with predictions from a set of LCDM models. In the models, a high mass spiral can have a bulge with a present-day low B/T<0.2 only if it did not undergo a major merger since z<2. The predicted fraction (~1.6%) of high mass spirals, which have undergone a major merger since z<4 and host a bulge with a present-day low B/T<0.2, is a factor of over thirty smaller than the observed fraction (~66%) of high mass spirals with B/T<0.2. Thus, contrary to common perception, bulges built via major mergers since z<4 seriously fail to account for the bulges present in ~66% of high mass spirals. Most of these present-day low B/T<0.2 bulges are likely to have been built by a combination of minor mergers and/or secular processes since z<4.

We report the discovery of a population of nearby, blue early-type galaxies with high star formation rates (0.5 < SFR < 50 Msun/yr). They are identified by their visual morphology as provided by Galaxy Zoo for SDSS DR6 and their u-r colour. We select a volume-limited sample in the redshift range 0.02 < z < 0.05, corresponding to luminosities of approximately L* and above, and with u-r colours significantly bluer than the red sequence. We confirm the early-type morphology of the objects in this sample and investigate their environmental dependence and star formation properties. Blue early-type galaxies tend to live in lower-density environments than `normal' red sequence early-types and make up 5.7 +/-0.4% of the low-redshift early-type galaxy population. We find that such blue early-type galaxies are virtually absent at high velocity dispersions above 200 km/s. Our analysis uses emission line diganostic diagrams and we find that ~25% of them are actively starforming, while another ~25% host both star formation and an AGN. Another ~12% are AGN. The remaining 38% show no strong emission lines. When present and uncontaminated by an AGN contribution, the star formation is generally intense. We consider star formation rates derived from Halpha, u-band and infrared luminosities, and radial colour profiles, and conclude that the star formation is spatially extended. Of those objects that are not currently undergoing star formation must have ceased doing so recently in order to account for their blue optical colours. The gas phase metallicity of the actively starforming blue early-types galaxies is supersolar in all cases. We discuss the place of these objects in the context of galaxy formation. A catalogue of all 204 blue early-type galaxies in our sample, including star formation rates and emission line classification, is provided.

We discuss the transitions of galaxy morphologies within the CDM paradigm under the assumption of bulge formation in mergers and disk growth via cooling of gas and subsequent star formation. Based on the relative importance of these two competing processes it is possible to make predictions on the expected morphological mix of galaxies. In particular we here discuss the generation of massive disk galaxies with low bulge-to-total mass ratios. Our results indicate that it is difficult to generate enough massive disk galaxies with B/T $< 0.2$ via major mergers and subsequent disk re-growth, if during the major merger progenitor disks get disrupted completely. On average low B/T galaxies must have had there last major merger at $z \ge 2$. The main limiting factor is the ability to re-grow massive disks at late times after the last major merger of a galaxy. Taking into account the contribution from minor mergers ($4 \ge M_1/M_2$, $M_1 \ge M_2$) to the formation of bulges, we recover the right fraction of massive low B/T disk galaxies, indicating that minor mergers play an important role in the formation of massive low B/T disk galaxies.

We discuss how the conditions at high redshift differ from those at low redshift, and what the impact is on the galaxy population. We focus in particular on the role of gaseous dissipation and its impact on sustaining high star formation rates as well as on driving star-bursts in mergers. Gas accretion onto galaxies at high redshifts occurs on a halo dynamical time allowing for very efficiently sustained star formation. In addition cold accretion flows are able to drive turbulence in high redshift disks at the level observed if at least 20% of the accretion energy is converted into random motion in the gaseous disk. In general we find that the fraction of gas involved in galaxy mergers is a strong function of time and increases with redshift. A model combining the role of dissipation during mergers and continued infall of satellite galaxies allows to reproduce the observed size-evolution of early-type galaxies with redshift. Furthermore we investigate how the evolution of the faint-end of the luminosity function can be explained in terms of the evolution of the underlying dark matter evolution.

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the HST/ACS images of the GOODS North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin2 with stellar masses above 3E10 Msun in the redshift range 0.4<z<1.2. Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems (>3E11 Msun) do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dFGRS, we find that the number density of massive early-type galaxies is consistent with no evolution between z=1.2 and 0, i.e. over an epoch spanning more than half of the current age of the Universe. Massive galaxies show very homogeneous **intrinsic** colour distributions, featuring red cores with small scatter. The distribution of half-light radii -- when compared to z=0 and z>1 samples -- is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers. However, in a more speculative fashion, the observations can also be interpreted as weak or even no evolution in comoving number density **and size** between 0.4<z<1.2, thus pushing major mergers of the most massive galaxies towards lower redshifts.

We demonstrate that the large scatter in the ultra-violet (UV) colours of intermediate-mass early-type galaxies in the local Universe and the inferred low-level recent star formation in these objects can be reproduced by minor mergers in the standard LCDM cosmology. Numerical simulations of mergers with mass ratios less than or equal to 1:4, with reasonable assumptions for the ages, metallicities and dust properties of the merger progenitors, produce good agreement to the observed UV colours of the early-type population, if the infalling satellites are assumed to have (cold) gas fractions greater than 20%. Early-types that satisfy (NUV-r) < 3.8 are likely to have experienced mergers with mass ratios between 1:4 and 1:6 within the last ~1.5 Gyrs, while those that satisfy 3.8<(NUV-r)<5.5 are consistent with either recent mergers with mass ratios < 1:6 or mergers with higher mass ratios that occurred more than ~1.5 Gyrs in the past. We demonstrate that the early-type colour-magnitude relations and colour distributions in both the UV and optical spectral ranges are consistent with the expected frequency of minor merging activity in the standard LCDM cosmology at low redshift. We present a strong plausibility argument for minor mergers to be the principal mechanism behind the large UV scatter and associated low-level recent star formation observed in early-type galaxies in the nearby Universe.

We present results of an attempt to include more detailed gas physics motivated from hydrodynamical simulations within semi-analytic models (SAM) of galaxy formation, focusing on the role that environmental effects play. The main difference to previous SAMs is that we include 'gravitational' heating of the intra-cluster medium (ICM) by the net surplus of gravitational potential energy released from gas that has been stripped from infalling satellites. Gravitational heating appears to be an efficient heating source able to prevent cooling in environments corresponding to dark matter halos more massive than $\sim 10^{13} $M$_{\odot}$. The energy release by gravitational heating can match that by AGN-feedback in massive galaxies and can exceed it in the most massive ones. However, there is a fundamental difference in the way the two processes operate. Gravitational heating becomes important at late times, when the peak activity of AGNs is already over, and it is very mass dependent. This mass dependency and time behaviour gives the right trend to recover down-sizing in the star-formation rate of massive galaxies. Abridged...

We combine deep UBVRIzJK photometry from the MUSYC survey with redshifts from the COMBO-17 survey to study the rest-frame ultraviolet (UV) properties of 674 high-redshift (0.5<z<1) early-type galaxies, drawn from the Extended Chandra Deep Field South (E-CDFS). Galaxy morphologies are determined through visual inspection of Hubble Space Telescope (HST) images taken from the GEMS survey. We harness the sensitivity of the UV to young (<1 Gyrs old) stars to quantify the recent star formation history of the early-type population. We find compelling evidence that early-types of all luminosities form stars over the lifetime of the Universe, although the bulk of their star formation is already complete at high redshift. Luminous (-23<M(V)<-20.5) early-types form 10-15 percent of their mass after z=1, while their less luminous (M(V)>-20.5) counterparts form 30-60 percent of their mass in the same redshift range.

We have studied ~2100 early-type galaxies in the SDSS DR3 which have been detected by the GALEX Medium Imaging Survey (MIS), in the redshift range 0 < z < 0.11. Combining GALEX UV photometry with corollary optical data from the SDSS, we find that, at a 95 percent confidence level, at least ~30 percent of galaxies in this sample have UV to optical colours consistent with some recent star formation within the last Gyr. In particular, galaxies with a NUV - r colour less than 5.5 are very likely to have experienced such recent star formation, taking into account the possibility of a contribution to NUV flux from the UV upturn phenomenon. We find quantitative agreement between the observations and the predictions of a semi-analytical LCDM hierarchical merger model and deduce that early-type galaxies in the redshift range 0 < z < 0.11 have ~1 to 3 percent of their stellar mass in stars less than 1 Gyr old. The average age of this recently formed population is ~300 to 500 Myrs. We also find that monolithically evolving galaxies, where recent star formation can be driven solely by recycled gas from stellar mass loss, cannot exhibit the blue colours (NUV - r < 5.5) seen in a significant fraction (~30 percent) of our observed sample.

We combine deep UBVRIzJK photometry from the Multiwavelength Survey by Yale-Chile (MUSYC) with redshifts from the COMBO-17 survey to perform a large-scale study of the rest-frame ultraviolet (UV) properties of 674 high-redshift (0.5<z<1) early-type galaxies, drawn from the Extended Chandra Deep Field South (E-CDFS). Galaxy morphologies are determined through visual inspection of Hubble Space Telescope (HST) images taken from the GEMS survey. We harness the sensitivity of the UV to young (<1 Gyr old) stars to quantify the recent star formation history of early-type galaxies across a range of luminosities (-23.5 < M(V) < -18). Comparisons to simple stellar populations forming at high redshift indicate that only ~1.1 percent of early-types in this sample are consistent with purely passive ageing since z=2. Parametrising the recent star formation (RSF) in terms of the mass fraction of stars less than a Gyr old, we find that the early-type population as a whole shows a typical RSF between 5 and 13% in the redshift range 0.5<z<1. Early-types on the UV red sequence show RSF values less than 5% while the reddest early-types are virtually quiescent with RSF values of ~1%. We find compelling evidence that early-types of all luminosities form stars over the lifetime of the Universe, although the bulk of their star formation is already complete at high redshift. This tail-end of star formation is measurable and not negligible, with luminous (-23<M(V)<-20.5) early-types potentially forming 10-15% of their mass since z=1, with their less luminous (M(V)>-20.5) counterparts potentially forming 30-60 percent of their mass in the same redshift range. (abridged)

We investigate the evolution of the faint-end slope of the luminosity function, $\alpha$, using semi-analytical modeling of galaxy formation. In agreement with observations, we find that the slope can be fitted well by $\alpha (z) =a+b z$, with a=-1.13 and b=-0.1. The main driver for the evolution in $\alpha$ is the evolution in the underlying dark matter mass function. Sub-L_* galaxies reside in dark matter halos that occupy a different part of the mass function. At high redshifts, this part of the mass function is steeper than at low redshifts and hence $\alpha$ is steeper. Supernova feedback in general causes the same relative flattening with respect to the dark matter mass function. The faint-end slope at low redshifts is dominated by field galaxies and at high redshifts by cluster galaxies. The evolution of $\alpha(z)$ in each of these environments is different, with field galaxies having a slope b=-0.14 and cluster galaxies b=-0.05. The transition from cluster-dominated to field-dominated faint-end slope occurs roughly at a redshift $z_* \sim 2$, and suggests that a single linear fit to the overall evolution of $\alpha(z)$ might not be appropriate. Furthermore, this result indicates that tidal disruption of dwarf galaxies in clusters cannot play a significant role in explaining the evolution of $\alpha(z)$ at z< z_*. In addition we find that different star formation efficiencies a_* in the Schmidt-Kennicutt-law and supernovae-feedback efficiencies $\epsilon$ generally do not strongly influence the evolution of $\alpha(z)$.

We investigate the hierarchical build-up of stars in bulges within the standard $\Lambda$-cold dark matter scenario. By separating the population into stars born during starbursts that accompany the formation of spheroids in major mergers ({\it starburst} component), and stars that are previously formed in discs of progenitor galaxies ({\it quiescent} component) and added to the spheroid by dynamical interaction. Our results are summarised as follows: bulges that form early have larger starburst fraction and hence should be smaller than their counter parts that form later. The quiescent fraction in bulges is an increasing function of bulge mass, becoming constant at $M_{\rm{q}}/M_{\rm{bul}} \sim 0.8$, mainly due to the infall of satellite galaxies that contribute disc stars to the bulge. Minor mergers are an order of magnitude more frequent than major mergers and must play a significant role in the evolution of bulges. Above the critical mass $M_{\rm{c}}\sim 3 \times 10^{10}$ M$_{\odot}$ most of the stars in the universe are in spheroids, which at high redshift are exclusively elliptical galaxies and at low redshifts partly bulges. Due to the enhanced evolution of galaxies ending up in high density environments, the starburst fraction and the surface mass densities of bulges below $M_{\rm{c}}$ should be enhanced with respect to field galaxies. Dissipation during the formation of massive bulges in present day early-type spirals is less important than for the formation of present day elliptical galaxies of the same mass thereby explaining the possible difference in phase-space densities between spiral galaxies and elliptical galaxies.

Detailed high-resolution observations of the innermost regions of nearby galaxies have revealed the presence of supermassive black holes1. These black holes may interact with their host galaxies by means of 'feedback' in the form of energy and material jets; this feedback affects the evolution of the host and gives rise to observed relations between the black hole and the host. Here we report observations of the ultraviolet emissions of massive early-type galaxies. We derive an empirical relation for a critical black-hole mass (as a function of velocity dispersion) above which the outflows from these black holes suppress star formation in their hosts by heating and expelling all available cold gas. Supermassive black holes are negligible in mass compared to their hosts but nevertheless seem to play a critical role in the star formation history of galaxies.

We use semi-analytical modelling of galaxy formation to predict the redshift-size-evolution of elliptical galaxies. Using a simple model in which relative sizes of elliptical galaxies of a given mass correlate with the fraction of stars formed in a star burst during a major merger event, we are able to reproduce the observed size redshif t evolution. The size evolution is a result of the amount of cold gas available during the major merger. Mergers at high redshifts are gas-rich and produce ellipticals with smaller sizes. In particular we find a power-law relation between the sizes at different redshifts, with the power-law index giving a measure of the relative amount of dissipation during the mergers that lead to the formation of an elliptical. The size evolution is found to be stronger for more massive galaxies as they involve more gas at high redshifts when they form compared to less massive ellipticals. Local ellipticals more massive than $5 \times 10^{11}$ M$_{\odot}$ will be approximately 4 times larger than their counterparts at $z=2$. Our results indicate that the scatter in the size of similar massive present day elliptical galaxies is a result of their formation epoch, with smaller ellipticals being formed earlier.

We investigate the stellar composition of bulges and elliptical galaxies as predicted by the CDM paradigm using semi-analytical modelling. We argue that spheroid stars are built up of two main components, {\it merger} and {\it quiescent}, according to the origin of the stars. The merger component is formed during major mergers by gas driven to the centre, while the quiescent component is formed in gaseous discs and added later to the spheroid during major mergers. Galaxies more massive than $M_C=3 \times 10^{10}$ M$_{\odot}$ have on average only a 15% merger component in their spheroids, while smaller galaxies can have up to 30%. The merger component increases with redshift due to mergers involving more gas. However we do not find mergers with gas fraction above $ \sim 40 %$ of the remnants mass. Generally the gas fraction is a decreasing function of the redshift at which the merger occurs and the mass of the remnant, with more massive remnants having smaller gas fraction and hence smaller merger components. This trend is independent of the environment of the galaxy with the only impact of the environment being that galaxies less massive than $M_C$ have slightly larger merger components in dense environments. The fraction of stars in bulges for galaxies more massive than $M_C$ is larger than 50%. We find that the majority of stars in galaxies more massive than $M_C$ reside within bulges and ellipticals independent of redshift and that the fraction increases with redshift. The most massive galaxies at each redshift are elliptical galaxies.

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disk galaxies that are too compact with respect to the observations. High-resolution N-body/SPH simulations in a cosmological context are presented including cold gas and dark matter. A halo with quiet merging activity since z~3.8 and with a high spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disk. We show that the gas and the dark matter have similar specific angular momenta until a merger event occurs at z~2 with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and falls quickly into the center. Dynamical friction plays a minor role,in contrast to previous claims. In fact, after this event a new extended disk begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disk galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered star burst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected as a result of these processes, late-type galaxies could form with a dominant extended disk component, resulting from late infall, a small bulge-to-disk ratio and a low baryon fraction, in agreement with observations.

We use \textit{GALEX} (Galaxy Evolution Explorer) near-UV (NUV) photometry of a sample of early-type galaxies selected in \textit{SDSS} (Sloan Digital Sky Survey) to study the UV color-magnitude relation (CMR). $NUV-r$ color is an excellent tracer of even small amounts ($\sim 1$% mass fraction) of recent ($\la 1$ Gyr) star formation and so the $NUV-r$ CMR allows us to study the effect of environment on the recent star formation history. We analyze a volume-limited sample of 839 visually-inspected early-type galaxies in the redshift range $0.05 < z < 0.10$ brighter than $M_{r}$ of -21.5 with any possible emission-line or radio-selected AGN removed to avoid contamination. We find that contamination by AGN candidates and late-type interlopers highly bias any study of recent star formation in early-type galaxies and that, after removing those, our lower limit to the fraction of massive early-type galaxies showing signs of recent star formation is roughly $30 \pm 3%$ This suggests that residual star formation is common even amongst the present day early-type galaxy population. We find that the fraction of UV-bright early-type galaxies is 25% higher in low-density environments. However, the density effect is clear only in the lowest density bin. The blue galaxy fraction for the subsample of the brightest early-type galaxies however shows a very strong density dependence, in the sense that the blue galaxy fraction is lower in a higher density region.

The luminosity dependence of kinematical and isophotal properties of elliptical galaxies is investigated using numerical simulations of galaxy merging, combined with semi-analytical models of hierarchical structure formation. Mergers of spiral galaxies as the only formation mechanism for elliptical galaxies can neither reproduce the kinematical and photometric properties of very massive elliptical galaxies nor the change from rotationally flattened disky to anisotropic boxy systems with increasing luminosity. We present numerical simulations showing that binary mergers of early-type galaxies open an additional channel for the formation of anisotropic, slowly rotating and boxy ellipticals. Including this channel in a semi-analytical model we can successfully reproduce the observed trend that more luminous giant ellipticals are more boxy and less flattened by rotation. This trend can be strengthened by suppressing residual gas infall and star formation for galaxies with stellar bulge masses $M_* \geq 3 \times 10^{10} M_{\odot}$. Hence we propose that mergers of early-type galaxies play an important role for the assembly of massive elliptical galaxies.

In order to specify cosmologically motivated initial conditions for major galaxy mergers (mass ratios $\leq$ 4:1) that are supposed to explain the formation of elliptical galaxies we study the orbital parameters of major mergers of cold dark matter halos using a high-resolution cosmological simulation. Almost half of all encounters are nearly parabolic with eccentricities $e \approx 1$ and no correlations between the halo spin planes or the orbital planes. The pericentric argument $\omega$ shows no correlation with the other orbital parameters and is distributed randomly. In addition we find that 50 % of typical pericenter distances are larger than half the halo's virial radii which is much larger than typically assumed in numerical simulations of galaxy mergers. In contrast to the usual assumption made in semi-analytic models of galaxy formation the circularities of major mergers are found to be not randomly distributed but to peak around a value of $\epsilon \approx 0.5$. Additionally all results are independent of the minimum progenitor mass and major merger definitions (i.e. mass ratios $\leq$ 4:1; 3:1; 2:1).

Using semi-analytical models of galaxy formation, the origin of boxy and disky elliptical galaxies is investigated. We find that the simple scenario, motivated by N-body simulations, in which the isophotal shape is only dependent on the mass ratio of the last major merger, is not able to reproduce the observation that the fraction of boxy and disky ellipticals depends on galaxy luminosity. The observations can however be reproduce with the following reasonable assumptions: (i) equal-mass mergers lead to boxy ellipticals and unequal-mass mergers produce disky ellipticals (as motivated by N-body simulations) (ii) major mergers between bulge-dominated galaxies result always in boxy ellipticals, independent of the mass ratio, (iii) merger remnants that subsequently accrete gas leading to a secondary stellar disk with more than $20 %$ the total stellar fraction are always disky. Boxy and disky ellipticals can be divided into two subclasses, depending on their formation history. Boxy ellipticals either formed by equal mass mergers of disk galaxies or by major mergers of early-type galaxies. We find that disky ellipticals are indicators of unequal-mass mergers or late gas infall. Disky ellipticals with high luminosities preferentially are 1-component systems that result from unequal mass mergers whereas low-luminosity disky ellipticals are more likely to harbor secondary disk components. In addition, the fraction of disky ellipticals with secular disk components should increase in regions with higher galaxy densities. Taking into account the conversion of cuspy cores into flat low-density cores by black hole merging we find that disky ellipticals should contain central density cusps whereas boxy ellipticals should in general be characterized by flat cores.

We use semi-analytical modelling techniques to investigate the progenitor morphologies of present day ellipticals. We find that, independent of the environment, the fraction of mergers of bulge dominated galaxies (early-types) increases with time. The last major merger of bright present day ellipticals with $M_{B} \lesssim -21$ is preferentially between bulge dominated galaxies, while those with $ M_{B} \sim -20$ have mainly experienced last major mergers between a bulge dominated and a disk dominated galaxy. Independent of specific model assumptions, more than 50% of present day elliptical in clusters with $M_B \lesssim -18$ had last major mergers which are not of spirals as usually expected within the standard merger scenario.

We use semi-analytical modelling of galaxy formation to study the redshift evolution of the galaxy merger fractions and merger ratesin Lambda CDM and quintessence (QCDM) cosmologies, their dependence on physical parameters as the environment, the merger timescale, the way major mergers are defined, and the minimum mass of objects taken into account. We find that for a given final halo mass the redshift dependence of the merger fraction F_mg and the resulting merger rate can be fitted well by a power law for redshifts z <= 1. The normalization F_{mg}(0) and the slope m depend on the final halo mass. For a given merger timescale t_merg and an assumed maximum mass ratio R_major for major mergers, F_mg(0) and m depend exponentially on each other. The slope m depends logarithmically on the ratio of the final halo mass and the minimum halo mass taken into account. In addition, the local normalization F_{mg}(0) increases for larger R_{major} while m decreases. We compare the predicted merger fractionwith recent observations and find that the model cannot reproduce both the merger index and the normalization at the same time. In general the model underestimates F_{mg}(0) and m by a factor of 2.