We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that
of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuation, about 2 mag in the B-band. Assuming the empirical bolometric corrections for active galactic nuclei (AGNs) by Marconi et al., we find that this dust attenuation is too strong to explain the B-band and X-ray quasar luminosity functions simultaneously. We consider two possible mechanisms that weaken the dust attenuation. As such a mechanism, we introduce a time delay for AGN activity, that is, gas fueling to a central black hole starts some time after the beginning of the starburst induced by a major merger. The other is the anisotropy in the dust distribution. We find that in order to make the dust attenuation of the quasars negligible, either the gas accretion into the black holes has to be delayed at least three times the dynamical timescale of their host bulges or the dust covering factor is as small as 0.1.
We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed
by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.
In this paper, we present a new implementation of feedback due to active galactic nuclei (AGN) in cosmological simulations of galaxy formation. We assume that a fraction of jet energy, which is generated by an AGN, is transferred to the surrounding g
as as thermal energy. Combining a theoretical model of mass accretion onto black holes with a multiphase description of star-forming gas, we self-consistently follow evolution of both galaxies and their central black holes. The novelty in our model is that we consider two distinct accretion modes: standard radiatively efficient thin accretion disks and radiatively inefficient accretion flows which we will generically refer to as RIAFs; motivated by theoretical models for jet production in accretion disks, we assume that only the RIAF is responsible for the AGN feedback. We find that, after an initial episode of bursting star formation, the accretion rate onto the central black hole drops so that the accretion disk switches to a RIAF structure. At this point, the feedback from the AGN becomes efficient and slightly suppresses star formation in the galactic disk and almost completely halts star formation in the bulge. As a result, the nucleus becomes a stochastically fuelled low-luminosity AGN (Seyfert galaxy) with recurrent short-lived episodes of activity after the star bursts. Our model predicts several properties of the low-luminosity AGN including the bolometric luminosity, jet powers, the effect on kpc-scale of the radio jet and the AGN lifetime, which are in broad agreement with observations of Seyfert galaxies and their radio activity. We also find that the mass ratios between the central black hole and the the host spheroid at z = 0 are ~10^{-3} regardless of the strength of either supernova feedback or AGN feedback. (abridged)
