We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, but it also remains uncertain whether the low-redshift LyC leakers are true "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission-line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, we find that identifying low-redshift galaxies based on [SII] deficiencies does not seem to produce true analogues. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158\mu m}$ and [OIII]$_{\rm 88\mu m}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149\_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. \ion{H}{i}, \ion{Si}{ii}, \ion{Mg}{ii}, \ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.

The physical origin of low escape fractions of ionizing radiation derived from Lyman-break galaxies (LBGs) at $z\sim3$--$4$ is not well understood. We perform idealised disc galaxy simulations to understand how galactic properties such as metallicity and gas mass affect the escape of Lyman Continuum (LyC) photons using radiation-hydrodynamic simulations with strong stellar feedback. We find that the luminosity-weighted escape fraction from a metal-poor ($Z=0.002$) galaxy embedded in a halo of mass $M_h\simeq10^{11}\,M_\odot$ is $\left<f_{\rm esc}^{\rm 3D}\right>\simeq7\%$. Roughly half of the LyC photons are absorbed within scales of 100 pc, and the other half is absorbed in the ISM ($\lesssim 2$ kpc). When the metallicity of the gas is increased to $Z=0.02$, the escape fraction is significantly reduced to $\left<f_{\rm esc}^{\rm 3D}\right>\simeq1\%$ because young stars are enshrouded by their birth clouds for a longer time. In contrast, increasing the gas mass by a factor of 5 leads to $\left<f_{\rm esc}^{\rm 3D}\right>\simeq4\%$ because LyC photons are only moderately absorbed by the thicker disc. Our experiments suggest that high metallicity is likely more responsible for the low escape fractions observed in LBGs, supporting the scenario in which the escape fraction is decreasing with increasing halo mass. Finally, negligible correlation is observed between the escape fraction and surface density of star formation or galactic outflow rates.

Using cosmological, radiation-hydrodynamic simulations targeting a rare $\approx \, 2 \times 10^{12} \, \rm M_\odot$ halo at $z \, = \, 6$, we show that the number counts and internal properties of satellite galaxies within the massive halo are sensitively regulated by a combination of local stellar radiative feedback and strong tidal forces. Radiative feedback operates before the first supernova explosions erupt and results in less tightly-bound galaxies. Satellites are therefore more vulnerable to tidal stripping when they accrete onto the main progenitor and are tidally disrupted on a significantly shorter timescale. Consequently, the number of satellites with $M_{\rm \star} > 10^{7} \, \rm M_\odot$ within the parent system's virial radius drops by up to $60 \%$ with respect to an identical simulation performed without stellar radiative feedback. Radiative feedback also impacts the central galaxy, whose effective radius increases by a factor $\lesssim 3$ due to the presence of a more extended and diffuse stellar component. We suggest that the number of satellites in the vicinity of massive high-redshift galaxies is an indication of the strength of stellar radiative feedback and and can be anomalously low in the extreme cosmic environments of high-redshift quasars.

The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values.

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA).

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction.

Accurate models of radiative cooling are a fundamental ingredient of modern cosmological simulations. Without cooling, accreted baryons will not efficiently dissipate their energy and collapse to the centres of haloes to form stars. It is well established that local variations in the amplitude and shape of the spectral energy distribution of the radiation field can drastically alter the cooling rate. Here we introduce deepCool, deepHeat, and deepMetal: methods for accurately modelling the total cooling rates, total heating rates, and metal-line only cooling rates of irradiated gas using artificial neural networks. We train our algorithm on a high-resolution cosmological radiation hydrodynamics simulation and demonstrate that we can predict the cooling rate, as measured with the photoionisation code CLOUDY, under the influence of a local radiation field, to an accuracy of ~5%. Our method is computationally and memory efficient, making it suitable for deployment in state-of-the-art radiation hydrodynamics simulations. We show that the circumgalactic medium and diffuse gas surrounding the central regions of a galaxy are most affected by the interplay of radiation and gas, and that standard cooling functions that ignore the local radiation field can incorrectly predict the cooling rate by more than an order of magnitude, indicating that the baryon cycle in galaxies is affected by the influence of a local radiation field on the cooling rate.

