We describe techniques for incorporating feedback from star formation and black hole (BH) accretion into simulations of isolated and merging galaxies.

We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak quasar luminosity, and we determine this dependence using a large set of galaxy merger simulations varying galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions and obtain the evolution of the formation rate of quasars with peak luminosity, (Lpeak, z). Quasars spend extended periods at luminosities well below peak, so (Lpeak) has a maximum corresponding to the break in the observed luminosity function. From (Lpeak) and our simulations, we obtain self-consistent hard and soft X-ray and optical luminosity functions and predict many observables at multiple redshifts, including column density distributions of optical and X-ray samples, the luminosity function of broad-line quasars in X-ray samples and broad-line fraction versus luminosity, active BH mass functions, the distribution of Eddington ratios, the mass function of relic BHs and total BH mass density, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, without invoking ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations.

We introduce a technique to filter out complex data sets by extracting a subgraph of representative links. Such a filtering can be tuned up to any desired level by controlling the genus of the resulting graph. We show that this technique is especially suitable for correlation-based graphs, giving filtered graphs that preserve the hierarchical organization of the minimum spanning tree but containing a larger amount of information in their internal structure. In particular in the case of planar filtered graphs (genus equal to 0), triangular loops and four-element cliques are formed. The application of this filtering procedure to 100 stocks in the U.S. equity markets shows that such loops and cliques have important and significant relationships with the market structure and properties.

We discuss a numerical model for black hole growth and its associated feedback processes that for the first time allows cosmological simulations of s