Despite their ubiquity, the origin of cosmic magnetic fields remains unknown. Various mechanisms have been proposed for their existence including primordial fields generated by inflation, or amplification and injection by compact astrophysical objects. Separating the potential impact of each magnetogenesis scenario on the magnitude and orientation of the magnetic field and their impact on gas dynamics may give insight into the physics that magnetised our Universe. In this work, we demonstrate that because the induction equation and solenoidal constraint are linear with $B$, the contribution from different sources of magnetic field can be separated in cosmological magnetohydrodynamics simulations and their evolution and influence on the gas dynamics can be tracked. Exploiting this property, we develop a magnetic field tracer algorithm for cosmological simulations that can track the origin and evolution of different components of the magnetic field. We present a suite of cosmological magnetohydrodynamical RAMSES simulations that employ this algorithm where the primordial field strength is varied to determine the contributions of the primordial and supernovae-injected magnetic fields to the total magnetic energy as a function of time and spatial location. We find that, for our specific model, the supernova-injected fields rarely penetrate far from haloes, despite often dominating the total magnetic energy in the simulations. The magnetic energy density from the supernova-injected field scales with density with a power-law slope steeper than 4/3 and often dominates the total magnetic energy inside of haloes. However, the star formation rates in our simulations are not affected by the presence of magnetic fields, for the ranges of primordial field strengths examined. These simulations represent a first demonstration of the magnetic field tracer algorithm (abridged).

We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}<M<10^{10}{\rm M_{\odot}}h^{-1}$, stars with metallicity $10^{-3}Z_{\odot}<Z<10^{-1.5}Z_{\odot}$, and stars with age $3\,{\rm Myr} < t < 10 \, {\rm Myr}$ dominate reionization by both mass and volume. We show that the sources that reionize most of the volume of the Universe by $z=6$ are not necessarily the same sources that dominate the meta-galactic UV background at the same redshift. We further show that in our simulations, the contribution of each type of source to reionization is not uniform across different gas phases. The IGM, CGM, filaments, ISM, and rarefied supernova heated gas have all been photoionized by different classes of sources. Collisional ionisation contributes at both the lowest and highest densities. In the early stages of the formation of individual HII bubbles, reionization proceeds with the formation of concentric shells of gas ionised by different classes of sources, leading to large temperature variations as a function of galacto-centric radius. The temperature structure of individual HII bubbles may thus give insight into the star formation history of the galaxies acting as the first ionising sources. Our explorative simulations highlight how the complex nature of reionization can be better understood by using our photon tracer algorithm.

The environmental effect is commonly used to explain the excess of gas-poor galaxies in galaxy clusters. Meanwhile, the presence of gas-poor galaxies at cluster outskirts, where galaxies have not spent enough time to feel the cluster environmental effect, hints for the presence of pre-processing. Using cosmological hydrodynamic simulations on 16 clusters, we investigate the mechanisms of gas depletion of galaxies found inside clusters. The gas depletion mechanisms can be categorized into three channels based on where and when they took place. First, 34$\%$ of our galaxies are gas poor before entering clusters (`pre-processing'). They are mainly satellites that have undergone the environmental effect inside group halos. Second, 43$\%$ of the sample became quickly gas deficient in clusters before the first pericentric pass (`fast cluster processing'). Some of them were group satellites that are low in gas at the time of cluster entry compared to the galaxies directly coming from the field. Even the galaxies with large gas fractions take this channel if they fall into massive clusters ($> 10^{14.5}\, \rm M_{\odot}$) or approach cluster centers through radial orbits. Third, 24$\%$ of our sample retain gas even after their first pericentric pass (`slow cluster processing') as they fall into the less massive clusters and/or have circular orbits. The relative importance of each channel varies with a cluster's mass, while the exact degree of significance is subject to large uncertainties. Group pre-processing accounts for a third of the total gas depletion; but it also determines the gas fraction of galaxies at their cluster entry which in turn determines whether a galaxy should take the fast or the slow cluster processing.

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies.

Using the Horizon-AGN simulation data, we study the relative role of mergers and environmental effects in shaping the spin of early-type galaxies (ETGs) after $z \simeq 1$. We follow the spin evolution of 10,037 color-selected ETGs more massive than 10$^{10} \rm \, M_{\odot}$ that are divided into four groups: cluster centrals (3%), cluster satellites (33%), group centrals (5%), and field ETGs (59%). We find a strong mass dependence of the slow rotator fraction, $f_{\rm SR}$, and the mean spin of massive ETGs. Although we do not find a clear environmental dependence of $f_{\rm SR}$, a weak trend is seen in the mean value of spin parameter driven by the satellite ETGs as they gradually lose their spin as their environment becomes denser. Galaxy mergers appear to be the main cause of total spin changes in 94% of central ETGs of halos with $M_{vir} > 10^{12.5}\rm M_{\odot}$, but only 22% of satellite and field ETGs. We find that non-merger induced tidal perturbations better correlate with the galaxy spin-down in satellite ETGs than mergers. Given that the majority of ETGs are not central in dense environments, we conclude that non-merger tidal perturbation effects played a key role in the spin evolution of ETGs observed in the local ($z < 1$) universe.

The abundance of nitrogen in the interstellar medium is a powerful probe of star for- mation processes over cosmological timescales. Since nitrogen can be produced both in massive and intermediate-mass stars with metallicity-dependent yields, its evolution is challenging to model, as evidenced by the differences between theoretical predictions and observations. In this work we attempt to identify the sources of these discrepancies using a cosmic evolution model. To further complicate matters, there is considerable dispersion in the abundances from observations of DLAs at redshift 2 - 3. We study the evolution of nitrogen with a detailed chemical evolution model and find good agreement with observations, including the relative abundances of N/O and N/Si ratios. We find that the principal contribution of nitrogen comes from intermediate mass stars, with the exception of systems with the lowest N/H, where nitrogen production might possibly be dominated by massive stars. This last result could be strengthened if stellar rotation which is important at low metallicity can produce significant amounts of nitrogen. Moreover, these systems likely reside in host galaxies with stellar masses below 10**8.5 solar mass. We also study the origin of the observed dispersion in nitrogen abundances using the cosmological hydrodynamical simulations Horizon-AGN. We conclude that this dispersion can originate from two effects: differ- ence in the masses of the DLA host galaxies, and difference in the their position inside the galaxy.

Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lya photon resonantly scatters and imparts 10-300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass loading factors in our metal-poor dwarf system are estimated to be $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lya feedback.

We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion.

We present cosmological, radiation-hydrodynamics simulations of galaxy formation during the epoch of reionization in an effort towards modelling the interstellar medium (ISM) and interpreting ALMA observations. Simulations with and without stellar radiation are compared at large (Mpc), intermediate (tens of kpc), and small (sub kpc) scales. At large scales, the dense regions around galaxies reionize first before UV photons penetrate the voids; however, considerable amounts of neutral gas remain present within the haloes. The spatial distribution of neutral gas is highly dynamic and is anti-correlated with the presence of stars older than a few Myrs. For our specific feedback implementation, most of the metals remain inside the virial radii of haloes and they are proportionally distributed over the ionized and neutral medium by mass. For our most massive galaxy with ${\rm M_h}\sim10^{11}$M$_{\odot}$, the majority of the CII and OI mass are associated with cold neutral clumps. NII is more diffuse and arises in warmer gas while OIII arises in hotter gas with a higher ionization parameter, produced by photo-heating and supernovae. If smaller pockets of high metallicity gas exist in the ISM, the emission from these ions may be observable by ALMA while the low metallicity of the galaxy may cause these systems to fall below the local [CII]-SFR relation. The presence of dust can cause spatial offsets between UV/Ly$\alpha$ and [CII] emission as suggested by the recent observations of Maiolino et al. [OIII] may be spatially offset from both of these components since it arises from a different part of density-temperature phase-space.

We compare the predictions of Horizon-AGN, a hydro-dynamical cosmological simulation that uses an adaptive mesh refinement code, to observational data in the redshift range 0<z<6. We study the reproduction, by the simulation, of quantities that trace the aggregate stellar-mass growth of galaxies over cosmic time: luminosity and stellar-mass functions, the star formation main sequence, rest-frame UV-optical-near infrared colours and the cosmic star-formation history. We show that Horizon-AGN, which is not tuned to reproduce the local Universe, produces good overall agreement with these quantities, from the present day to the epoch when the Universe was 5% of its current age. By comparison to Horizon-noAGN, a twin simulation without AGN feedback, we quantify how feedback from black holes is likely to help shape galaxy stellar-mass growth in the redshift range 0<z<6, particularly in the most massive galaxies. Our results demonstrate that Horizon-AGN successfully captures the evolutionary trends of observed galaxies over the lifetime of the Universe, making it an excellent tool for studying the processes that drive galaxy evolution and making predictions for the next generation of galaxy surveys.

Reionisation in the early Universe is likely driven by dwarf galaxies. Using cosmological radiation-hydrodynamic simulations, we study star formation and the escape of Lyman continuum (LyC) photons from mini-haloes with $M_{\rm halo} \le 10^8\,M_\odot$. Our simulations include a new thermo-turbulent star formation model, non-equilibrium chemistry, and relevant stellar feedback processes (photoionisation by young massive stars, radiation pressure, and mechanical supernova explosions). We find that feedback reduces star formation very efficiently in mini-haloes, resulting in the stellar mass consistent with the slope and normalisation reported in Kimm \& Cen and the empirical stellar mass-to-halo mass relation derived in the local Universe. Because star formation is stochastic and dominated by a few gas clumps, the escape fraction in mini-haloes is generally determined by radiation feedback (heating due to photo-ionisation), rather than supernova explosions. We also find that the photon number-weighted mean escape fraction in mini-haloes is higher ($\sim20$-$40\%$) than that in atomic-cooling haloes, although the instantaneous fraction in individual haloes varies significantly. The escape fraction from Pop III stars is found to be significant ($\ge10\%$) only when the mass is greater than $\sim$100\,\msun. Based on simple analytic calculations, we show that LyC photons from mini-haloes are, despite their high escape fractions, of minor importance for reionisation due to inefficient star formation. We confirm previous claims that stars in atomic-cooling haloes with masses $10^8\,M_\odot\le M_{\rm halo} \le 10^{11}\,M_\odot$ are likely to be the most important source of reionisation.

We compare 5 sub-grid models for supernova (SN) feedback in adaptive mesh refinement (AMR) simulations of isolated dwarf and L-star disk galaxies with 20-40 pc resolution. The models are thermal dump, stochastic thermal, 'mechanical' (injecting energy or momentum depending on the resolution), kinetic, and delayed cooling feedback. We focus on the ability of each model to suppress star formation and generate outflows. Our highest-resolution runs marginally resolve the adiabatic phase of the feedback events, which correspond to 40 SN explosions, and the first three models yield nearly identical results, possibly indicating that kinetic and delayed cooling feedback converge to wrong results. At lower resolution all models differ, with thermal dump feedback becoming inefficient. Thermal dump, stochastic, and mechanical feedback generate multiphase outflows with mass loading factors $\beta \ll 1$, which is much lower than observed. For the case of stochastic feedback we compare to published SPH simulations, and find much lower outflow rates. Kinetic feedback yields fast, hot outflows with $\beta\sim 1$, but only if the wind is in effect hydrodynamically decoupled from the disk by using a large bubble radius. Delayed cooling generates cold, dense and slow winds with $\beta> 1$, but large amounts of gas occupy regions of temperature-density space with short cooling times. We conclude that either our resolution is too low to warrant physically motivated models for SN feedback, that feedback mechanisms other than SNe are important, or that other aspects of galaxy evolution, such as star formation, require better treatment.

We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{\rm halo} \sim4\times10^7 M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5 M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$\alpha$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

Simulations have indicated that most of the escaped Lyman continuum photons escape through a minority of solid angles with near complete transparency, with the remaining majority of the solid angles largely opaque, resulting in a very broad and skewed probability distribution function (PDF) of the escape fraction when viewed at different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy observed along a line of sight merely represents the properties of the interstellar medium along that line of sight, which may be an ill-representation of true escape fraction of the galaxy averaged over its full sky. Here we study how Lyman continuum photons escape from galaxies at $z=4-6$, utilizing high-resolution large-scale cosmological radiation-hydrodynamic simulations. We compute the PDF of the mean escape fraction ($\left<f_{\rm esc,1D}\right>$) averaged over mock observational samples, as a function of the sample size, compared to the true mean (had you an infinite sample size). We find that, when the sample size is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, an observational sample of (2,10,50) galaxies at $z=4$ would have have 2.5% probability of obtaining the sample mean lower than $\left<f_{\rm esc,1D}\right>=$(0.007%, 1.8%, 4.1%) and 2.5% probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our simulations suggest that at least $\sim$ 100 galaxies should be stacked in order to constrain the true escape fraction within 20% uncertainty.

Analysis is performed on ultra-high resolution large-scale cosmological radiation-hydrodynamic simulations to, for the first time, quantify the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high density ($\ge 10^4$cm$^{-3}$) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density ($\sim 10^{-2.5}$cm$^{-3}$) high temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature and metallicity of intervening absorbers, vary strongly from sightline to sightline. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium of the host galaxy based on a very small number of (usually one) sightlines would be precarious. Utilizing high-z GRBs to probe interstellar medium and intergalactic medium should be undertaken properly taking into consideration the physical diversities of the interstellar medium.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

The fraction of hydrogen ionizing photons escaping from galaxies into the intergalactic medium is a critical ingredient in the theory of reionization. We use two zoomed-in, high-resolution (4 pc), cosmological radiation hydrodynamic simulations with adaptive mesh refinement to investigate the impact of two physical mechanisms (supernova feedback and runaway OB stars) on the escape fraction (f_esc) at the epoch of reionization (z>7). We implement a new, physically motivated supernova feedback model that can approximate the Sedov solutions at all (from the free expansion to snowplow) stages. We find that there is a significant time delay of about ten million years between the peak of star formation and that of escape fraction, due to the time required for the build-up and subsequent destruction of the star-forming cloud by supernova feedback. Consequently, the photon number-weighted mean escape fraction for dwarf galaxies in halos of mass 10^8-10^10.5 Msun is found to be <fesc>~11%, although instantaneous values of f_esc>20% are common when star formation is strongly modulated by the supernova explosions. We find that the inclusion of runaway OB stars increases the mean escape fraction by 22% to <fesc>~14%. As supernovae resulting from runaway OB stars tend to occur in less dense environments, the feedback effect is enhanced and star formation is further suppressed in halos with Mvir>10^9 Msun in the simulation with runaway OB stars compared with the model without them. While both our models produce enough ionizing photons to maintain a fully ionized universe at z>7 as observed, a still higher amount of ionizing photons at z>9 appears necessary to accommodate the high observed electron optical depth inferred from cosmic microwave background observations.

Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range 5x10^8-3x10^10 Msun from a cosmological hydrodynamic simulation (resolved at 29pc/h) at z=7. In a companion paper (Kimm & Cen), the stellar mass and UV luminosity functions, and UV-optical and FUV-NUV colors are shown to be in good agreement with observations, if an SMC-type dust extinction curve is adopted. Here we make useful predictions, self-consistently, of the infrared properties of these z=7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS 70 micron band is in the range of 80-400pc proper for z=7 galaxies with L_FIR=10^{11.3-12}Lsun. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 micron, depending on the dust-to-metal ratio. (3) For star formation rate in the range 3-100 Msun/yr the median FIR to bolometric luminosity ratio is 60-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of -3.1 +- 0.4, instead of the usual exponential decline.

Hubble Space Telescope observations with Wide Field Camera 3/IR reveal that galaxies at z~7 have very blue ultraviolet (UV) colors, consistent with these systems being dominated by young stellar populations with moderate or little attenuation by dust. We investigate UV and optical properties of the high-z galaxies in the standard cold dark matter model using a high-resolution adaptive mesh refinement cosmological hydrodynamic simulation. For this purpose, we perform panchromatic three-dimensional dust radiative transfer calculations on 198 galaxies of stellar mass 5x10^8-3x10^{10} Msun with three parameters, the dust-to-metal ratio, the extinction curve, and the fraction of directly escaped light from stars (\fesc). Our stellar mass function is found to be in broad agreement with Gonzalez et al., independent of these parameters. We find that our heavily dust-attenuated galaxies (A_V~1.8) can also reasonably match modest UV-optical colors, blue UV slopes, as well as UV luminosity functions, provided that a significant fraction (~10%) of light directly escapes from them. The observed UV slope and scatter are better explained with an Small Magellanic Cloud-type extinction curve, whereas Milky Way-type dust predicts too blue UV colors due to the 2175A bump. We expect that upcoming observations by ALMA will be able to test this heavily obscured model.

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

We study stellar assembly and feedback from active galactic nuclei (AGN) around the epoch of peak star formation (1<z<2), by comparing hydrodynamic simulations to rest-frame UV-optical galaxy colours from the Wide Field Camera 3 (WFC3) Early-Release Science (ERS) Programme. Our Adaptive Mesh Refinement simulations include metal-dependent radiative cooling, star formation, kinetic outflows due to supernova explosions, and feedback from supermassive black holes. Our model assumes that when gas accretes onto black holes, a fraction of the energy is used to form either thermal winds or sub-relativistic momentum-imparting collimated jets, depending on the accretion rate. We find that the predicted rest-frame UV-optical colours of galaxies in the model that includes AGN feedback is in broad agreement with the observed colours of the WFC3 ERS sample at 1<z<2. The predicted number of massive galaxies also matches well with observations in this redshift range. However, the massive galaxies are predicted to show higher levels of residual star formation activity than the observational estimates, suggesting the need for further suppression of star formation without significantly altering the stellar mass function. We discuss possible improvements, involving faster stellar assembly through enhanced star formation during galaxy mergers while star formation at the peak epoch is still modulated by the AGN feedback.

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

State-of-the-art hydrodynamical simulations show that gas inflow through the virial sphere of dark matter halos is focused (i.e. has a preferred inflow direction), consistent (i.e. its orientation is steady in time) and amplified (i.e. the amplitude of its advected specific angular momentum increases with time). We explain this to be a consequence of the dynamics of the cosmic web within the neighbourhood of the halo, which produces steady, angular momentum rich, filamentary inflow of cold gas. On large scales, the dynamics within neighbouring patches drives matter out of the surrounding voids, into walls and filaments before it finally gets accreted onto virialised dark matter halos. As these walls/filaments constitute the boundaries of asymmetric voids, they acquire a net transverse motion, which explains the angular momentum rich nature of the later infall which comes from further away. We conjecture that this large-scale driven consistency explains why cold flows are so efficient at building up high redshift thin discs from the inside out.

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive mesh refinement, we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r=0.1rvir. In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/rvir>0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons remains close to that of dark matter halos. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

[Abridged] Cold gas flowing within the "cosmic web" is believed to be an important source of fuel for star formation at high redshift. However, the presence of such filamentary gas has never been observationally confirmed. In this work, we investigate in detail whether such cold gas is detectable using low-ionisation metal absorption lines, such as CII \lambda1334 as this technique has a proven observational record for detecting gaseous structures. Using a large statistical sample of galaxies from the Mare Nostrum N-body+AMR cosmological simulation, we find that the typical covering fraction of the dense, cold gas in 10^12 Msun haloes at z~2.5 is lower than expected (~5%). In addition, the absorption signal by the interstellar medium of the galaxy itself turns out to be so deep and so broad in velocity space that it completely drowns that of the filamentary gas. A detectable signal might be obtained from a cold filament exactly aligned with the line of sight, but this configuration is so unlikely that it would require surveying an overwhelmingly large number of candidate galaxies to tease it out. Finally, the predicted metallicity of the cold gas in filaments is extremely low (\leq 0.001 Zsun). Should this result persist when higher resolution runs are performed, it would significantly increase the difficulty of detecting filamentary gas inflows using metal lines. However, even if we assume that filaments are enriched to Zsun, the absorption signal that we compute is still weak. We are therefore led to conclude that it is extremely difficult to observationally prove or disprove the presence of cold filaments as the favorite accretion mode of galaxies using low-ionisation metal absorption lines. The Ly-alpha emission route looks more promising but due to the resonant nature of the line, radiative transfer simulations are required to fully characterize the observed signal.

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.

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation.

Using Sloan Digital Sky Survey Data Release 5, we have investigated the intrinsic axis ratio distribution (ARD) for early-type galaxies. We have constructed a volume-limited sample of 3,922 visually-inspected early-type galaxies at $0.05 \leq z \leq 0.06$ carefully considering sampling biases caused by the galaxy isophotal size and luminosity. We attempt to de-project the observed ARD into three-dimensional types (oblate, prolate, and triaxial), which are classified in terms of triaxiality. We confirm that no linear combination of $randomly$-distributed axis ratios of the three types can reproduce the observed ARD. However, using Gaussian intrinsic distributions, we have found reasonable fits to the data with preferred mean axis ratios for oblate, prolate, and triaxial (triaxials in two axis ratios), $\mu_o=0.44, \mu_p=0.72, \mu_{t,\beta}=0.92, \mu_{t,\gamma}=0.78$ where the fractions of oblate, prolate and triaxial types are \textrm{O:P:T}=0.29^{\pm0.09}:0.26^{\pm0.11}:0.45^{\pm0.13}$. We have also found that the luminous sample ($-23.3 < M_r \leq -21.2$) tends to have more triaxials than the less luminous ($-21.2 < M_r <-19.3$) sample does. Oblate is relatively more abundant among the less luminous galaxies. Interestingly, the preferences of axis ratios for triaxial types in the two luminosity classes are remarkably similar. We have not found any significant influence of the local galaxy number density on ARD. We show that the results can be seriously affected by the details in the data selection and type classification scheme. Caveats and implications on galaxy formation are discussed.